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ABSTRACT 


In  this  volume  is  presented  an  analysis  of  the  static  anu  dynamic  characteristics 
of  the  epical-grooved  journal  bearing  operating  with  incompressible  lubricant  in 
both  laminar  and  turbulent  regimes-  Both  single  film  and  floating  ring  bearing 
configuratlona  are  considered.  Extensive  design  data  are  presented  giving  load 
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two  computer  programs  accompany  the  volume,  and  instructions  and  listings  of  the 
programs  are  included.  These  programs  may  be  used  to  obtain  data  for  cases  not 
covered  by  the  presented  design  data 


This  abstract  is  subject  to  special  export  controls 
and  each  transmittal  to  foreign  governments  or  foreign 
nationals  may  be  made  only  with  prior  approval  of  the 
Air  Force  Aero  Propulsion  Laboratory  (APFL),  Wright- 
Patterson  Air  Force  Base,  Ohio  45433- 


ill 


TABLE  OF  CONTENTS 


"ant 

INTRODUCTION . I 

DESCRIPTION  OF  SINGLE  FILM  AND  FLOATING-RING 

SPIRAL  GROOVED  JOURNAL  BEARINGS . 3 

ANALYSIS . 7 

RESULTS . 27 

SUMMARY  AND  CONCLUSIONS  . . 39 

TABUS .  41 

ILLUSTRATIONS .  47 

APPENDIX  I:  Derivation  and  Solution  of  Equation  for  Turbulent 

Spiral-Grooved  Journal  Bearing  _ _ _ — _  89 

APPENDIX  II:  Derivation  of  Relationshipe  for  Calculating 

Stiffness  and  Damping  Coefficients  ............ — ........  101 

APPENDIX  III:  Computer  Program  PN  412 

Performance  of  a  Herringbone  Journal  Bearing 
Operated  in  the  Turbulent  Regime  ........................  107 

APPENDIX  IV:  Computer  Program  PN  406 

Static  Performance  of  a  Spiral- Grooved,  Floating 

Ring  Journal  Bearing  Operated  in  the  Turbulent  Regime  ...  131 

REFERENCES .  175 


ILLUSTRATIONS 


Figure 


Page 


1  Schematic  of  Spiral-Grooved  Journal  Bearing  _ 

2  Various  Configurations  of  Spiral-Grooved  Journal  Bearings  - 

3  Coordinate  System  for  Forces  and  Displacements  _ 

4  Geometry  of  Floating  Ring  Bearing  _ 

5  G  and  G  vs.  Local  Reynolds  Number  R,  _ 

x  z  h 

6  Coefficient  of  Friction  vs-  Local  Reynolds  Number  R,  _ 

i  h 

7  Load  vs.  Eccentricity,  Single  Film  Bearing  ...... _ .......... _ 

8  ;  Load  vs.  Eccentricity,  Single  Film  Bearing  _ 

9  Load  v 8 .  Eccentricity,  Single  Film  Bearing - - - - - 

10  Load  vs.  Eccentricity,  Effect  of  Pressurization  _ 

11  Attitude  Angle  va.  Eccentricity,  Single  Film  Bearing  .... - 

12  Attitude  Angle  va.  Eccentricity,  Single  Film  Bearing  - 

13  Attitude  Angle  vs.  Eccentricity,  Single  Film  Bearing  _ .... 

14  Single  Film  Bearing  Flow  Due  to  Self  Pumping  Effect 

of  Spiral  Grooves  ... _ ............... 

15  Single  Film  Bearing  Flow  Due  to  Pressurized  Supply _ 

16  Torque  vs.  Eccentricity,  Single  Film  Bearing  - 

17  Torque  vs.  Eccentricity,  Single  Film  Bearing  ---- _ - _ 

18  Torque  to  Load  Ratio  vs.  Reynolds  Number,  Single  Film  Bearing  _ 

19  Load  vs.  Inner  Film  Eccentricity  Ratio,  Floating  Ring  Bearing  .... 

20  Attitude  Angle  vs.  Inner  Film  Eccentricity  Ratio,  Floating  Ring 

Bearing - - - - .... _ 

21  Attitude  Angle  vs.  Inner  Film  Eccentricity  Ratio,  Floating  Ring 

Bearing  - 

22  Torque  to  Load  Ratio  vs.  Inner  Film  Eccentricity,  Floating  Ring 

Bearing  - 


47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 


67 


68 


vii 


Illustrations  (Continued) 


F  1-KU1  t:  raxc 

23  Torque  to  Load  Ratio,  Comparison  between  Single  Film  and 

Floating  Ring  Bearing - 69 

24  Floating  Ring  Bearing  Flow  Due  to  Self-Pumping  Effect  of  Spiral 

Grooves _  70 

23  Ring  Speed  Ratio  vs.  Inner  Film  Eccentricity  Ratio,  Floating 

Ring  Bearing -  71 

26  Outer  Film  tccenrricity  Ratio  vs.  Inner  Film  Eccentricity  Ratio, 

Floating  Ring  Bearing _ - _  72 

27.  Stiffness  Coefficients,  Floating  Ring  Bearing  -  73 

28  Damping  Coefficients,  r loa ting  Ring  Bearing  ......... - - .....  74 

29  Stiffness  Coefficients,  Floating  Ring  Bearing  ..... - ..........  73 

30  Damping  Coefficients,  Floating  Ring  Bearing - -  76 

31  Stiffness  Coefficients,  Floating  Ring  Bearing  - - 77 

32  Damping  Coefficients,  Floating  Ring  Bearing - .... -  78 

33  Critical  Journal  Mass,  Floating  Ring  Bearing  _ 79 

34  Critical  Journal  Mass,  Floating  Ring  Bearing - ...... - ...  80 

33  Critical  Journal  Mass,  Comparison  between  Plain  and  Spiral- 

Grooved  Floating  Ring  Bearings  - 81 

36  Critical  Journal  Mass,  Comparison  between  Plain  and  Spiral- 

Grooved  Floating  Ring  Bearings  _  82 

37  Critical  Journal  Mass,  Comparison  between  Single  Film  and 

Floating  Ring  Spiral-Grooved  Bearings _ _  83 

38  7]  and  £  Coordinates  for  Grooved  Surface _  84 

35  Pressure  Distribution  around  Spiral-Grooved  Journal  Bearing  _  85 

40  Continuity  of  Mass  Flow  across  Groove-Ridge  Interface  and 

Pressure  Variation  across  Groove-Ridge  Interface - ............  86 

41  Control  Volume  for  Mass  Flow  Continuity  Analysis  _ 87 


viii 


SYMBOLS 


B  ,  B  , 
xx  xy 


B  ,B 
yx’  yy 


Bearing  damping  coefficients,  Lb-sec/in. 


B  ,B 

xx  xy 

B  ,B 
yx  yy 


Dimensionless  damping  for  overall  floating  ring  bearing, 

d  n  n 


XX 


vs  ,s.  t 
<*?  ct“- 


<■-->  *<v 


xx 


(B  )  ,  (B  ) 

UV  *  ' 


yx 


yy 


Dimensionless  damping  for  inner  film, 
B  ,C.  C.  2 

2*  etc 

HLD  '  etCl 


(B  ) 
xx  , 


(B  )  ,(B  ) 

xx  2,x  xy  2 

(B  )  ,(B  ) 

•  r  '  '  IMt  r 


yx 


yy 


Dimensionless  damping  for  outer  film, 


(B  ) 
xx  , 


B  C  C 

fa?  .?  etc 

tiLD,  CC 


C 

D 


e 

F_ 


F  ,F 

x  y 


F 


t 


6  ,  G 
x’  z 


Bearing  mean  radial  clearance,  in. 

Bearing  diameter,  in. 

Eccentricity,  in. 

Radial  component  of  bearing  film  force  -  WcosO,  lb. 

Tangential  component  of  bearing  film  force  *  Wsinrf,  lb. 

Bearing  film  forces  in  x  and  y  directions,  lb. 

Dimensionless  radial  (cosine)  component  of  bearing  film 

force  ■  Wcos»i  C^/n(N,  +  N  )R^ 

i  o 

Dimensionless  tangential  (sin)  component  of  bearing  film 

force  *  Wsin(S  C^/n(N.  +  N  )R^ 

l  o 

Turbulent  flow  correction  factors 


ix 


SYMBOLS  Crnn  t~  i  nnpH  1 


h 

h 

g 

h 

r 

h' 


K  ,K  , 
xx  xy 

K  ,K 
yx’  yy 


Local  film  clearance,  In. 

Local  film  clearance  in  groove  region,  in. 
Local  film  clearance  in  ridge  region,  in. 
Dimensionless  clearance  =  h/C 
Bearing  stiffness  coefficients,  lb/in. 


K  ,K  , 
xx  xy 

K  ,K 
yx  yy 


Dimensionless  stiffness  for  overall  floating  fing  bearing, 


xx 


KXX  C1  ^ 

“W  V 


etc  ■ 


(U  » <K  )  , 


xx 


1 


xy 


(K  )  ,(K  ) 

in*  ' 


yx 


yy 


1 


Dimensionless  stiffnesses  for  inner  film, 
K  ,  C,  C  2 

(K  )  *>  ■  (— -L\  etc 

Q  xx;j  n(Nx  +  N2)LD1  ^  etc- 


(K  )  , (K  )  ,  Dimensionless  stiffness  for  outer  film, 


xx 


_xy 

(K  )  ,(K  .) 


yx 


yy 


K  „  C„  C„  * 

(K  )  ■  — SS2 — —  (—2.)  e  , 
'»'2  nN2  ld2  vr2; 


L 

L. 

M 

M 


c 


cl 


M 


c2 


Bearing  length,  in. 

Length  of  grooved  region,  in. 

Shaft  mass,  lb. 

Critical  mass  for  threshold  of  whirl  instability  in  floating 
ring  bearing,  lb. 

Inner  film  critical  mass. 

Outer  film  critical  mass. 


x 


SYMBOLS  (continued) 


Radial  component  of  restoring  moment,  in-lb. 

Tangential  component  of  restoring  moment,  In- lb. 

Dimensionless  critical  mass  = 

Dimensionless  critical  mass  - 

Dimensionless  critical  mass  - 

"  enslonless  radial  (cosine)  component  of  bearing  film 
force  moment  »  Mr  C^/p(N^  +  Nq)R5 

Dimensionless  tangential  (sin)  component  of  bearing  film 
•>  5 

force  moment  ■  M.  C*7p(N.  +  N  )R 
t  l  o 

Ring  mass,  lb. 

Speed  «  +  N2>  rps 

Speed  of  journal,  rps 

Ring  speed,  rps 

Speed  of  inner  member,  rps 

Speed  of  outer  member,  rps 

Ring  speed  ratio  -  N2^N1 

2 

Pressure,  lb/in 

o 

Supply  pressure,  lb/in 


M  N,C. 
c  1  1 

^r12ld1 


MciO*i  +  yy 
uRt2  LDt 


Mc2M2C2 
“» 2 2  “2 


SYMBOLS  {continued) 


bimehsionless  pressure  «  Hlj^/pN^Cj2 

Smoothed,  "overall"  pressure  distribution  around  spiral 

grooved  journal,  lb/ in2 

Dimensionless  supply  pressure  «  P  R  2/liN  C  2 

8  1  1  I 

Dimensionless  supply  pressure,  inner  film  -PR  2/m(n  +  Nlr 

s  l  1  T  1 

Dimensionless  supply  pressure,  outer  film  -PR  2/un  c  2 

s  2  2  2 

Total  lubricant  flow  rate,  in3/sec. 


Lubricant  flow  due  to  self  pumping  of  spiral  grooves,  in3/sec. 
Lubricant  flow  due  to  pressurization  of  supply,  in3/sec. 


Dimensionless  total  lubricant  flow  rate 
Bearing  radius,  in. 

Overall  Reynolds  number  -  2*  NjRjC^/v 

Local  Reynolds  number  -  2 «  (N.  -  N  )Rh/v 

1  o 

Inner  film  Reynolds  Number 
Outer  film  Reynolds  Number 


Q/k  C(N  +  N  ) 
l  o 


2*P(Nq  -  N2)RlCl/n 

2jcP  n2R2c2/m 


uN.d.l  R,  * 

Overall  Sommerfeld  Number  -  — = — —  I _ =■'. 

W 


Inner  film  Sommerfeld  Number  * 


M(NL  +  N2)DjL 

w 


R1 

% 


2 


Outer  film  Sommerfeld  Number 


W 


Bearing  torque,  in- lb. 


SYMBOLS  (continued) 


t 

U 


V 

w 

w*.  ^ 

v 

u,y 

Y 


z  ,z  , 

xx  xy 

z  ,z 

yx’  yy 


Z' 


Journal  Torque,  in- lb- 
Dimensionless  bearing  torque  «*  T^/UC 
Dimensionless  Journal  torque  -  T^/WC 
Time,  sec. 

Surface  velocity  of  bearing,  in/scc . 

Mean  flow  velocity  in  direction  of  rotation,  in/sec. 

Mean  flow  velocity  due  to  pressure  gradient  ■  u  -  V/2,  in/sec. 
Surface  speed  of  journal,  ln/sec. 

Bearing  load,  lb. 

Mass  flow  components,  lb/ (sec- in  ) 

Mean  flow  velocity  in  axial  direction,  ln/sec. 

Coordinates  in  direction  of  and  normal  to  load  vector,  in. 
Ratio  of  length  of  grooving  to  total  length  of  bearing  “  L^/L 
Complex  dimensionless  bearing  lmpedences 


Dimensionless  Z  coordinate  ■  Z/L 


Greek 

a 

e 

r 

Y 


Ratio  of  groove  width  to  groove  plus  land  width  m  a^/(a^  +  af) 
Groove  angle,  deg.  (radians  in  equations) 

Ratio  of  groove  clearance  to  ridge  clearance  ■  h^/h^. 
Misalignment  angle,  deg.  (radians  in  equations) 


xiii 


SYMBOLS  (continued) 


Y 

e 


n.  « 

9,  Z 

\ 

4 

v 

v 

P 

i 

6' 


<u 


Whirl  frequency  ratio  m  v/u> 

Eccentricity  ratio  -  e/C 

Skewed  coordinates 

Cylindrical  coordinates 

Dimensionless  parameter  ■  (c/R)2 

* 

Viscosity,  lb-sec/in2 
Kinematic  viscosity,  in  /sec. 

Whirl  frequency,  radians/sec. 

Density,  lb/in3 

Attitude  angle,  deg.  (radians  in  equations) 

Moment  attitude  angle,  deg.  (radians  in  equations) 
Total  rotational  speed  -  (o^  +•  ci^) ,  radians/sec. 
Rotational  speed  Journel,  radians/sec. 

Rotational  speed  ring,  radians/sec. 


Subocrlpta 

1 

2 


Refers  to  inner  film 
Refers  to  outer  film 


xiv 


INTRODUCTION 


In  many  applications  of  rotating  machinery  it  is  desirable  to  lubricate  bearings 
with  the  process  fluid  in  order  to  avoid  the  complication  of  a  separate  lube 
system  with  accompanying  seal  problems.  In  many  instances,  process  fluids  may 
have  low  kinematic  viscosities  which  result  in  operation  of  the  bearings  in  the 
turbulent  flow  regime. 

Bearing  power  loss  rises  rapidly  with  Reynolds  number  in  turbulent  film  bearings 
In  a  number  of  prototype  systems  Involving  rotary  machines  operating  In  the  tur¬ 
bulent  regime,  the  bearing  power  lose  has  bean  sn  appreciable  percentage  of  the 
net  eyetem  output.  In  addition  to  the  effect  on  eyetem  efficiency,  high  bearing 
power  losaes  mean  that  large  quantities  of  lubricant  must  be  circulated  through 
the  bearing  for  cooling. 

A  very  attractive  bearing  for  use  in  applications  where  power  loae  la  important 
is  the  floating  ring  bearing.  This  la  a  journal  bearing  in  which  a  loose  ring 
la  fitted  between  the  aheft  end  the  bearing  housing.  Thle  ring  is  free  to  rotate 
when  the  journal  rotates,  end  by  eo  doing,  can  reduce  the  rate  of  ahsar  between 
adjacent  bearing  surfaces  thereby  reducing  power  loei. 

In  plain  bearing  form,  the  principal  disadvantage  of  the  floating  ring  bearing  la 
the  rather  poor  etabillty  characteristics  of  plain  journal  bearings,  particularly 
when  lightly  loaded.  There  arc  a  variety  of  configurations  of  Journal  bearings 
which  offer  improved  stability  over  plain  journal  bearings.  These  Include  the 
tilting  pad  Journal  bearing,  the  Rayleigh  step  journal  bearing,  multi-lobed  journ 
al  bearings,  and  the  spiral-grooved  Journal  bearing.  Of  these,  the  most  suitable 
for  use  in  a  floating  ring  configuration  ia  the  spiral-grooved  bearing.  In  ad¬ 
dition  to  improved  stability  characteristics  over  plain  bearings,  the  spiral- 
grooved  bearing  also  has  several  other  attractive  advantages.  These  Include  the 
ability  to  self-pump  lubricant  axially  through  the  bearing  film  and  the  ability 
to  operate  without  cavitation  at  high  eccentricity  ratios  without  using  a  pres¬ 
surized  lubricant  supply. 
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Because  ol  the  potential  advantages  of  spiral -grooved  floating  ring  journal  bear¬ 
ings  for  many  applications,  an  analysis  vaa  performed  of  this  bearing  for  both 
laminar  and  turbulent  floe  regimes.  This  analysis  and  the  results  obtained  there¬ 
of  form  the  subject  of  this  report.  Extensive  performance  data  are  presented  for 
both  floating  ring  and  single  film  configuration  of  the  spiral-grooved  bearing. 
These  data  Include  stiffness  and  damping  coefficients  and  evaluations  of  critical 
mass  for  fractional  frequency  whirl  instability. 
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DESCRIPTION  OF  SINGLE  FILM  AND  FLOATING-RING  SPIRAL  GROOVED  JOURNAL  BEARINGS 


A  typical  apiral-grooved  journal  bearing  is  shown  schematically  in  Fig.  1.  The 
configuration  shown  is  that  In  which  the  grooving  on  the  ends  of  the  journal 
tends  to  pump  Inward,  thereby  pressurizing  the  interior  o£  the  bearing.  Because 
of  the  symmetry  of  this  configuration,  there  is  no  net  flow  of  lubricant  through 
the  bearing . 

Other  possible  configurations  for  the  single  film  spiral-grooved  bearing  are  shown 
in  Fig.  2.  Configuration  1  is  as  described  above.  Configuration  2  is  an  arrange¬ 
ment  wherein  spiral  grooving  is  inscribed  on  one  end  of  the  bearing  in  a  manner 
so  as  to  pump  fluid  axially  through  the  bearing  toward  the  smooth  or  seal  end. 

The  axial  flow  of  lubricant  can  remove  heat  from  the  bearing.  The  flow  through 
the  bearing  can  be  increased  by  pressurizing  the  grooved  end  of  the  bearing  to 
a  supply  pressure  Ps  as  shown. 

Configuration  3  shows  a  symmetrical  arrangement  wherein  lubricant  is  pumped  by 
spiral  grooves  outward  from  the  center  of  the  bearing.  Lubricant  may  be  supplied 
to  the  center  of  the  bearing  at  an  elevated  pressure  PB .  Configuration  3  con¬ 
sists,  essentially,  of  two  bearings  of  configuration  2  placed  "back  to  back". 

In  all  the  configurations  illustrated  in  Pigs.  1  and  2,  the  spiral-grooving  1b 
shown  Inscribed  on  the  journal  which  is  considered  to  be  the  rotating  member. 

One  could, Instead,  inscribe  the  grooving  on  the  stationary  member  and,  for  an 
incompressible  lubricant,  the  performance  of  the  bearing  would  be  essentially 
the  same.  In  this  report,  we  shall  only  consider  the  case  where  the  grooves 
are  Inscribed  on  the  rotating  member,  this  being  the  most  common  situation  in 
practice . 

The  performance  of  a  spiral- grooved  bearing  depends  on  the  values  of  the  groove 
parameters  3  (groove  angle),  a  (ratio  of  groove  width  to  total  width), P  (ratio 
of  groove  clearance  to  ridge  clearance), and  Y  (ratio  of  length  of  grooving  to 
total  length  of  bearing).  These  groove  psrameters  are  defined  in  Fig.  1.  By 
a  suitable  choice  of  these  parameters,  one  can  optimize  various  performance 
characteristics  of  the  bearing.  In  this  report,  we  have  chosen  to  present  per¬ 
formance  data  for  the  case  where  the  groove  parameters  are  set  at  those  values 
which  yield  maximum  radial  component  of  bearing  stiffness  at  zero  eccentricity. 
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The  radial  component  of  bearing  stiffness  at  aero  eccentricity  is  defined  as 
the  limiting  value  of  W  cos  d/e  as  e  approaches  ?ero  where  Vf  is  bearing  load, 
j.6  the  journal  displacement  from  the  center  of  the  bearing  and  d  is  the  at> 
tlj-uhe  angle  for  the  bearing  (see  fig.  3).  In  determining  the  optimum  values 
of  groove  parameters  for  maximum  radial  stiffness.  It  turned  out  that  optimum 
values  of  a,  9  and  Y  did  not  vary  appreciably  with  Reynolds  number  or  with  L/D 
ratio  in  the  range  0.5  <  L/D  <  1.0.  For  simplicity,  therefore,  a  fixed  set 
of  optimum  values  for  these  parameters  was  settled  upon  as  valid  for  all  Reynolds 
numbers  and  all  L/D  ratios  between  0.5  and  1.0.  The  optimum  value  of  groove 
depth  ratio  T  however,  did  vary  significantly  with  Reynolds  number.  Optimum 
values  of  a,  3,  Y  and  F  selected  for  configuration  1  and  configuration  2  type 
bearings  are  given  in  Tables  1  and  2  below. 


TABLE  1 

OPTIMUM  VALUES  OF  GROOVE  PARAMETERS  SELECTED  FOR  CONFIGURATION  1  WAR  INS 
Reynolds  No.  a  3  Y  F 


laminar,  Re  <  500 

0.5 

151.5* 

0.75 

2.1 

1,000 

ii 

•1 

ii 

2.7 

5,000 

ti 

II 

ii 

3.0 

9,000 

ii 

II 

ii 

3.0 

TABU!  2 

-PTIMUM  VALUES  OF  GROOVE  PARAMETERS  SELECTED  TOR  CONFIGURATION  2  BEARING 
Reynolds  No.  a  &  Y  F 


laminar,  Re  <  500  0.55 

149° 

0.67 

2.4 

1,000  " 

ii 

ii 

3.1 

5,000  " 

ii 

ii 

3.8 

1,000 

5,000 

9,000 


•i 
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Fig.  4  shows  a  schematic  drswing  of  a  floating-ring  bearing.  A  floating-ring 
bearing  is  a  journal  bearing  in  which  there  is  a  loose  ring  between  the  shaft 
snd  the  bearing  housing.  In  this  way,  the  fluid  film  is  separated  into  an 
inner  film  and  an  outer  film.  The  quantities  connected  with  the  inner  film 
are  identified  by  subscript  1,  whereas  subscript  2  refers  to  the  outer  film. 

The  floating-ring  bearing  shown  in  Fig.  4  is  shown  with  spiral- grooving  on 
the  journal  and  on  the  outer  surface  of  the  floating  ring.  The  configuration 
of  the  grooving  is  such  as  to  pump  lubricant  outward  from  the  axial  midplane 
of  the  bearing.  Lubricant  is  supplied  to  the  two  lubricant  films  at  supply 
pressure  via  supply  holes  in  the  bearing  and  in  the  floating  ring.  Circum¬ 
ferential  grooves  machined  at  the  midplane  of  the  Journal  and  the  floating 
ring  distribute  the  lubricant  at  uniform  pressure  around  the  journal  and 
around  the  outeide  of  the  floating  ring.  The  ring  is  free  to  rotate,  and 
under  the  influence  of  shear  stress  from  the  revolving  journal,  turns  at  some 
ring  speed  H2  less  than  the  speed  of  the  journal.  Since  in  journal  bearings, 
load  capacity  is  essentially  proportional  to  +  N^)  whereas  power  loss  is 
roughly  proportional  to  (N^  -  N^) ,  it  follows  that  rotation  of  the  ring  will 
improve  the  load  capacity  of  the  inner  film  while  reducing  the  power  loss.  This 
is  the  principle  of  operation  of  the  floating  ring  bearing. 

Calculation  of  the  performance  characteristics  of  the  floating  ring  bearing 
requires  calculation  and  matching  of  the  performance  characteristics  of  the 
individual  lubricant  films.  These  individual  performance  characteristics 
depend  on  the  values  selected  for  the  groove  parameters,  with  the  possibility 
of  having  different  groove  parameters  for  the  inner  and  outer  films.  For  the 
calculations  presented  in  this  report,  the  groove  parameters  were  taken  to  be 
the  same  for  the  inner  and  outer  film  and  were  selected  to  be  those  which  pro¬ 
vide  maximum  radial  component  of  bearing  stiffness  for  each  individual  bearing 
film.  (See  Table  2 ) . 

A  description  of  the  analysis  of  the  turbulent  single  film  and  floating  ring 
spiral- grooved  journal  bearings  is  given  in  the  next  section. 
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ANALYSIS 


Analysis  of  the  performance  characteristics  of  the  turbulent,  spiral-grooved, 
single  film  and  floating-ring  journal  bearing  is  based  on  the  concept  of  solving 
for  the  "overall",  "smoothed"  pressure  distribution  around  the  bearing,  neglect¬ 
ing  the  local  sig-eag  details  of  the  pressure  profiles  which  arise  due  to  the 
discontinuous  groove-ridge  geometry,  "he  theoretical  basis  for  this  analytical 
approach  is  discussed  in  detail  in  Reference  1.  Essentially,  this  analytical 
approach  is  valid  in  the  limit  as  the  number  of  grooves  approaches  an  infinite 
number,  but  practically  speaking,  the  analysis  proves  to  be  quite  accurate  when 
applied  to  bearings  having  a  reasonable  number  of  grooves.  Experimental  verifi¬ 
cation  of  this  analytical  approach  has  been  provided  by  a  number  of  investigations 
(References  2  and  3 ) . 

The  differential  equation  that  had  been  derived  in  Reference  1  for  the  smoothed, 
overall  pressure  distribution  around  a  spiral-grooved  journal  bearing  was  re- 
derived  in  the  present  study  to  take  account  of  the  effects  of  turbulence  in 
the  bearing  film.  This  derivation  is  presented  in  Appendix  1.  The  effects  of 
turbulence  in  the  bearing  film  are  accounted  for  by  means  of  the  linearised 
turbulent  lubrication  theory  developed  by  Mg  and  Fan  1  (Reference  4).  In  thia 
theory,  which  is  based  on  the  concept  of  a  turbulent  eddy  vlacosity,  there  ere 
developed  turbulent  flow  correction  factors  6  and  G  which  relate  the  mean 
pressure  flow  In  the  direction  of  rotation  (x  direction)  and  the  axial  direction 
(t  direction)  to  the  pressure  gradients  in  these  respective  directions.  The 
relationships  developed  are 
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where 


—  —  V 

u  ■  u  -  j  ■  the  mean  flow  velocity  in  the  x  direction  minus  1/2  the 
surface  velocity  of  rotation 

w  •  the  mean  flow  velocity  in  the  axial  direction 

,h  ■  local  film  clearance 

4  “  viscosity 

In  the  genera Used  theory  for  turbulent  fluid  films  developed  by  Elrod,  Mg  and 
fan  (Reference  5),  G  and  C  are  functions  of  the  pressure  gradient  in  the  film, 
the  angle  between  the  pressure  gradient  snd  the  direction  of  rotation,  and  the 
Reynolds  number  based  on  rotational  velocity.  In  the  linearised  theory  of  Ng 

and  Fan,  C„  and  G.  are  functions  only  of  the  local  Reynolds  number  R.  ■  p  V  h/p. 

x  *  n 

Values  of  C  and  G  ,  plotted  as  a  function  of  R.  ,  are  shown  In  Figure  5. 
x  s  n 

A  discussion  of  the  theoretical  basis  of  the  linearised  theory  of  turbulence  is 
beyond  the  scope  of  this  present  report.  For  such  a  discussion,  the  reader  can 
consult  Reference  4.  In  this  report,  we  have  simply  applied  the  results  of  this 
theory  to  derive  the  differential  equation  for  a  apiral- grooved  journal  bearing 
with  turbulent,  incompressible  lubricant.  This  differential  equation,  obtained 
in  Appendix  1,  is  given  below  in  dimensionless  form. 
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A1  ”  Al/A3 
Aj  "  A^ /A^ 

®1  *  Bl/A3 

b2  -  b2/a3 

and  where  ,  G^,  Q^r>  G^j  etc.  are  lumped,  turbulent  flow  correction  factora 
defined  by  Eqa.  (65) ,(66)  and  (67)  in  Appendix  X. 

Eq.  (3)  vaa  aolved  numerically  on  a  digital  computer  uaing  the  method  of  column* 
viae  Influence  coefflcienta  developed  by  Caatelll  and  Shapiro  (Ref.  6)  and 
Caatelli  and  Pirvic-  (Ref.  7).  Two  separate  computer  programs  were  developed, 
one  to  obtain  reaulta  for  the  static  and  dynamic  characteristics  of  a  single 
film  spiral-grooved  bearing,  and  the  second  to  calculate  the  overall  static  per¬ 
formance  characteristics  of  a  spiral -grooved  floating-ring  bearing. 

Calculation  of  Performance  of  Single  Film.  Spiral-Grooved  Journal  Bearing 

The  computer  program  for  calculating  the  performance  characteristics  of  a  single- 
film  spiral -grooved  journal  bearing  requires  that  the  following  quantities  be 
specified.  (Symbols  are  defined  in  the  nomenclature). 

Reynolds  number  based  on  mean  radial  clearance  in  the  seal  region  of  the 


bearing  " 


2«p(N  -  N  )RC 
l  o 


L/D  ratio 
C/R  ratio 
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ni5“S2ier.lcs5  pr cuan  ai  uui.ii  ends  or  cne  Dearing 

Nt>  '  Ni 

Speed  ratio  factor  *  — +n 

"i  o 

Dimenaionleas  rate  of  change  of  eccentricity  ratio 

Dimensionless  whirl  speed  ratio  ■  ^/2rt(R^  +  N^) 

Dimensionless  rate  of  change  of  miaalignment  angle 

Eccentricity  ratio,  e 

Angle  of  misalignment,  J 

Groove  geometry  parameters  a,  &,  I\  Y 


p _ 

+  V 


2 


dt/2n<Ni  +  No> 


&L 

dt 


/2*(»i  +  No) 


In  addition,  one  must  specify  whether  the  bearing  Is  of  configuration  1  or 
configuration  2  (see  Fig.  2).  Subject  to  the  above  input  conditions,  Eq.  (3) 
is  solved  by  the  computer  program  to  determine  the  dimensionless  pressure 
distribution  P.  From  this  pressure  distribution,  the  computer  program  then 

I 

determines  the  following  performance  characteristics  of  the  bearing. 


Dimensionless  tangential  (sin)  component  of  bearing  fils  force, 
-  '  W  tin  t  /C,2 


P 


n(Nt  +  No)R' 


2  <f> 


(aee  Fig.  3) 


Dimensionless  radial  (cos)  component  of  bearing  film  force, 
?.  -  4>‘  «...  Fig.  3) 


n(Nt  +  No)E 

Dimensionless  tangential  (sin)  component  of  the  moment  exerted  by  the 
bearing  film  force  about  an  axis  through  the  initial  end  of  the  bearing 
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The  problem  of  cevltetioa  of  Che  bearing  film  la  handled  by  the  approximate 
method  of  netting  all  sub-ambient  fluid  film  preaeure  equal  to  sero  before 
integrating  for  loads  and  flows.  Experience  with  plain  journal  bearings  in* 
dlcates  that  thla  approach  yields  values  for  load  which  are  on  the  order  of 
51  to  10%  conservative  when  compared  to  a  more  exact  treatment  (Ref.  8).  For 
spiral-grooved  bearings,  the  extent  of  cavitation  is  much  less  than  for  plain 
bearings.  Thus,  one  would  expect  that  the  error  introduced  by  the  approximate 
method  of  handling  cavitation  would  not  be  significant  in  the  case  of  spiral- 
grooved  bearings. 

In  the  calculation  of  bearing  torque,  it  is  assumed  that  regions  of  sub- 
ambient  pressure  are  cavltated  and  therefore  do  not  contribute  to  the  shear 
stress  on  the  journal  or  bearing. 

» 

The  program  for  the  eingle  film  bearing  calculates  values  of  the  radial  end 
tangential  components  of  fluid  film  force  Fr  and  Ffc  as  functions  of  the  steady 
state  eccentricity  of  the  Journal,  e,  the  instantaneous  rate  of  change  of 
eccentricity  of  the  Journal,  de/dt,  and  the  instantaneous  whirl  velocity  of 
the  Journal  dd/dt.  Let  us  now  see  how  this  program  nay  be  used  to  obtain 
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atiffneaa  and  damping  coefficienta  for  the  bearing.  Conaider  the  reference  axea 
x  and  y  ahown  in  Fig.  3  where  x  ia  taken  aa  the  direction  of  the  ateady  at*.*  load 
V,  and  y  ia  normal  to  thia.  The  atiffneaa  end  damning  coefficients  arc  defined  by 


dF  •  -  K  x  -  B  -  K  v-B 

x  xx  xx  dt  xy  xy 


(4) 


dF  »  -  K  x-B  ^7  -  It  y-B  ^ 
y  yx  yx  it  yy  yy  ot 


(5) 


It  ia  ahovn  in  Appendix  II  that  for  bearinga  poaaeaaing  rotational  aymmetry, 
etc.  may  be  determined  from  derivativea  of  Ff  and  F£  with  reapect  to  e  and 
e  by  meana  of  the  following  expreaaiona 
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where  *  »  ^  t*  the  instantaneous  rate  of  change  of  eccentricity  and  i> 
ia  the  whirl  velocity  of  the  Journal. 
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The  determination  of  the  derivatives  dFr/de,  dF^/de  etc.  from  the  alngle  film 
program  la  tedloua  but  esaentlally  straight  forward.  Care  must  be  taken,  however, 
that  the  changes  in  e,  e  and  i  used  to  evaluate  theae  derivativea  be  chosen  suffi¬ 
ciently  small  such  thst  the  results  accurately  apply  to  Infinitesimally  amall 
amplitude  motions  of  the  journal  center  about  a  steady  state  position*.  It  is 
recommended  that  one  keep  Ae/C  <  .05, 4e/2x(Mi  +  No)C  <  .05  and  +■  1^)0 

<  .05  for  evaluation  of  the  above  mentioned  derivativea. 

A  detailed  description  of  how  to  prepare  input  for  the  alngle  film  spiral-grooved 
journal  bearing  program  is  given  in  Appendix  ZZI  together  with  a  listing  of  the 
program. 

Calculation  of  Steady  State  Performance  of  Floating  Ring  Bearing 

The  program  for  calculating  the  performance  of  a  floating  ring,  spiral-grooved 
bearing  consists,  essentially,  of  two  parta.  The  first  part  contains  the  program 
for  a  single  film,  spiral  grooved  journal  bearing  described  above.  This  la  used 
to  calculate  the  individual  performance  characteristics  of  the  inner  and  outer 
films  of  the  floating  ring  bearing.  The  second  part  of  the  program  consists  of 
the  logic  required  to  determine  the  correct  ring  speed  and  eccentricity  ratio  of 
the  outer  film  such  that  the  load  capacity  of  the  outer  film  is  equal  to  the  load 
capacity  of  the  inner  film  and  the  torque  exerted  by  fluid  shear  stresses  on  the 


it 

The  linear  formulation  represented  by  Eqa.  (4)  and  (5)  implicitly  carries  the 
assumption  that  the  motions  x,  y,  dx/dt,  dy/dt  are  vanishingly  small. 
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inner  surface  of  the  i«  balanced  by  the  fluid  snear  ■  tresses  on  the  outer 

surface  of  the  ring.  In  detail,  the  computational  procedure  works  ss  described 
below 

1.  To  determine  the  overall  preformance  of  the  floating  ring  bearing  requires 
that  the  following  quantities  be  specified  as  input  to  the  program. 

Radius  ratio  of  ring  and  Journal  *  R^/R^  where  R^  1b  the  outside  radius 
of  the  ring  and  R^  is  the  radius  of  the  journal. 

Overall  Reynolds  number  under  which  the  bearing  is  to  be  operated 
Re  -  RJ^RjCj^ 

Eccentricity  ratio  of  inner  film  to  be  examined  ■ 

Clearance  to  radiua  ratio  for  inner  film  ■  C^/R^ 

Clearance  to  radius  ratio  for  outer  film  »  C2^2 

2 

Dimensionless  supply  pressure  to  center  of  bearing  a(P(/dN^)  (C^/R^) 

Length  to  inner  diameter  ratio  for  bearing  ■  L/D^ 

Groove  geometry  parameters  a,  0,  T  and  Y 


2.  Given  this  input,  the  first  thing  the  program  does  is  to  calculate  and  store 
in  tabular  form  the  following  performance  data  for  the  inner  and  outer  films  of 
the  floating  ring  bearing. 

Inner  film 

Sommer f eld  Mo. 


S1  " 


U(NX  +  M2)D1L 


C1 


Dimensionleas  torque  on  Journal  (T,)  ■ 

3  1 
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Attitude  angle 

These  performance  date  ere  stored  in  tables  as  functions  of  the  eccentricity 
ratio  and  Reynolds  number  for  the  film  concerned.  Values  are  calculated  for 
three  predetermined  values  of  eccentricity  ratio  and  three  predetermined  values 
of,  Reynolds  number.  Xhua,  9  separate  calculations  of  performance  character¬ 
istics  must  be  made  for  each  film  (18  calculations  altogether).  The  Reynolds 
number  for  the  inner  film  is  defined  as 


2«p(M1  -  W2)>1C1 


while  the  Reynolds  number  for  the  outer  film  is  defined  as 


(15) 


Note  that  aince  the  overall  Reynolda  number  Re  ia  specified,  Rc^  and  Re^  *r* 
determined  uniquely  by  the  ring  speed  ratio  ^2^1  *  Typically,  the  three  values 
of  Re^  and  Re^  for  which  film  characteristics  sre  calculated  are  those  correa- 
ponding  to  ^/N^  -  0.25,  0.35  and  0.45.  Typical  values  of  eccentricity  ratio 
for  which  the  inner  film  characteristics  are  determined  are  ■  0.2,  0.3  and 
0.5.  The  three  values  of  selected  for  calculation  of  outer  film  data  are 
chosen  in  accordance  with  anticipated  operating  eccentricities  of  the  outer 
film. 


3.  Once  the  tables  of  inner  and  outer  film  characteristics  sre  prepared,  the 
program  next  considers  an  initial  guess  for  the  ring  speed  ratio  This 

initial  gueas  is  read  in  as  input  to  the  program.  From  this  ring  speed  ratio, 
the  program  then  calculates  a  value  for  Re^.  Corresponding  to  this  value  for 
Re^,  and  the  value  of  tj  read  into  the  program,  there  will  be  specific  values 
of  S.  and  (T„)  which  the  program  will  determine  from  the  tabular  data  for  the' 
inner  film  characteristics.  The  program  will  interpolate  within  the  tables  if 
necessary. 


4.  The  program  next  determines  &2  from  the  condition  that  the  load  capacities 
of  the  inner  and  outer  film  must  be  equal.  This  condition  la  expressed  by 


*1  + 
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Ths  program  also  detsrmlnes  Re2  corresponding  to  the  guessed  value  of  ring  speed 
ratio. 

5.  With  and  R*g  calculated,  ths  program  next  determines  the  corresponding 
values  of,  and  from  ths  tabular  data  for  ths  outer  film  performance 

characteristics.  The  program  then  checks  to  see  if 

C1 

<TJ  -  (T.)  (17) 

i  2  B  1  C2 
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1.*.  the  program  checks  to  aee  if  the  torque  oa  the  outaide  of  the  floating  ring 
matches  the  torque  on  the  inaide.  Xf  they  match,  the  aolution  it  complete.  If 
the  torque  on  the  outaide  of  the  ring  ia  lower  (higher),  then  a  slightly  higher 
(lower)  value  for  ring  a peed  la  guessed  and  the  process  is  repeated  until  a  con¬ 
vergent  solution  la  obtained. 

Vhen  a  convergent  aolution  Is  obtained  for  the  floating  ring  bearing,  the  program 
prints  out  various  performance  data  for  the  individual  films  and  for  the  overall 
floating  ring  bearing.  This  output  ia  described  fully  in  Appendix  IV.  Detailed 
instructions  for  preparing  the  Input  for  the  floating  ring  program  are  also  pro¬ 
vided  in  this  appendix  along  with  a  listing  of  the  program. 

Itiffneaa  and  Damping  Coefficients  for  floating  Ring  flmiai 

In  order  to  determine  the  overall  stiffness  end  damping  coefficients  for  the 
floating  ring  bearing,  It  la  neceaaary  to  first  determine  the  atlffneaa  and  damp¬ 
ing  coefficients  for  each  Individual  btarlng  film,  for  each  steady  state  aolution 
for  the  floating  ring  bearing,  the  steady  state  operating  conditions  for  each 
film  are  established.  Stiffness  and  damping  coefficients  for  each  film  can  there¬ 
fore  be  determined  by  the  single-film,  spiral- grooved  journal  bearing  program  ea 
described  earlier. 


Let  us  denote  the  atlffneaa  and  damping  coefficients  for  the  inner  film  by  the 
subscript  1  and  those  for  the  outer  film  by  the  subscript  2,  i.e.,  Kxxl*  *xx2’ 
ate.  The  overall  atlffneaa  and  damping  coefficients  for  the  floating  ring  bearing 
are  denoted  simply  as  K^,  K  ,  ctc*  Consider  that  the  ehaft  moves  in  sychronous 
whirl  with  e  frequency  -  2flB^  radlana/eec.  end  with  amplitude  component* 
x^a^l*.  The  overall  damping  and  stiffness  coefficients  are  defined  by  the  fol¬ 
lowing  relationships . 


r  elv 

_* _ 


0. 
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or  in  short: 


Z  2 
xx  xy 


Z  Z 
yx  yy 


where: 


C.ai.B 
-1  1  » 


(analogous  for  Z  ,  Z  ,  Z  ) 
xy  yx’  yy  ■ 


1  c. 


*1  C, 


Let  the  center  of  the  ring  have  whirl  amplitude*  x2ei“lt  and  and  let: 


2  C, 


r2  C, 


Ihe  dimeneionle.a  dynamic  coefficient*  for  the  inner  and  outer  film  are  obtained 
in  the  form: 
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Knlvln*  for  2.2.  etc.  and  us inn  Eq. 
xx  xy 

flcients  for  the  float: ig-ring  bearing. 


(21)  yields  the  overall  dynamic  coef 


In  general  the  masa  of  the  ring  la  relatively  email  such  that: 


cl«l 


<  < 


C.K 
.  2  2US 
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in  which  case  it  can  be  ignored  in  the  calculations.  This  condition  applies 
to  all  the  numerical  results  given  in  the  present  report. 


Stability  Calculation  for  Floating  Rin*  tearing 

The  overall  stability  of  the  floating  ring  journal  bearing  to  self-excited 
whirl  may  be  calculated  from  the  overall  dynamic  coefficients  described  above. 
In  this  present  study,  the  stability  calculations  were  performed  by  stem  of 
the  computer  program  daveloped  under  USAF  contract  No.  AF  33(615)-3238  and 
described  in  pert  V  of  the  flnai  documentation  issued  under  this  contract 
(Xcf.  9).  A  brief  description  of  the  analysis  upon  which  the  atsbility  cal¬ 
culations  are  baaed  is  provided  below. 

At  any  given  rotor  speed  and  with  a  known  static  load  on  the  bearing,  the 

journal  center  occupies  a  certain  unique  equilibrium  position  relative  to 

the  bearing  center.  When  the  journal  whirls  around  this  equilibrium  in  a 

small  orbit,  the  dynamic  forces  F  and  F  generated  in  the  bearing  fluid 

x  y 

film  can  be  expressed  in  linearized  form  os: 


F  - 

x 


F  - 

y 


Kx-B  -  K  y  *  B 

xx  xx  dt  xy  xy  dt 


K  x-B  -  K  y  -  B  ^ 

yx  yx  dt  yy  yy  dt 


(34) 

(35) 
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where  x  and  v  are  the  whirl  mmiiMiH*.  ae-rurcd  frc=  the  ; title  e4wX1.iu1.iu1a 
position,  t  is  tine,  and  the  four  spring  coefficients  (the  K  -  coefficients) 
and  the  four  damping  coefficients  (the  B  -  coefficients)  would  be  deterained 
for  the  floating  ring  bearing  from  the  analysis  described  above.  For  a  given 
bearing  geometry  and  known  lubricant  properties,  the  8  coefficients  are  func¬ 
tions  of  the  bearing  lead  and  the  rotor  speed  and,  if  the  lubricant  la  compress¬ 
ible  like  a  gas,  they  are  also  functions  of  the  whirl  frequency.  In  the  latter 
case,  Eqs.  (34)  and  (35)  are  only  valid  for  harmonic  motions  such  that: 

x  ■  xc  cos  (vt)  -  xB  sin  (vt)  t36) 

y  *  yc  cos  (vt)  -  yB  sin  (vt)  (37) 


where  v  is  the  angular  whirl  frequency.  These  equations  can  also  be  written: 


x 


y 


Re  (<xc  +  ix<)eivt J  - 

R*  ^(yc  +  iys)«lvtJ 


(38) 


(39) 


where  i  *  V-l  and  'fce  (  means  that  only  the  real  part  of  the  bracketed 
expression  applies.  For  convenience  the  "Re  ^  and  the  eivt  are  dropped 
whereby  Eqs.  (38)  and  (39)  are  written 


x  ■  x  +  ix 

c  s 


(40) 


y  -  yc  +  iys 


(41) 
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When  these  equations  srs  used  in  ths  analysis  their  complete  meaning  is  defined 
through  Kqs.  (38)  and  (39). 

/ 

With  this  notation  Eqa.  (34)  and  (35)  can  be  written: 


(42) 


where: 


Z 

xx 


+  1  (£)  CUB  -  K  +  iycuB 
No  xx  xx  a 


(43) 


(44) 


and  similarly  for  Z  ,  2^  and  Zyy.  Here,  cu  la  the  angular  speed  of  rotation 
and  7  gives  the  ratio  between  the  whirl  frequency  and  the  rotational  frequency. 
In  this  form,  the  aquations  are  equally  valid  for  an  inconq. resslble  and  a  com¬ 
pressible  lubricant. 


To  illustrate  the  procedure  for  calculating  the  threshold  of  instability,  assume 
for  simplicity  that  the  rotor  is  rigid  and  symmetric  such  that  the  two  bearings 
support  an  equal  mass  M  which  equals  half  the  mass  of  the  rotor.  Then  the  equa¬ 
tions  of  motion  for  a  journal  become: 


(46) 
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By  substitution  from  Eqs ■  (41)  through  (43),  these  equations  can  be  written  in 
matrix  form: 


These  two  equations  must  be  satisfied  simultaneously  at  the  threshold  of  insta¬ 
bility.  They  contain  two  unknowns,  namely  the  whirl  frequency  ratio,  y,  and 
the  angular  speed  of  rotation,  o>.  In  the  general  case,  the  8  dynamic  fluid 
film  coefficients  are  functions  of  both  y  and  o>,  making  a  closed  form  solution 
impossible,  and  the  solution  is  most  conveniently  obtained  graphically.  For 
any  fixed  value  of  to,  A  and  A  can  be  plotted  as  functions  of  y  to  find  their 

_  C  0 

aero  points.  With  y  >  0  it  is  seen  that  A  has  one  sero  point  and  A  has  up 

f  c 

to  two  sero  points  (only  true  In  this  simple  case).  The  calculation  is  repeated 
for  several  values  of  o>  and  the  results  may  be  plotted  as  shown: 
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RESULTS 


Single  Flic  Bearing 

Results  of  calculations  of  the  performance  of  single  film  spiral-grooved  Journal 
bearings  are  shown  in  Figs*  7  through  18.  In  all  cases,  the  results  shown  are 
for  bearings  having  groove  geometry  optimised  for  maximum  radial  component  of 
stiffness  at  €  >0.  These  optimum  values  of  groove  geometry  parameters  are  given 
in  Tables  1  and  2  presented  earlier  in  the  text. 

Load  capacity  of  the  single  film  spiral-grooved  journal  bearing  is  shown  in  Riga. 
7,  8,  9  and  10,  in  terms  of  dimensionless  load  V/pN^LD  (C/R)  vs.  eccentricity 
ratio  €.  Fig.  7  shows  results  for  a  bearing  of  configuration  1  i.e.  a  bearing 
having  no  flow-through  of  lubricant  (see  Fig.  2).  L/D  ratio  ia  taken  to  be  1.0. 

It  is  assumed  that  only  the  grooved  journal  ia  rotating.  In  this  figure,  it  can 
be  noted  that  the  existence  of  turbulence  does  not  significantly  effect  the  line¬ 
arity  of  the  load  vs.  eccentricity  curve,  but  only  serves  to  increase  the  load 
capacity  over  that  which  would  be  obtained  if  flow  remained  laminar. 

In  Fig.  8  are  shown  load  va.  eccentricity  curves  for  a  bearing  of  configuration 
2  with  L/D  -  1.0,  in  which  there  is  net  flow  of  lubricant  pumped  through  the 
bearing  entirely  by  action  of  the  spiral  grooving  (F^  ■  0).  The  load  capacity 
of  this  bearing  is  seen  to  be  slightly  less  than  that  for  the  configuration  1 
bearing.  ,  On  the  other  hand,  the  through-flow  of  lubricant  is  useful  in  removing 
heat  from  the  bearing. 

Fig.  9  shows  load  capacity  of  a  configuration  2  type  bearing  with  L/D  -  0.5.  As 
can  be  seen,  unit  load  capacity,  V/LD,  decreases  significantly  with  decrease  in 
L/D  ratio.  As  a  rough  guide,  it  is  found  that  in  the  range  0.5  <  L/D  <  2.0,  unit 
load  capacity  is  very  nearly  proportional  to  L/D  ratio. 

To  provide  a  greater  flow  of  lubricant  through  a  configuration  2  bearing,  one 
can  supply  the  lubricant  to  the  grooved  end  at  an  elevated  pressure  P  .  The 
effect  of  thia  on  load  capacity  is  relatively  slight  as  shown  by  the  curves  in 
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Vlg.  10.  The  result*  are  shown  for  Re  ■  5000  but  are  typical  of  those  obtained 
at  all  values  of  Reynolds  number.  The  degree  of  pressurization  considered  is 
Indicated  by  the  dimensionless  parameter  P^  G  D/L  where  is  the  dimensionless 
supply  pressure  defined  ac  7  -  Pp  /u(N.+  H  )  1  (C/R)2  and  G  la  a  turbulent 

I  V.  |  1  O  j  X 

viscosity  correction  factor  corresponding  to  the  mean  Reynolds  No.(G  is  obtslned 

mm  X 

from  Fig.  5).  When  the  supply  pressure  parameter  P  G  (D/L)  Is  maintained  at 

S  X 

0.35,  the  net  flow  of  lubricant  through  the  bearing  due  to  prejaurlzatlon  is 
approximately  equal  to  that  due  to  self-pumping  of  the  grooves,  independent  of 
Reynolds  number.  When  this  parameter  doubles,  flow  duo  to  pressurization  doubles. 

In  general  one  nay  conclude  that  for  values  of  P  G  (D/L)  less  than  0.7,  the 

8  X 

affect  of  preasurisation  on  load  capacity  may  be  neglected.  In  any  case,  it  is 
conservative  to  do  so  since  pressurise tion  tends  to  increase  load  capacity*. 

Ylgs.  11,  12,  and  13  show  curves  of  attitude  angle  i  vs.  eccentricity  ratio  for 
configuration  1  and  configuration  2  bearings  at  different  Reynolds  numbers  and 
L/D  ratios.  The  effect  of  preasurisation  on  attitude  angle  is  also  shown  by  the 
dashed  curves  in  Figs.  12  and  13. 

Referring  to  Fig.  11,  which  gives  attitude  angle  for  a  configuration. 1  bearing, 
we  see  that  for  laminar  flow  i  decreases  slightly  with  eccentricity  ratio.  This 
decrease  is  mostly  due  to  effects  of  cavitation  in  the  bearing  film.  At  higher 
values  of  Reynolds  mnber,  very  little  or  no  cavitation  occurs  In  the  bearing  film 
out  to  t  ■  0.7  and,  as  a  consequence,  attitude  angle  shows  vary  little  dependence 
on  «.  The  extent  of  cevltation  that  doee  occur  in  the  bearing  will  be  discussed 
later  In  connection  with  predicted  bearing  torque. 


*  Preasurisation  of  the  grooved  end  of  ■  spiral -grooved  journal  bearing  produces 
hydrostatic  load  capacity  whether  the  journal  is  rotating  or  not.  Such  prea¬ 
surisation  will  help  to  promote  rotation  of  the  ring  in  a  spiral-grooved  float¬ 
ing  ring  bearing.  Getting  the  ring  to  rotate  cen  be  somewhat  of  a  problem  In 
plain,  cylindrical  floating  ring  bearings.  Whether  this  is  true  for. spiral- 
grooved  floating  ring  bearings  remains  to  be  seen. 
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For  conf't  snrjiHnn  7  bearlnga.  we  find  again  that  for  laminar  flow.  6  decreases 
with  e  due  to  cavitation  whereaa  thia  effect  la  leaa  pronounced  at  higher  values 
of  Reynolds  number.  We  also  find  that  pressurization  tends  to  decrease  attitude 
angle.  In  most  instances,  this  effect  is  not  great,  although  for  laminar  flow 

and  L/D  >0.5,  a  pressurization  of  P  G  D/L  >0.7  produces  approximately  a  10 

0  s 

degree  reduction  in  attitude  angle. 

Comparing  Figs.  11  and  12,  we  see  that  a  configuration  2  bearing  of  L/D  •  1.0  has 
a  slightly  higher  attitude  angle  than  a  configuration  1  bearing  of  the  seam  L/D 
ratio.  Comparing  Figs.  12  and  13  we  see  that  for  configuration  2  spiral-grooved 
journal  bearings,  attitude  angle  decreases  quite  substantially  as  L/D  ratio  is 
decreased  from  1.0  to  0.5.  Ihia  Latter  effect  does  not  occur  in  plain  cylindrical 
bearings . 

| 

In  general,  we  can  observe  that  development  of  turbulence  in  spiral- grooved  Journal  ' 

bearings  reduces  the  attitude  angle  by  a  significant  amount  l.e.  approximately  * 

10  degrees  at  low  eccentricity  ratios  for  both  configuration  1  and  configuration 
2  geometries. 

The  bearing  through-flow,  Q^,  that  la  generated  by  the  self-pumping  of  spiral  ; 

grooving  is  shown  in  Fig.  14.  Results  are  plotted  in  terms  of  Q^/R^C(Nj-  Nq) 
vs.  e.  The  results  shown  were  obtsined  neglecting  effects  of  cavitation.  This 
was  done  because  the  way  in  which  cavitation  is  handled  in  the  present  analysis 
does  not  provide  an  accurate  calculation  of  flow  rate  when  cavitation  appears. 

t 

This  is  not  a  serious  deficiency  since  spiral-grooved  bearings  usually  operate 
with  a  full  fiuid  film. 

As  Fig.  14  indicates,  the  "self-pumping"  flow  of  an  optimized  bearing  increases 
aa  turbulence  develops  in  the  bearing  film.  This  is  probably  due  to  the  fact 
that  the  optimum  value  of  groove  depth  increases  as  turbulence  develops.  When 
turbulence  is  fully  developed  (Re  ■  5000)  both  flow  and  optimum  pocket  depth 
approach  asymptotic  values. 

It  should  be  kept  in  mind  that  the  results  shown  in  Fig.  14  pertain  to  a  bearing 
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with  Y  “  0.67  l.e.  a  bearing  with  grooving  on  67%  of  ita  length.  Flow  through 
the  bearing  can  be  increaaed  by  increaaing  Y  • lightly  with  very  little  penalty 
paid  in  tcr^s  of  loss  of  load. 

The  flow  of  lubricant,  Q  ,  through  a  spiral- grooved  bearing  that  results  from 

•  3 

pressurising  either  end  of  the  bearing  la  directly  proportional  to  F  C  G  /fi  (D/L) 

fl  2 

for  any  given  eccentricity  ratio  and  Reynolds  number.  Curves  of  dimensionless 

3 

pressure  flow,  Q  u/G  C  (L/D)  are  plotted  vs.  c  in  Fig-  15 .  Use  of  the  turbulent 

2  2 

viscosity  correction  factor  G(  in  forming  this  dimensionless  flow  takes  account 
of  the  influence  of  Reynolds  number  quite  well,  although  some  slight  dependence 
on  Reynolds  number  still  remains. 

3 

The  dimensionless  torque,  TjC/nN^R  L,  on  the  journal  of  a  spiral- grooved  bearing 
Is  plotted  In  Fig.  16  for  a  configuration  1  bearing  and  in  Fig.  17  for  a  con¬ 
figuration  2  bearing.  The  eolld  lines  show  the  torque  for  a  bearing  with  a  com¬ 
plete  fluid  film  while  the  dashed  curve  shows  the  torque  taking  account  of  cav¬ 
itation  that  would  be  expected  to  occur.  Cavitation  Is  accounted  for  by  assuming 
that  all  regions  of  subamblent  pressure  ere  cavltated  and  that  shear  stresses  In 
these  regions  are  negligible.  The  dlscrepencles  between  the  dashed  and  solid 
curves  In  Figs.  16  and  17  provide  an  Indication  of  the  extant  of  cavitation  that 
develops  as  eccentricity  increases. 

Looking  st  Fig.  16,  we  see  that  for  a  configuration  1  bearing,  cavitation  does 
not  set  In  until  ?  0.3  for  laminar  flow,  until  e  *  0.3  for  Re  >  1000,  and  does 
not  occur  at  all  below  c  ■  0.7  for  Re  ■  5000  and  9000.  For  an  unpressurlsed 
configuration  2  bearing  (Fig.  17),  cavitation  occurs  at  lower  values  of  c  than 
for  a  configuration  1  bearing  and  does  occur  at  Re  ■  3000  and  9000.  Cavitation 
can  easily  be  eliminated,  however,  by  modest  pressurising  of  the  bearing. 

The  solid  curves  shown  In  Fig.  17  also  apply  with  reasonable  accuracy  to  bearings 
with  L/D  >0.5.  The  dashed  curves  do  not,  however, because  there  Is  less  tendency 
for  a  spiral-grooved  bearing  to  cavitate  as  L/D  ratio  is  decreased.'  For  configur¬ 
ation  2  bearings  with  L/D  ■  0.5,  no  cavitation  occurs  for  Re  >  1000  and  c  <  0.7. 


An  important  quantity  to  consider  is  how  the  dimensionless  ratio  of  friction 
torque  to  load,  XjVC.  variea  with  Reynolds  number.  This  <9  shown  in  Fig,  18 
for  e  ■  0.2,  0.5  and  e  a  0.7.  As  can  be  seen,  torque  increases  more  rapidly 
than  does  load  when  turbulence  develops  in  the  bearing  film.  One  can  note  at 
this  point  that  one  of  the  advantages  of  the  floating  ring  bearing  is  that  the 
Reynolds  number  in  each  separate  bearing  film  is  less  than  the  Reynolds  number 
that  would  be  obtained  if  the  bearing  had  only  a  single  film.  Hence,  due  to 
this  effect  alone,  the  floating  ring  bearing  can  operate  with  a  more  favorable 
torque  to  load  ratio  than  an  equivalent  bearing  with  only  a  aingle  film. 

Floating  Ring  Bearing 

The  floating  ring  configuration  choaen  for  analysis  was  one  having  an  overall 
length  to  diameter  ratio  of  L/D  -  1.0.  (See  Figure  4).  This  means  that  each 
half  of  the  floating  ring  bearing  had  an  effective  length  to  diameter  ratio 
of  0.5.  Grooving  parameters  for  each  half  of  the  floating  ring  bearing  were 
those  presented  in  Table  2.  Grooving  on  the  outside  surface  of  the  ring  was 
the  same  as  that  on  the  shaft. 

Two  values  of  the  ratio  of  inner  clearance  to  outer  clearance  were  considered 
i.e.  Cj/Cj^  -1.2  and  0.8.  Results  were  obtained  for  no  pressurization  of  the 

bearing  (P  -  0)  and  for  a  degree  of  pressurization  corresponding  to  F  G  D/L 

8  8  8 

-  0.35. 


The  static  performance  data  for  the  floating  ring  journal  bearing  are  given 
in  Table  3.  Dynamic  performance  data  are  given  in  Tables  4  and  5.  Much  of 
this  data  la  presented  in  graphical  form  in  Figs.  19  through  36.  These  fig¬ 
ures  are  discussed  below. 

Curves  of  dimensionless  load,  W,  vs  the  eccentricity  ratio  of  the  inner  film 
€ ,  are  presented  in  Fig.  19.  The  curves  are  for  Cg/Cj.  *  1*2  although  they 
apply  within  a  few  percent  accuracy  to  the  case  of  C^/C^  -  0.8.  As  can  be 
seen  the  degree  of  pressurization  considered  for  the  floating  ring  bearing 
(dashed  curves)  results  in  a  substantial  increase  in  load  capacity  over  the 
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unpreasurlsed  cate  (solid  curves).  For  the  single  film  bearing,  the  same  appar¬ 
ent  degree  of  pressurisa tion  resulted  in  only  a  slight  increase  in  load  capacity 
(See  Figure  10) .  The  reasons  for  the  greater  apparent  effect  of  pressurization 
on  the  floating  ring  bearing  are  as  follows. 

First,  the  parameter  P  G  D./L  for  the  floating  ring  bearing  is  based  on  the  over- 
all  L/D^.  The  effective  value  of  this  parameter  for  each  half  of  the  floating 
ring  bearing  is  actually  twice  the  overall  value.  Second,  the  value  for  Gf  used 
in  establishing  the  parameters  P  G^D^/L  is  based  on  the  overall  Reynolds  number 
Re  "  2*N1RjC1/v •  Because  of  rotation  of  the  ring,  the  individual  Reynolds  numbers 
for  the  inner  and  outer  films  are  each  less  than  the  overall.  Hence,  the  effect 
of  pressurization  of  each  film  is  greater  than  is  indicated  by  the  parameter 
G^  D^/L  because,  in  turbulent  flow,  the  lower  the  Reynolds  number  the  greater 
the  amount  of  flow  for  a  given  supply  pressure.  Third,  the  dimensionless  supply 
pressure  F#  is  based  on  the  shaft  speed  H^.  For  spiral* grooved  bearings,  it  turns 
out  that  rotation  of  the  ring  decreases  rather  than  increases  the  load  capacity 
of  the  inner  film.  Hence  pressurization  of  the  inner  film  becomes  relatively  more 
significant  with  respect  to  load  capacity  of  the  film. 

Roughly  speaking,  the  increase  in  load  capacity  resulting  from  pressurization  of 
a  floating  ring  bearing  is  linearly  proportional  to  the  guage  supply  pressure. 

One  can  therefore  linearly  Interpolate  between  the  curves  shown  in  Fig.  19  to 
deterorine  load  capacity  at  supply  pressures  different  from  thst  considered. 

Curves  of  overall  attitude  angle  4  of  the  floating  ring  bearing  are  shown  in  Figs- 
20  and  21.  Overall  attitude  angle  is  defined  in  Fig.  4.  Values  for  the  attitude 
angles  of  the  inner  and  outer  films  are  given  in  Table  3.  Due  to  rotation  of  the 
ring,  which  decreases  the  spiral-grooved  pumping  effect  in  the  inner  film,  atti¬ 
tude  angles  for  the  inner  film  are  considerably  greater  than  for  the  outer  film. 

Valuea  of  dimensionless  journal  torque,  T^  •  T^/VC,  are  given  in  Fig.  22  for  a 
floating  ring  bearing  with  C^/C^  ■  1.2.  Similar  curves  for  a  bearing  with 
Cg/C^  ■  0.8  would  run  about  71  higher. 
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One  of  the  primary  advantages  of  the  floating  ring  bearing  is  that,  for  a  given 
eccentricity  of  the  inner  film,  the  ratio  of  torque  to  load  is  much  lower  than 
for  a  comparable  aingle  film  bearing.  This  is  particularly  true  for  plain  bear¬ 
ings  for  which  load  capacity  of  the  inner  film  is  proportional  to  the  sum  of  the 
speeds  of  the  shaft  and  ring.  In  spiral  grooved  bearings,  however,  the  pumping 
effect  of  the  grooves  is  proportional  to  the  difference  in  speed  between  the 
shaft  and  ring.  Load  capacity  of  these  bearings  is  due  partly  to  this  pumping 
effect  and  partly  to  the  usual  hydrodynamic  effect  which  is  proportional  to  the 
sum  of  the  shaft  and  ring  speed.  The  net  effect  is  that  as  ring  speed  increasea, 
load  capacity  of  the  inner  film  for  a  spiral-grooved  bearing  decreases  slightly. 
Consequently,  spiral -grooved  floating  ring  bearings  do  not  enjoy  the  same  torque 
to  load  advantage  possessed  by  plain  floating  ring  bearings.  Nonetheless,  the 
spiral-grooved  floating  ring  bearing  does  have  a  torque  to  load  ratio  better 
than  that  of  a  single  film  bearing  operating  at  the  same  eccentricity  ratio.  This 
is  evidenced  by  the  curves  shown  in  Fig.  23.  The  single  film  bearing  used  for 
comparison  in  this  figure  is  a  configuration  2  bearing  with  L/D  *  0.5-  This  pro¬ 
vides  a  fair  comparison  because  each  side  of  the  centrally  fed  floating  ring 
configuration  we  are  considering  consists,  essentially,  of  an  isolated  bearing 
with  L/D  *0.5. 

Total  flow  pumped  through  the  floating  ring  bearing  by  the  self-pumping  effect 
of  the  spiral  grooves  is  plotted  in  Fig.  24.  The  increase  in  lubricant  flow 
that  would  result  from  pressurization  of  the  bearing  can  be  calculated  from 
the  aingle  film  curvea  plotted  in  Fig.  IS.  These  single  film  curves  can  be 
applied  directly  to  each  Individual  film  of  the  floating  ring  bearing  on  either 
side  of  the  central  feeding  groove.  In  applying  these  curvea  to  calculate  pres¬ 
sure  flow,  one  must  be  careful  to  use  the  appropriate  Reynolds  number  and  L/D 
ratio  corresponding  to  the  individual  film  being  considered.  The  pressure  flow 
that  is  calculated  can  be  added  directly  to  the  aelf-pumping  flow  calculated  from 
Fig.  24. 

Ring  speed  ratio,  N^/N^,  is  plotted  aa  a  function  of  inner  film  eccentricity 
in  Figs.  25a  and  25b.  Results  for  laminar  flow  and  Re  *  9000  are  shown.  Results 
for  Re  *  5000  are  nearly  the  same  as  for  Re  ■  9000  while  results  for  Re  *  1000 
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11^  between  the  Re  “  9000  and  laminar  curves.  Ring  speed  rstio  for  lamlnsr  flow 
rend*  to  decrease  yirh  while,  £ nr  turKijlent  flow,  Ting  speed  ratio  remains 

relatively  constant  with  eccentricity.  Note  thst  ring  speed  rstio  is  not  strongly 
affected  by  the  clearance  ratio  Cj/C^. 


fc2 ,  the  eccentricity  ratio  fer  the  outer  film,  is  plotted  as  a  function  of  in 
Fig.  25.  In  general,  e2  increases  approximately  linearly  with  e^.  As  would  be 
expected,  considerably  greater  than  for  02/0^  *  1.2  and  more  nearly 

equal  to  for  C^/C^  «  0  8. 


Dynamic  data  for  the  floating  ring  bearing  are  given  in  tables  4  and  5.  These 
data  consist  mostly  of  stiffness  and  damping  coefficients  for  the  individual 
bearing  films  and  for  the  overall  bearing.  The  data  also  include  values  of 
the  dimensionless  critical  mass  for  the  threshold  of  whirl  instability. 


Data  for  the  individual  bearing  films  are  given  in  Table  4.  Stiffness  and  damp¬ 
ing  for  the  inner  film  pertain  to  the  fluid  film  forces  that  develop  in  the  inner 
film  due  to  relative  motion  between  the  shaft  and  the  ring.  Stiffness  and  damping 

v 

for  the  outer  film  pertain  to  fluid  film  forces  developed  in  that  film  due  to 
relative  motion  between  the  ring  and  the  outer  bearing. 

Stiffness  and  damping  coefficients  for  the  overall  floating  ring  bearing  are 
given  in  Table  5.  These  pertain  to  the  forces  developed  on  the  shaft  due  to  a 
relative  motion  between  the  shaft  and  the  outer  bearing  including  the  effect  of 
motion  of  the  ring. 

Overall  stiffness  and  damping  coefficients  for  the  floating  ring  bearing  are 
plotted  vs.  for  a  r umber  of  representative  situations  in  Figs.  27  through 
32.  Qualitatively,  the  behavior  of  the  overall  stiffness  and  damping  coeffici¬ 
ents  as  eccentricity  ratio  increases  is  similar  to  that  for  single  film  bearings. 

The  stiffness  and  damping  coefficients  presented  in  Tables  4  and  5  may  be  used 
to  calculate  a  critical  mass  at  the  threshold  of  whirl  instability.  Tne  critical 
masses  determined  from  the  inner  or  outer  film  damping  and  stiffness  coefficients) 
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while  not  of  great  physical  significance,  are  still  of  interest  to  calculate. 
The  critical  mass  Mcl,  determined  from  the  inner  film  coefficients,  la  the 
critical  mass  of  the  shaft  within  the  rotating  ring  assuming  that  the  ring 
were  restricted  from  any  translational  motion  i.e.  assuming  that  the  outer 
film  were  infinitely  stiff.  M^t  the  critical  mass  calculated  from  the  outer 
film  coefficients,  is  the  critical  mass  for  the  ring  rotating  within  the  bear¬ 
ing  neglecting  any  effect  of  the  inner  bearing  film.  Essentially,  repre¬ 
sents  the  critical  mass  for  a  single  film  bearing  operating  at  Reynolds  number 
Re^  and  eccentricity  ratio  e^. 

M  ,  the  critical  mass  determined  from  the  stiffness  and  damping  coefficients 
for  the  overall  floating  ring  bearing,  represents  the  critical  mass  for  the 
shaft  rotating  within  the  composite  floating  ring  structure.  It  is  the  ap¬ 
propriate  value  of  critical  mass  for  the  shaft  of  the  floating  ring  bearing. 

The  critical  masses  M  , ,  M  .  and  M  are  presented  in  Tables  4  and  5  in  the 

Cl  Ci  c 

dimensionless  form 


Mcl  C1  <W1  +  V 

^<»1/C1)2  LDl 


Mc2  C2  N2 

u(r2/c2)2  ld2 


M  C.  N. 

M  -  — c.  -LJ — 

c  mCRj/CjT  ID1 

In  attempting  to  compare  the  critical  massea  in  these  tables  one  should  keep 
in  mind  that  each  is  made  non-dimensionlecs  on  a  slightly  different  basis. 


Values  of  overall  critical  mass  vs.  are  plotted  in  Figs.  33  and  34  for 


Cj/C^  -  1.2  and  0.8  respectively .  In  general,  critical  mass  for  hydrodynamic 
bearings  increaaea  *’ith  eccentricity  ratio.  However,  when  operating  in  the 
turbulent  regime,  *«l-grooved  bearings  exhibit  the  anomaloua  characteristic 
that  critical  mass  increases  with  eccentricity  in  certain  ranges.  This  char¬ 
acteristic  may  be  related  to  the  fact  that  the  attitude  angle  of  spiral- grooved 
bearings  sometimes  increaaea  with  eccentricity  ratio  in  the  turbulent  regime. 
This  increase  in  attitude  angle  can  be  qualitatively  explained  on  the  basis 
that  as  eccentricity  increases,  the  usual  plain-bearing- type  of  hydrodynamic 
actioa  becomes  relatively  more  significant  compared  with  the  self-pressurising 
hydrodynamic  action  of  the  spiral  grooves,  particularly  in  turbulent  flow. 

Hence,  there  ia  a  tendency  for  attitude  angle  to  increase  with  eccentricity 
in  spiral" grooved  hearings. 

It  la  interesting  to  note  in  Figs.  33  and  34  that  preaeurlsatlon  of  the  floating 
ring  bearing  has  a  more  algnlf leant  effect  on  bearing  stability  than  it  does  on 
bearing  load  capacity.  It  la  obvious  that  pressure  feeding  of  the  bearings  ia 
an  effective  means  of  stabilisation. 

An  important  point  to  investigate  la  whether  the  spiral-grooved,  floating  ring 

bearing  configuration  ia  more  or  leas  stable  than  an  equivalent  spiral-grooved, 

single  film  bearing.  Comparisons  of  the  critical  masses  of  the  floating  ring 

and  single  film  configurations  are  shown  in  Figs.  35  and  36  for  laminar  and 

turbulent  flow  respectively.  The  single  film  vslues  of  critical  mass  were  taken 

from  the  values  for  M  _  calculated  for  the  outer  film.  The  L/D  ratio  of  the 

c2 

outer  film  ia  about  20%  less  than  that  for  the  inner  film  but  this  fact  should 
not  greatly  effect  the  comparison. 

For  laminar  flow,  the  single  film  bearing  is  more  stable  than  the  floating  ring 
bearing  with  Cj/C^  ■  0.8  and,  at  low  eccentricities,  is  also  more  stable  flhan 
the  floating  ring  bearing  with  C^/C^  >1.2.  However,  at  high  eccentricities, 
the  critical  mass  for  the  floating  ring  bearing  with  Cg/C^  -  1.2  exceeds  that 
of  the  single  film  bearing. 

For  turbulent  flow,  the  stability  of  the  single  film  configuration  is  clearly 
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“•jpericr  tc  thst  cZ  the  floating  ring  oearing. 

The  floating  ring  configuration  of  spiral-grooved  bearing*  *uffera  in  regard  to 
stability  for  the  following  reaaon.  The  good  atability  of  apiral-grooved  bearing* 
is  due  principally  to  the  self-pressurizing  effect  of  the  apiral  groove*.  This 
self-pressurizing  effect  is  proportional  to  the  difference  in  speeds,  (N^  -  N^) . 
Consequently,  in  a  floating  ring  configuration,  rotation  of  the  ring  reduces 
this  self-pressurizing  effect  in  the  inner  film  and,  hence,  reduces  the  stability 
of  the  inner  film.  The  stability  of  the  overall  bearing  suffers  a a  a  consequence. 

Stability  of  the  apiral-grooved  configuration  is  compared  with  the  stability  of 
a  plain  floating  ring  bearing  in  Fig.  37.  At  low  eccentricities,  stability  of 
the  spiral-grooved  configuration  is  better  but,  at  high  eccentricity,  stability 
of  the  plain  bearing  increases  rapidly  due  to  the  effect  of  cavitation  and  begins 
to  surpass  that  of  the  spiral-grooved  bearing. 

An  important  point  to  keep  in  mind  when  comparing  the  stability  of  plain  and  . 
spiral-grooved  bearings  is  that  the  plain  bearings  achieves  atability  only  as 
a  result  of  cavitation  in  the  bearing  film.  A  plain  bearing  operating  with  a 
full  fluid  film  is  Inherently  unstable  for  any  speed  or  mass.  Often,  particul¬ 
arly  with  liquid  metals  aa  lubricants,  it  la  undesirable  to  operate  with  cavit- 
ated  bearing  films  because  of  the  problem  of  cavitation  damage.  If  a  pressur¬ 
ized  supply  is  used  with  plain  bearings  to  suppress  cavitation,  the  stability 
of  the  bearings  suffers  drastically. 

/ 

On  the  other  hand,  spiral-grooved  bearings,  even  when  unpressurized,  operate 
without  cavitation  in  the  bearing  film  out  to  quite  large  eccentricity  ratios 
due  to  the  self-pressurising  action  of  the  spiral  grooves.  Stability  of  these 
bearings  is  achieved  through  this  self-pressurizing  action.  Moreover,  use  of 
a  pressurized  lubricant  supply  further  enhances  the  stability  of  these  bearings 
as  is  evidenced  by  the  performance  charts  presented  in  this  report. 
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SUMMARY  AND  CONCLUSIONS 


An  analysis  has  bean  performed  of  the  turbulent,  spiral-grooved.  Journal  bearing 
with  incompressible  lubricant.  Optimum  values  of  groove  parameters  have  been 
deteradLned  to  provide  maximum  radial  stiffness.  Performance  charts  are  presented 
for  load  capacity,  attitude  angle,  lubricant  flow  rate  and  bearing  frictional 
torque  for  both  single  film  and  floating  ring  configurations  of  the  spiral -grooved 
journal  bearing.  In  addition,  stiffness  and  damping  coefficients  are  presented 
for  the  floating  ring  configuration  and  values  of  the  critical  maaa  for  threshold 
of  whirl  instability  are  determined.  Performance  data  are  presented  for  bearings 
operating  with  and  without  a  pressurised  supply  of  lubricant. 

The  following  general  conclusions  can  be  drawn  concerning  the  performance  of 
spiral-grooved  single  film  and  floating  ring  Journal  bearings: 

1)  Compared  with  plain,  cylindrical  Journal  bearings,  spiral-grooved  bear¬ 
ings  offer  the  following  advantages:  (a)  They  possess  greater  stability 
under  lightly  loaded  conditions,  (b)  They  tend  to  operate  without 
cavitation  due  to  the  self-pressurising  effect  of  the  spiral-grooving. 

(This  is  an  important  consideration  when  operating  the  bearing  with 
liquid  metal  lubricants  where  the  problem  of  cavitation  damage  is  sig¬ 
nificant.)  (c)  The  spiral  grooving  provides  self-pumping  of  lubricant 
through  the  bearing  eliminating  the  need  for  a  pressurized  supply,  (d) 

All  performance  characteristics  of  spiral-grooved  bearings  can  be  easily 
enhanced  by  use  of  a  pressurized  supply  of  lubricant. 

2)  The  floating  ring  configuration  of  spiral-grooved  bearing  operates  with 
lower  torque  to  load  ratio  than  is  achieved  with  a  similarly  loaded 
tingle  film  bearing.  On  the  other  hand,  stability  of  the  single  film 
spiral -grooved  bearing  la  generally  better  than  that  for  the  floating 
ring  bearing. 

3)  In  spiral-grooved  bearings,  development  of  turbulence  results  in  an 
Increase  in  frictional  torque,  an  increase  in  load  capacity,  an  increase 
in  stability  and  a  decrease  in  attitude  angle.  In  general,  the  ratio  of 
frictional  torque  to  load  capacity  increases  with  development  of  turbulence. 
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Configuration  2  (Vented  Bearin.) 
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Configuration  3  (Vented  Bearing) 


Figure  2 .  Various  Configurations  of  Spiral-Grooved  Journal  Bearings 


Coefficient  of  Friction  vs.  Local  Reynolds  Number 
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Figure  7.  Load  V6 .  Eccentricity,  Single  Film  Bearing 
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Figure  8.  Load  vs.  Eccentricity,  Single  Film  Bearing 
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Figure  9.  Load  vs.  Eccentricity,  Single  Film  Bearing 
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Figure  11.  Attitude  Angle  vs-  Eccentricity,  Single  Film  Bearing 
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Figure  22*  Torque  to  Load  Ratio  vs .  Inner  Film  Eccentricity,  Floating  Ring  Bearing 
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Figure  23*  Torque  to  Load  Ratio,  Comparison  between  Single  Film  and  Floating  Ring  Bearing 
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Figure  27.  Stiffness  Coefficients,  Floating  Ring  Bearing 
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Figure  30.  Damping  Coefficients, 
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Figure  31.  Stiffness  Coefficients,  Floating  Ring  Bearing 


Figure  32 •  Damping  Coefficients ,  Floating  Ring  Bearing 
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Figure  33.  Critical  Journal  Mass,  Floating  Ring  Bearing 
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Figure  36.  Critical  Journal  Mass,  Comparison  between  Plain  and  Spiral- Grooved 
Floating  rang  Beatings 
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Figure  38.  T|  end  £  Coordinates  for  Grooved  Surface 
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APPENDIX  I 

DERIVATION  AND  SOLUTION  OF  EQUATION  VCR 


i 


TURBULENT  SPIRAL-GROOVED  JOURNAL  BEARING 


i 

1 


The  differential  equation  for  the  pressure  distribution  around  a  spiral* grooved 
Journal  bearing  has  been  derived  previously  by  the  present  authors  for  the  case 
of  laminar  flow  (Ref.  1).  In  this  appendix  will  be  presented  a  re-derivation 
of  this  equation  in  which  the  effects  of  turbulence  in  the  bearing  film  will  be 
accounted  for  by  means  of  the  linearized  turbulent  lubrication  theory  developed 
by  Ng  sad  Pen  (Ref.  4).  The  derivation  will  Involve  four  major  steps: 

1.  Express  the  local  mass  flows  in  the  ridge  and  groove  regions  in  terms 
of  the  local  pressure  gradients  in  those  two  regions. 

2.  Define  an  ‘"overall"  pressure  profile  around  the  bearing  neglecting  the 
zig-zag  ripples  in  the  profile  which  arise  due  to  the  dlscountinuous 
groove-ridge  geometry. 

3.  Use  the  requirement  thet  the  flow  normal  to  a  groove-ridge  interface  be 
continuous  across  the  interface  in  order  to  solve  for  local  groove-ridge 
pressure  gradients  in  terms  of  the  "overall"  pressure  gradient. 


4.  Apply  the  principle  of  conservation  of  mass  to  obtain  a  differential 
equation  for  the  overall  pressure  profile. 

According  to  the  linearized  turbulent  lubrication  theory  of  Ng  and  Pan,  the  local 
turbulent  mass  flows  in  a  bearing  film  can  be  expressed  as 


puh 


ph3G 

_ _ X 


(WO 

2 


pwh 


ph3G 


& 

dz 


(51) 


(52) 
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where  x  end  z  ere  the  coordinates  in  the  direction  of  rotation  and  the  axial 
direction  respectively  and  where  U  and  V  are  the  surface  velocities  of  the 
bearing  and  journal  respectively.  The  factors  G  and  G  can  be  considered, 

X  X 

essentially,  to  be  turbulent  viscosity  correction  factors.  In  the  linearised 

turbulent  lubrication  theory  G  and  G  are  considered  to  be  functions  of  the 

local  Reynolds  number,  in  the  bearing  film  where  *  p(V-U)h/u.  A  plot 

of  G  and  G  vi.  L  ii  shown  in  Fig.  5. 
x  z  n 

Next,  we  introduce  the  "skewed*1  {,  T)  coordinate  system  shown  in  Fig.  38  in 
which  lines  parallel  to  groove«ridge  interfaces  are  lines  of  constant  £  while 
planes  perpendicular  to  the  axis  of  the  bearing  are  planes  of  constant  t). 

The  relationships  between  £  and  T|  and  the  cylindrical  coordinates  9  and  z 
are  given  below 


s  cot 


z/sin3 


,n  «UL  6 


s  ■  tj  sin  p 


note  also  that 


hi  _L_ 

(Jt)  sin  P 


St  i 


(60) 


In  the  £,  rj  coordinate  system,  the  pressure  gradient  dP/dn  Is  continuous 
everywhere  in  the  bearing  film  because  the  geometry  has  no  discontinuities 
in  the  T]  direction.  On  the  other  hand,  the  gradient  c)P/d£  is  discontinuous 
at  groove-ridge  interferes  du#  to  the  discontinuity  in  film  height.  Con¬ 
sequently  the  pressure  profile  in  the  £  (clrcumferentlal)direction  will  hsve 
tne  "sig-sag"  appearance  as  shown  schematically  by  the  solid  line  in  Fig.  39  , 
In  this  figure  the  symbols  £  ,  £. . ,  etc.,  denote  the  interfaces  at  the  begin- 
ning  of  ridge  region  while  ^n+3/2’  etc'  denote  t*ie  lnter^*ces  at  the 

beginning  of  groove  regions.  Now, by  neglecting  the  saw-toothed  ripples  in  the 
actual  pressure  dictribution,  one  can  conceive  of  an  approximate,  smoothed 
"overall"  pressure  distribution  around  the  Journal  shown  in  Fig.  39  as  the 
dashed  line  through  the  pressures  pn.3/2»  Pn-l/2*  etc*»  at  ^n-3/2’  ^n-l/2* 
etc.  Siuce  the  local  pressure  gradients  within  each  groove  and  ridge  region 
are  bounded  in  magnitude £ha  saw-toothed  fluctuations  in  pressure  due  to  al¬ 
ternating  groove  and  ridge  regions  will  reduce  to  a  negligible  magnitude  as 
the  width  of  the  groove-ridge  pair  becomes  very  amall,  i.e.,  aa  the  number  of 
grooves  becomes  very  large.  In  the  limit,  as  the  width  of  each  groove-ridge 
pair  becomes  very  small,  the  smooth  "overall"  pressure  distribution  through 

the  discrete  points  P  P  ,  etc.,  will  approach  r  continuous  distri- 

_  n-J/2  n-i/i 

bution  P(£,t)>  which  should  provide  a  very  good  approximation  to  the  actual 
"saw-toothed"  pressure  distribution  around  the  journal.  Formally,  one  can 
define  the  slope  of  P(£,q)  in  the  £  direction  at  the  point  £r,  as 


where  a£  ■  £n+i/2  '  ^n-1/2  f£iC*  vhere  d*g/d£  ,nd  refer  to  the  local, 

pressure  gradients  within  the  groove  and  ridge  region  respectively. 

Having  defined  the  overall  gradient  d?/d£i  we  next  consider  writing  an  ex¬ 
pression  for  the  mass  flux  normal  to  a  moving  groove-ridge  interface 
(see  Fig.  39).  In  terma  of  the  mass  fluxes  puh  and  pwh,  is  expressed  as 


|^(u-V) 


■in  {)  -  w  cos  £  h 


Substituting  for  u  snd  v  in  Eq.  (62)  by  Msns  of  Eqi.  (51)  «nd  (52)  and  noting 
that  dP/dx  ■  dP/KdS  we  obtain 


groove  region 


■-4^ 


h„  c  a?  „v 

-8 - 2t&  —a  .  lU^O.  h  aln  ft 


h  G _ dP_ 


u  da 


coa  0 


^Idg.  region-  •  p  ([“V^  I*  “  ^  hr],tn  * 


h  3  G 


aT co#  e  > 


Rote  that  the  factors  G^  and  G^  are  aubacripted  r  and  g  because  they  have  dif¬ 
ferent  values  in  the  ridge  and  groove  regions. 


Rext  we  nake  the  following  useful  definitions 


G.  ■  12  (G  +  G  cot' 

i  X  C 


2a.  L2 


G2  * 


12Gg  coa  P 
R  a In2 0 


G3  ‘ 


2 

sinp 


If,  In  Eqs.  (63)  and  (64),  we  transform  the  derivatives  of  P  with  respect  to 
9  and  z  into  derivatives  with  respect  to  £  sn<i  «  (»<>»  Ena.  (57)  and  (58)) and 
collect  terms,  we  obtain 

it  h«3  [°18  If  +  *28  I?] 

(68) 


ridge  region 

+ 

Note  that  the  derivative  dP/dn  doea  not  have  to  be  identified  by  a  subscript 
g  or  r  since  it  la  continuous  everywhere  within  the  bearing  film. 

Next  we  note  that  by  continuity  of  maas  flow 


-  It  sin  0 

groove  region 


■ 


+  ph 


(U-V)  1 
g  2R 


«  ||^ 

groove  region  ridge  region 


(70) 


Eqs.  (61)  and  (70)  constitute  two  linear  equations  in  the  "unknowns" 
and  3P r/3{.  Eqa.  (61)  and  (70)  may  therefore  be  solved  to  yield 


(l-o)  Bx 


(l-ad  b2 


(71) 
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(72) 


N1» 


Eqa.  (71)  and  (72)  may  be  aubatltuted  in  either  Eq.  (68)  or  (67)  to. yield  an 
expression  for  the  mass  flux  in  terms  of  the  overall  pressure  gradients 
and  dP/dTi  (note  that  oP/d>1  may  be  written  as  dP/dl  because  of  the  con¬ 
tinuous  nature  of  this  derivative).  The  expression  obtained  is 


94- 


(74) 


“»!  $ 


a  b2) 


Next,  consider  the  control  volume  AV,  for  one  ridge-groove  pair  shown  in  Fig.  41. 
We  take  the  control  volume  to  move  with  the  grooved  surface  at  velocity  V  and  ve 
consider  that,  at  the  instant  shown  in  Fig.  41,  the  various  ridge-groove  inter¬ 
faces  are  located  at  £n»  5^/2*  ^n+l  etc'»  ••  depicted •  The  mass  flows  entering 
and  leaving  the  moving  control  volume  through  the  surfaces  we  denote  as 

Hrid  e  ^  whereas  the  mass  flows  entering  and  leaving  the  control  volume  through 

®  n  n 


the  surfaces  AS  1  we  denote  as  W  (£ 
the  mass,  AM  contained  in  AV  la  given  by 


h+1 


«n>‘ 


The  total  time  rate  of  ch  nge  of 


a_ 

at 


(AM)  +  V  .  V  (<2M) 


(75) 


By  the  conservation  of  masa  we  obtain  that 


w; 


ridge 
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•n+l 


ridge 


An  + 


w 


-  w 


n+An 


^n+l  ^n 


+  It  +  V  .  7  (AM)  -  0 


(76) 


Dividing  Eq.  (76)  through  by  A$An  where  A^  -  $n+1  -  and  taking  the  limit  as 
A£  and  An  go  to  zero  we  have 


+ 


(AM) 
St  a4  An 


V 


3— toB)  ■  - 
7  a€an 


(77) 
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where  aw*/#  i»  defined  to  be 


A 


lim 

dr o 


A  l 


»tH»l 


aT 


(78) 


The  subscript  ridge  haa  been  dropped  from  W ^  became  of  the  equality  condition 
expressed  by  Eq.  (70).  The  expression  for  la  given  by  equation  (74). 

The  expreaalon  fox  the  mass  flux  W1*  la 


(l-ot)h  dP  ah  /  dP, 

- - £-  r-1  c  +  vJ 

M  da  sr  u 


a*  G*g 


(79) 


If  we  tranaform  d^/da  and  dP^/ds  Into  f  -  tj  coordinates  by  means  of  Eq.  (58)  we 
will  obtain  an  expreaalon  for  V1  in  terms  of  dt/di)>  d?r/d6  and  d*g/d|.  He  then 
can  substitute  for  d*r/d{  and  £P  /d$  by  means  of  Eqa.  (71)  and  (72)  to  obtain, 
finally 


t  ain  d 


12m 


{“  ^  [°3.  !5  '  °2g  +  ^  \  +  ll-“>  V] 


+  a-o)ht3  [o3r  fj,  -  0Jt  a 


2  4'oii  “-“V 


(80) 


where  A^,  and  ware  defined  previously  by  Eqa.  (73) 
The  quantity  AM/A6&)  in  Iq.  (77)  la 


■pi  sin  p 


£(i-a)hr  +  afrg  J 


(81) 
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while 


5  • 9  -  *  If  «» 

Substitution  of  Eqs.  (74),  (80),  (81)  end  (82)  into  Sq.  (77)  yields  a  second- 
order  differential  equation  in  P({,T))  the  overall  pressure  distribution  around 
the  spiral-grooved  journal*  To  solva  this  equation,  we  must  first  transform 
it  back  into  the  orthogonal  <3^s  coordinate  system.  This  is  done  by  means  of 
Eqs.  (59)  and  (60).  The  result  is 

4*  SP  i  £ 

+  (tr  +  R  |«)  (Ro  ,in*  [«V  *°  <83) 

Expressions  for  and  in  the  8-z  coordinate  system  may  be  obtained  by 
substituting  Eqs.  (59)  and  (60)  into  Eqs.  (74)  and  (80). 


The  Numerical  Solution 

Eq.  (83)  is  solved  for  any  combination  of  groove  and  seal  arrangement  with  or 
without  feed  or  vent  in  the  middle  of  the  Journal  (Fig.  2).  Thus,  it  can  pro¬ 
vide  the  solution  for  a  seal-groove,  groove-seal,  fully  grooved  or  herringbone 
journal  containing  up  to  three  distinct  sections.  The  basic  technique  is  first 
to  divide  the  bearing  into  a  numerical  grid  with  dimensions  mxn.  Then  write 
the  differential  equation  into  a  finite  difference  form  of  three  columns.  Thus, 
each  point  in  concern  ia  related  to  the  five  neighboring  points.  At  the  bound¬ 
aries,  the  pressures  are  given.  Therefore  mxn  equations  are  established  for 
mxn  unknowns.  By  means  of  the  columnwise  matrix  inversion  solution  routine 
developed  by  Castelli,  Pirvics  and  Shapiro,  the  pressure  field  is  obtained. 
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This  technique  is  fully  discussed  in  Esf  sr*rsces  6  and  7 . 

Once  the  pressure  field  is  obtained,  loads,  moments,  torques,  etc.  ere  obtained 
by  numerically  evaluating  the  following  expressions. 

Radial  component  of  force  (cosine  component) 


r 


r 


(P  -  P^)  R  cos  6  d0  ds 


(84) 


Tangential  component  of  force  (sine  component) 


(P  -  Fa)  R  sin  0  d0  ds 


(85) 


Radial  component  of  moment  (cosine  component) 


M_ 


(P  -  Pa)  Re  cos  0  d0  ds 


(86) 
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Tangential  component  of  moment  (sine  component) 
2* 


0 

Attitude  angles 


(P  -  P  )  Ra  sin  0  d0  da 
a 


where  C£  la  the  Couette  friction  factor  which  la  plotted  In  Pig.  6  against 
&P/d0  la  the  local  pressure  gradient.  By  Iqs>  (71),  (72),  (59)  and 
(60)  the  second  term  of  Equation  (89),  the  Poiaeullle  torque,  can  be  written 
in  terms  of  the  overall  pressure  gradients. 


(l-o)B, 
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(l-o)h 


*a 


+  cpA_a 
^  & 


|§>  Sin  9 


Journal  torque 


T 


J 


Ip  +  cCU  sin  i 


(91) 


Where  W 


load  - 


Flow 


0 


where  U1*  ie  defined  by  Eq.  (80). 
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DERIVATION  OF  RELATIONSHIPS  FOR 
CALCULATING  STIFFNESS  AND  DAMPING  COEFFICIENTS 

Consider  Che  reference  axes  shown  in  Fig.  3.  The  relationship  between  the  forces 

F  and  F  and  the  forces  F  and  L  can  be  written  in  matrix  form  as 
x  y  r  t 


For  an  infinitesimally  small  motion  around  the  steady  state  position  the  dynamic 
forces  become 


The  infinitesimal  dynamic  motion  of  the  journal  center  is  described  by  the 
coordinates  (x,y) : 

x  ■  d(e  cos  4)  y  ■  d(e  sin  i) 

or 


‘de 


cos  4  sin  4* 


ed4  J  vrsln  i  cos  4 


x 


(95) 
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The  velocities  transform  similarly: 


(96) 


If,  in  the  differential  equation  for  pressure  P  obtained  in  Appendix  I,  we  intro 
duce  the  following  diaensionlesa  variables 
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we  obtain  the  result  that  the  dimensionless  pressure  P  in  a  spiral- grooved  bear 

* 

lug  of  fixed  geometry  is  a  function  only  of  ths  dimensionless  variables  «,  c/m 
and  d/m  i.e. 


P 


P 


(98) 
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The  resulting  fluid  film  forces  in  radlsl  snd  tangential  directions  are 


F 
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cos  9  dx'  dr' 


(99) 


(100) 


where: 
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Differentiating  Eqs.  (99)  and  (100)  we  obtain 
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r  __  —  1 

X  3F  .  3f.  .  ,  dr,  -I 

r  .  hS  de  +  —  — 7  de  +  ~  — r  e  d  4 

t  C  ae  a>  ;  «®d 

L  ae  d0  J 


(103) 


substitution  of  Eqa.  (102  and  (103)  Into  Eq.  (94) 
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The  stiffness  and  damping  coefficienta  are  defined  by: 


",  •  '  V '  B»‘  *  W 


K  X  -  B  X  -  K  y  -  B  y 

yx  y*  yy  yy 


To  determine  the  8  coefficienta,  aubatitute  Eq.  (95)  and  (96)  into  Eq.  (104)  and 
collect  the  terms  in  accordance  with  Eq.  (105)  to  get: 


•  f  btr  2  »t  ,  ,  !v 

*  c  ^  bf"  C0,Z  *  +  ^  c0*  *  *in  *  * 

,  r  2  3?t 

•  £  \<o  — 7  cos2  d  +  — »  coa  d  ein  d 

»«> 

.  3?  3F.  1 

„  llfl-i  (—1  coa  4  +  -4  ain  d) 

ad)  ad)  J 


•jL  sin  4 
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where,  in  the  coordinate  system  selected. 


_JL 

jiNLD 


(1X4) 


(115) 


and  all  force,  and  derivatives  are  calculated  for  the  given  steady  state  position, 
defined  by  (€q,  6q). 
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APPENDIX  III  -  COMPUTER  PROGRAM  PH  412 


PERFORMANCE  OF  A  HERRINGBONE  JOURNAL  BEARING 
OtagAtiu  i»  IriE  IuR5uLZnI  REGI.-m 


All  input  data  appears  in  the  form  of  name  list.  A  full  description  of  name 
Hat  can  be  found  in  Ref.  10.  The  input  are  contained  in  the  namelist  "INPUT". 
The  input  data  are : 

1.  REN  ■  Reynolds  number 

■  2 « I »,  - » J  M 

v 


where  N^  and  Nq  are  the  rotational  velocities  of  the  Journal 
and  bearing  in  cycles/sec.,  respectively. 

R  ■  Radius  of  Journal ,  in. 

C  ■  Nominal  radial  clearance,  in. 

2 

v  >»  Kinematic  viscosity,  in  /sec. 


2.  ALOVD  -  Length-diameter  ratio 
-  L/D 

where  L  ■  Length  of  Journal,  in. 

D  -  Diameter  of  Journal,  in. 


3.  COR  ■  Film  clearance  ratio 
-  C/R 


4.  PFIX  »  Dimensionless  gage  pressures  at  the  ends  and  the  middle  of 

the  bearing. 

P  (Pressure.  PSIG) 

U  <!>  (st  +  H0> 

2 

where  n  ■  dynamic  viscosity,  lb-sec/in  . 

The  pressure  at  the  initial  end  appears  first,  then  st  the 
middle  and  finally  at  the  final  end. 

5.  11  ■  Number  of  axial  grid  points  in  the  first  region  of  the  bearing 

counted  from  the  initial  end  of  journal  (see  Fig.  2).  An  odd 
number  is  required  and  also  it  must  be  at  least  2  less  than  12. 
The  minimum  permissible  value  for  11  is  3. 


-107- 


6.  12  -  Number  of  axial  grid  point*  counted  from  the  initial  end  of 

•  •  •  *  •  -  .  •  J  .  ^  «-t _ _  .  _ _ J _ _  i _ _  n _ _  _ - _ ] 

klta  jwuuwii  tu  UHB  cuu  vi  LUC  acvvitv  iwa  a  ««nw«w 

bearing,  12  -eprerents  the  grid  points  of  one  half  of  a  bearing. 
For  a  thre  glon,  non-vented  bearing,  12  representa  the  grid 
pointa  for  tna  initial  end  up  to  the  Interface  of  the  aecond 
and  third  region  (aaa  Fig.  2).  The  minimum  permissible  value 
for  12  ia  5  and  the  maximum  la  17.  Again,  12  must  be  an  odd 
number. 


7.  MORI  »  Indicator  to  a pacify  whether  or  not  there  la  another  set  of 
input  carda  to  follow  tha  present  sat. 

MORI  ■  F,  last  aat  of  input 

MORI  ■  T,  another  aat  of  input  to  follow 


».  SIGN 

9.  N 


10  MOOT 


11.  PHOOT 


"0  "1 
"  »o+Ni 

-  Number  of  grid  points  In  the  circumferential  direction.  The 
maximum  permissible  value  for  N  is  19. 

-  S/2*  <Nt  +  SQ) 

where  t  «  time  rate  change  of  eccentricity  ratio  of  the  journal. 

■  Dimensionless  whirl  velocity  ratio 
4/2*  (Nt+  Nq) 

where  $  ■  journal  whirl  velocity,  rad/sec. 


12 .  OAMD0T  -  7/2*  (N£  +  Nq) 

where  y  -  time  rate  change  of  angular  misalignment. 


13.  IPS  ■  Eccentricity  ratio 
-  e/c 

where  e  »  eccentricity  of  journal,  in. 


14.  GAM  "  Misalignment  of  Journal,  degrees. 

15.  VENT  ■  Indicator  to  specify  whether  or  not  the  middle  of  bearing  is 

vented.  Set  VINT  ■  T  only  for  a  bearing  vented  at  the  middle. 

For  non-vented  bearings,  set  VENT  “  F. 
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16.  BETA 


«  Groove  angle  In  degrees.  Beta  must  be  specified  for  the  first 
region  of  the  bearing  (first  value 1  and  the  second  r#»(nn  of 
the  bearing  (second  value).  If  VENT  ■  F,  the  value  of  BETA 
in  the  third  region  uniat  be  Specified.  For  a  pump-in  design, 
with  grooving  in  the  first  region,  and  «  smooth  seal  in  the 
second  region,  BETA  should  be  read  in  as  an  obtuse  angle, 
the  same  value  for  BETA  being  read  In  for  the  second  region  as 
the  first  region.  For  a  pump-out  design,  with  grooving  in 
the  second  region  of  the  bearing  and  a  seal  in  the  first  region 
BETA  should  be  read  in  as  acute  angle  with  the  same  value  of 
BETA  being  read  in  for  the  first  region  at  for  the  second 
region.  If  VENT  -  F,  the  value  of  BETA  in  the  third  region 
should  be  given.  The  third  value  of  BETA  should  be  a 
conjugate  of  the  first  value  of  BETA,  i.e.  BETA(3)-180°-BBTA(1) 
Never  set  BETA  -  0. 

17.  DBF  •  Groove  recess  ratios  in  two  or  three  regions  of  the  journal 
with  the  first  value  referring  to  the  first  region,  etc. 

-  B/c 

where  &  ■  groove  recess,  in. 

To  impose  the  condition  of  a  smooth  bearing  (no  grooving)  in 
either  region,  sat  DEP  ■  0  for  that  region. 

18.  ALPHA  -  Fractional  groove  width  in  two  or  three  regions  of  the  Journal 
with  the  first  value  referring  to  the  first  region,  etc. 


where  a  and  «  are  the  widths  in  inches  in  the  groove  and 

8  r 

ridge  portion  respectively.  To  Impose  the  conditions  of  a 
smooth  portion  in  either  region,  set  ALPHA  «  0  for  that  region. 

19.  PPOUT  ■  Indicator  to  specify  whether  the  pressure  distribution  in  the 
bearing  is  required  as  a  part  of  output. 

PPOUT  ■  T,  pressure  distribution  is  printed  out;  no  pressure 
is  printed  out  when  PPOUT  ■  F. 
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Output 


1.  Under  the  heeding  of  INPUT,  the  complete  aet  of  input  in  the  namelist  "INPUT" 
IB  printed  out .  seen  quantity  !■  luBui.ifi.su  ujr  un  >am  cyicbcl  iS  used  it! 

the  namelist  "INPUT". 

2.  In  the  case  of  PPOUT  ■  T,  a  heading  of  "Final  Preasure  Diatr ibution"  is 
printed  out.  Below  the  heading,  if  VSNT  ■  T,  there  ere  a  total  of  12  number 
of  lines  of  pressures  with  the  first  line  referring  to  the  initial  end  of 
the  bearing.  If  N  is  less  or  equal  to  10,  each  line  contains  N  number  of 
pressures  starting  at  ^  ■  0  (see  Fig.  1).  For  the  case  of  N  >  10,  the  number 
of  lines  are  double.  If  VINT  -  F,  there  are  a  total  of  (12  +  Il-l)  lines 
when  N  £  10;  The  number  of  lines  will  be  double  for  the  caae  of  N  >  10. 

3.  Regular  Output: 

There  are  a  total  of  nine  quantities  in  one  line  under  the  heading  of 
ESN  NO.,  ICC.,  TORQUI  J.,  TORQU*  B.,  RADIAL  LOAD,  TANG.  LOAD,  FLOW,  COS. 
MOMENT,  SIN  MOMENT,  which  are  defined  below. 

a.  UN  NO. 

b.  IPS 

c.  TORQUE  J. 


d.  TORQUI  B. 


a.  RADIAL  LOAD  ■  Dimensionless  radial  component  of  load 


*  <N1  +  V  <C>  r2 

where  F;  •  radial  component  of  load,  lb. 


•  Same  as  the  input. 

*  Same  as  the  input. 

-  Dimensionless  torque  on  bearing 
T. 

"  wT 

where  Tj  ■  torque  on  Journal,  in-lb. 
V  -  total  load,  lb. 

■  Dimensionless  torque  on  bearing 


*b 

VC 


where  T^  -  torque  on  bearing,  in-lb. 
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f.  TANG.  LOAD  -  Dimensionless  tangential  component  of  load 

Ft 

"  - 

it  /W  4-  W  ^  f—\  V 

-  '  i  o'  'C'  " 

where  Ff  -  tangential  component  of  load,  lb . 

g.  FLOW  -  Dimensionless  flow 

■  “J - ^ - - - 

R  C  (Nt  +  No) 

where  Q  ■  flow,  cu.  in/sec. 

h.  COS.  MOMENT  -  Dimensionless  radial  (cosine)  component  of  moment, 

about  the  initial  end  of  journal. 

M_ 


- r^r 

H  (Nt  +  No>  <§)  R 


where  -  radial  component  of  moment,  lb- in. 


1.  SIN  MOMENT  -  Dimensionless  tangential  (sin)  component  of  moment, 
about  the  initial  end  of  Journal . 

M 


- - 

H  (N,  +  No)  (§)  R 

where  ■  tangential  component  of  moment,  lb-in. 

A  Fortran  listing  of  Program  PN  496  is  provided  in  the  next  few  pages. 
Typical  listinga  of  input  and  output  are  also  given. 
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4 


or»n 


KtNKiNUOWnic  JWwSikAl 

WITH  LARGE  ECCENTRICITY  AND  MISALIGNMENT 
IN  TURBULENCE  REGIME 

COMMON  MDlAG.DELX.DELi.OELiS.M.MS.MGl ,M62 *N .NP 117 .19 ) .AH 2. 5 ,1 
19)»Bi(2* 19 I.CC (2*19)* AF 1(17.19). AF2(17*19J .AF3 ( 17 . 19) »AF4< 17 ,19 ) *A 
2F6(17*19)»AF 7117*19)* AF 5(17*19) .PHI < 17.19) 

DIMENSION  BETA! 3 ) *DEP O ) .ALPHA ( 3 ) *PF1 X( 3 ) »  QQQ ( 19 ) *H9( 17 

1*19) tOQQQI 19 ) *PP(17)*PPP(17) *PPX( 17) .PPPXt I  7 ) , AX1 ( 17 . 19 ) *AX2 ( 1 7. 19 
2)  ,AX3( 17.19) *AX6 Cl 7.19). AX7( 17,19) ,AX6( 17,19) ,AX9( 17,19) ,WS11( 17.1 
39)*  WSM  (17*19) ,WS1 3 ( 17. 19). WS 14(17.19). XX (17) , AXL( 2  «19).CZ(2.19) 
DIMENSION  KUPT(17,19) 

LOGICAL  PPOUT.VENT. TORO. MORE 

NAMELIST/ INPUT/REN, ALOVD.COR.PFIX.il. 12. 

1M0RE.SIGN*  N,  EDOT .PHDOT i&AMDOT *EPS*GAM. VENT .bETA.OfcP. ALPHA  * 

1PP0UT 

2  FORMAT  t  7E14»7 ) 

4  FORMAT (/10C1X.F11.7)) 

9  FORMAT (29H0F INAL  PRESSURE  DISTRIBUTION*  //> 

11  FORMAT (6H1 INPUT  I 
10  READ( 5* INPUT ) 

WR1TE(6*11 I 
WRI TE ( 6  » INPUT ) 

MDIAG-0 

MD*C 

RATLD*2.«AL0VD 
MG1-I1-1 
AMG-MG1 
MS* I 2“I 1 

IF (VENT )  GO  TO  102 
M«12-*ll-l 
BMM1-M-1 
MG2*MG1 
GO  TO  101 
102  M«I2 

RATL0*2*»RATL0 
BMMl*  M-l 
BMM1-BMM1*2. 

MG2«0 

101  DO  30  1  *  1  *M 
DO  30  J“1 »N 
30  NP  ( I  *  J )  *0 

DO  20  J-l.N 
NP( 1 . J ) *1 
20  NP(M. J)*l 
35  KX*1 
100  AN*N 

PI*3. 14159265358979 

DTHETA»2.*PI/AN 

DTHE2-0.5/DTHETA 

TPS«2.»PI*SIGN 

C0N1«6.*TPS 

RADIAN* *01 745 3292 5 19943 
GA1*GAM#RA0I AN 
EDT*0«0 
C0N2-C0N1/REN 
TORQ*. FALSE. 

_TPSl*TPS/8.0  .  _ 

TRQ«TPS1*REN 
OELX-DTHETA 
0EL2-RATLD/BMM1 
TPS1*(2.»P I ) **?  *8  « 
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31 

32 

33 

34 
3E 
36 
3*/ 
36 

35 

43, 

42 

43 
4^ 
45 
4tl 

46, 

4C.I 

5C 

51 

52 
5? 
5<J 
5: 
56 
51 


r' 


I 


i 


v 


AM5»HS  _  _  _ 

DELZS-DELZ 

ANG-0.0  _ 

PHlMOPHDOT-l  l.-SIGN) *0.5 . 

1  ND*1  _ 

ER02*1. 

Ll*l  _ 

L2“  1 1 

38  DO  40  J«1.N  _ 

CC  <  LI  * J ) M0.0 
AXL ( L 1  * J ) *0  « 

CZ(L1» J>*0. 

A1 <L1.2»J)*0.0 
A1(L1*3»J>«0.0 

40  AF7 ( L2 * J ) *0.0  _ _ 

1 F ( IND.EQ.2 )G0  TO  48 
I F ( VENT )  GO  TO  48 
Llr2 
L2“  I  2 
I  NO-2 
GO  TO  38 
48  NR* 1 
MN-MG1 
MM*  1 1 
ISS*1 
OEL-DELZ 
D2*0. 5/DEL2 
13*1 
2*0.0 

200  BET«BETA(NR)*RADIAN 
DEPH-OEP ( NR  1 
ALPH*ALPHA(NR) 

ALM1-1.-ALPH 
ALTAL1*ALPH*ALM1 
SINB-SINCBET ) 

COSB-CQS ( BET ) 

SINB2*SINB*SINB 

COSB2»COSB*COSB 

COT2*COSB2/SINB2 

0LZ»DEL*UTHETA 

I F ( TORO )  GO  TO  1001 

DO  2000  I  *  I S  »MM 

xxm*z 

ZCOR-Z/COR 
EPZGH*EP$+GA1*ZC0R 
EDZR*ED0T+GAMD0T*ZC0R 
EPZG“EP2GH*PH IMC 
DO  2001  J*1.N 
SI*SIN( ANG) 

CO*COS ( ANG ) 

H«l.+CO*EPZGH 
H9( I t J ) *H 
H3*H*H*H 
HG-H+DEPH 
HGHR-HG/H 

AX9 ( I  .J)-(,EDZR*CO+EPZG*SI  )*SINB»P1*24^ 
HG3»HG**3* ' 

HGHR3«HG3/H3 

51- 1./S1NB2 

52- -COSB/S I NB2 

RENR*REN*H  _ 
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“S 

ai 

91 

« 

91 

9; 

9; 

9< 

“51 

9« 

_9I 

91 

IOC 

To: 

10; 

10: 

10* 

10! 

io« 

10‘ 

101 

To- 

IK 

ii' 

n; 

ll: 

IK 

11: 

U- 

ir 

11. 
ii* 
121 
12 

12. 
12 
12- 
T2 
12' 
J2 
12 
T2 

II' 

13 

13 

TJ 

13 

JJ 

13 

T5 


14 

Vi 

14 


I¥ 

14 

7* 


t 

( 


i 


•>*»{.  4*.  u 


R£NG«REN»HG 

GXR«GTCF4RENR)*12< 

GXG«GTCF ( RENG J*12< 

GZR*GZCF  4  RENR ) *1 2  « 

62G«G2CF4RENG)*12< 

S«l.<fiXR*GZR*COT?» 

S6MGXG+GZG_*COT2  I 
~S7"G2G*S2 
_S<L«GZR*S2 
S9-G2G«S1 
§10»GZR*S1 
$11  —  S5/S6/HGHR3 

512- 1.-S11*AIPH-ALPH 

513- 4S7-S8/HGHR3I/S6/S12 
_S14*RENR/HG3/S12*(1<-HGHR>/S6 

ALH3«ALPH*HG3 

A1S«SINB*4ALH3*$6*S11-ALM1*S5*H3)/S12 
A2— ALMi*4ALPH*S5*S13*$6)*H3 


A2-( A? +ALH3* ( ALM1*S6*S 1 3-S7 ) ) »S I NB 

A3«ALTALi*Sl4*4S6*HG3-$5»H3)/REN 

A3-C0N1*  4 A3+H*ALM1+HG»ALPH ) *S I NB 

0MY>H3«ALM1 

S7H«HGHR3*S7 

BIS— H3/S12*<  S8*ALM1-S7H*ALPH*S11 > «SI NB 
B2«ALH3»4$9-S7*$13*ALM1»+DMY*4$10+ALPH*S8*S13) 

B3"C0N1*ALTA11*S14* (S8-S7H) 

B2— B2«S1NB 

AX1 ( I • J) «A1$+A2*C0SB 

B3— B3*H3#SINB/REN 

AX2 ( 1 1 J ) "A2*5INB 

WS11 (I »J)aS 11/512 

WS12  4 I t J )*1 »/S12 

WS134I.JI-S13 

WS144 1 *  J ) «$14 

AX3 4 1  >  J) "A3 

AX64 I»J)-B1S+B2*C0SB 

AX7( l»J)-B2*$INB 

AX8 4 1  *  J ) "93 

IF<M0.NE.2>  GO  TO  2001 
IF  4  I. EG. IS.  AN0<J<EU<1 ) 

I  WRITE(6i2 )  AX 14 I»J) »AX2( I t  J I »AX3  4 I »  J ) « AX6 ( I » J ) * AX7 ( 

II  •J)tAXB41»J)*AX9( I iJt 
2001  ANG-ANG+OTHETA 

Z-2+DEL 
2000  CONTINUE 

IF4NR.EQ.il  GO  TO  2010 
IF<NR<EQ<3I  GO  TO  2014 
IF(VENT)  GO  TO  2014 
GO  TO  2020 
2010  U«1 

G.Q  TO  3800 
2014  IK-MM 

_ _ _ _ - . 

COMPUTE  THE  DIFFERENCES  IN  XS  AND  COEFFICIENT 


3800  I-IK _ 

DO  4010  J«1 »N 


AF1 4 1 » J>  »0<0 
AF24 I »J)*0<0 


MlUi.4L-Q.0__ 
AF44 I • J) «0<0 
MUixJJ-flafl _ 
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149 

ISQ 
1  ?ll 

152 

153 

154 

155 

156 

157 
156 
15? 
16C 
161 
162 
162 

164 

165 

166 
167 
166 
16? 
17C 

171 

172 
172 
174 
17* 

176 

177 
176 
17? 
16C 
181 
182 
18  2 
18* 
16' 
186 
167 
186 
16? 
19C 

191 

192 
192 

194 

195 

196 

197 
196 
19? 
20C 
201 
202 
20: 
20* 
205 
204 
201 
204 

„2P? 


AF6 ( I *J)=1.0 
AF7|  I  ..II--PF  IXtNRl 
4010  CONTINUE 
2020  I  ST* I S+ 1 
MM1*MM-1 

00  4000  I  *  I  ST  *MM1 
10=1-1 
17=1+1 

00  4000  J= 1 *N 
0THsUTHE2 

1 F  C J.EQ.l.OR. J.EQ.N)  GO  TO  4004 
J0= J-l 
J1*J+1 
GO  TO  4008 

4004  IF(J.EQ.l)  GO  TO  4006 
J0*N-1 
Jl  =  l 

GO  TO  4008 
4006  JO=N 
Jl=2 

4008  AF1< J  ,J)=SlNB*AX7< I ,J) 

AF2( ' .  J  1  *  < AX2 ( I ♦ J1 ) -AX2 I  1 » JO ) +COSB* ( AX7 ( 1 tJl )-AX7( 1 .JO)  > >»DTH 
1  +SI NB*  t  AX7( I 7. J 1-AX7 ( 1 0 i J ) ) *DZ 
AF3 ( 1 ,J)=AX2 ( I  * J)+C0S8*AX7 ( I «J)+SINB*AX6l I .J) 

AF4( 1 • J I  * ( AX  1 ( I ,J1 »-AXl ( 1 i JO ) +C0S6* ( AX6 ( I »J1 1-AX6  < I »J0) ) )*DTH 
1  +S1NB*CAX6< I7,J1-AX6( I0.J))*DZ 
AF  5  ( 1 • J 1 =AX 1 (  1 1 J ) +C0Sb*AX6 I  i  .  J ) 

AF6( I  .Jl=0.0 

AF7 ( 1 < J 1 = ( AX3 ( 1 »J1 )-AX3( 1  * JO ) +COSB* ( AX8 ( I » J1 )-AX8( l.JO) > )*DTH 
1  +S1NB*(AX81I7,J1-AX8( 1 0 . J ) I *DZ+AX9 < I .J) 

IF  (MD.EQ.2  .ANO.  J.  EQ.  1» 

1  WRITE16.2)  AFl(  I  .J)  .AF2U.J)  »AF3<  l  .Jl  »AF4<  I.J)  .AF5( 

11 »J) ,AK7< 1 .JJ 
4000  CONTINUE 

IFlNR.EU.il  GO  TO  4020 
IF  INK. EO. 2)  GO  TO  4040 
11=2 

I E= I  2  t 

GO  TO  4090 
4020  11=1 
1  E=MM 

GO  TO  4090 
4040  11=1 
I  E=  I S 

4090  DO  4100  J* 1 *N 
All  1 1 »1.J1=0.0 
A1 ( 1 1 *5  * Jl =0.0 
TRS=AX7« IE  t J  1 
TR$Z=TRS*DZ*2.0 
A1U1.4»JI=-TRSZ 
All  II»2.J)=AUII.2.J)-TRSZ*ER02 
DOT=AX&( 1E.J 1+UTHE2 
CCl 1 1 .J)*CZ< 1  I » J ) —DDT 
C2( 1 1 ♦ J1=0DT 

ami.ji—cccii.ji 

A1(1I.3,J)*A1<I1.3»J)+TRSZ  . 

AF7( IE.J)*AXLI 1 1 .JI-AX81 IE.J) 

AXL (  1 1  »J)  =AX8(  IE.J)  - - 

AF1( IE, J)«0.0 
AF2 ( IE  * J ) =0. 0 

AF3( IE.J) =0.0  _  _ 
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_2id 

21Ii 

212 

212 

214 

215 

214 
217 
21E 

21 5 
22C 
221 
222 
222 
224 
222 
22t 
221 
224 
229 
23C 

231 
2  32 

232 

234 

235 

236 

237 
23E 
239 
24C 

241 

242 
242 

244 

245 
244 
247 
244 
249 
25C 
25] 
252 
252 

254 

255 
254 
257 
254 
259 
26C 
261 
262 
262 

264 

265 
264 

TF3 

264 

27C 


_ _ Af-Ai.L£  0 _ 

AF5I IE. JI-0.0 
AF6  ( IE* JI-0.0 

“ - IF  I  MO  .  K.'2  .'AW5;^TEi57Kn - 

1  WB!te»6*2)  Al(!Is2:J)sAl<ll»4*J). 

1 »  »CC 1 1 l  i J I • AF 7 1 r E.J I 

4100  CONTINUE  _ 

iFINR.EQ.il  GO  tO  4102 
GO  TO  41^0 
4102  IS-tl 

_  NR*2  _  _  _ 

MM*  1 2 
ER02-G. _ 

DEl-DELZS  "  .  .  '  ' 

DZ-0.5/DELZS _ 

GO  TO  200"  ‘  . . 

4160  IFi  VENTI_  GO_TO  4190_ _ 

iFIKK.EQ.il  GO  TO  4180 

IS*  12 _ _ 

MM*M 

_ NR*  3  _ _ : _ _ 

ER02*0. 

_ PEL-DEiZ _  _ _ _ _ _ 

DZ-0.5/DELZ 

_ 60  JO  200 _ _ _  _ _ _ 

4180  IE-12 

_ 0*02*1,  _ _ _ _ _ 

KK-2 

_ 1.1*2. .  ..  ..  . . . 

GO  TO  4090 

JUL50..  JFINR.EQ.2 1  GO  TO  4160  . . 

4190  MTSR*MD 

_ MP._*MT_SR _ _ _ _ _ _ 

CALL  CIC1  CMOS 
M0*Q 

IF (PPOUT 1  WRI TE ( 6.9  I 
_  00  575  1*1. M 

IF  I PPOUT 1  WRITE! 6.4) (PHI  I  I ,J)*J*1.N) 

_  00  575  J*1_»N _ _  _  _ _ 

KUPflf  »J)*0 

.  !F<PHli| .JI.GE.O.OJ  GO  TQ  575  _ 

PHI ( I. J I -0.0 

_ KUPT«1.JI»1 _ _ _ _ _ 

575  CONTINUE 

C _  DIMENSIONLESS  FLOW  Q/C ( N1+N2 1 R*R _ 

IP-I1-1 

A204  AFL5«!?.0xO_ _ _ 

DO  4200  J*1 »N 

_ A1J  jLt£S.*±<LQR±J.tJUlt  Nl  00  TQ  4 US _ _ 

VF1-CPHI I  IP. J+l I -PM I ( IP  * J-l I )*DTHE2 
4210  AFLOW-  (  YF1*AX6I  I  P.JI  I+AX81  1P.J  )+AX7(  IP.  J)»  (PH!  I 
1 I *DZ+AFLOW 

__  1FIHDIAG.EQ.2)  WRITE(  6.2  I  YF1  .AFLOW _ 

GO  TO  4200 

4220  IFIJ.EQ.N)  GO  TO  4230 _ 

YF1*(PHI< I  P.2) -PHI  I IP.N) l*DTHE2 

GO  TO  4210 _ 

A2S0 ,YF1*IPHI I 1P»1 I-PHI C IP.N-l I )*0THE2 

—■  ...'1W  TO  4110 _ 

■HMCPiTINUE 

..  Ji»L  OM*Ari_Oy»DTHETA/12. 


1 1 .5.  J)  »t?l  i  1 1  *  J 


,  . J I -PHI ( I  P—1 » J I 


C  TCiRGuE  DiviGEu  or  "iuaiukhkxk^n/i. 

NR*  1 

TORO*. TRUE. 

TU=0.0 
I  Tl*l 
IT2-11 
TC1*0.0 
GO  TO  200 

1001  00  1000  I=1T1.IT2 
DO  1000  J*1 iN 

I F(KUPT< I *J> .EO. 1 1  GO  TO  1000 
I F ( I.EQ. IT1.0R.I .EQ.IT2)  GO  TO  1800 
0PD2=(PHI ( 1+1 .J) -PHI ( 1-1 ,J) ) /DZ 


33  2 
33* 
33! 
336 
331 
33£ 
33* 
34C 

341 

342 
342 
34* 


AFTR-1.0 

1200  IF( J.EQ.l.OR.J.EQ.NJ  GO  TO  1600 

DPOT* I  PHI 1 1  . J+1J-PHI (  I tJ-1 ) )*OTHE2 
1240  H*H9 ( I »JI 

BQ*H*ALTAL1*DEPH 

AQ=(-BQ*C0SB*W513( I . J ) -H*ALM1*WS12 ( I » J>-ALPH*WS11 < I .J)*(H+DEPH» )> 
1DPDT 

AQ=AQ+BQ*(-C0N2*WS14< 1 .J1-WS13I I . J ) *DPDZ*S I NB ) 

TQ=TQ+AQ*DLZ*0.5*AFTR 

GO  TO  1010 

1600  IF(J.EQ.N)  GO  TO  1610 

DPOT* (PHI ( I .2) -PHI ( I.N) >#DTHE2 
GO  TO  1240 

1610  DPDT=1PHI(I.1)-PHI( UN-11 )*DTHE2 
GO  TO  1240 
1800  AFTR*.5 

I F ( I.EQ.IT2)  GO  TO  1810 

DPDZMPHI  (  I  T  1  + 1  *  J  )  -PH  I  ( IU.JI  J/DELZ 

GO  TO  1200 

181G  DPDZ=( PHI ( IT2.JI-PHI ( 1T2-1.J) I/DELZ 
GO  TO  1200 
1  ('  10  R  E=REN*H 

TC2=TCC(REI*ALM1 
RE*REN*(H+OEPH ! 

T  Cl «  TC1  +  ( TC2+TCC I  RE  1 *ALPH I *TRQ*DLZ*AFTR 
IF  (J.NE.l)  GO  TO  1000 
I FIM0.E0.2.AND.I .EQ.IT1) 

1  WRITEI6.2IAQ.TQ. RE. TCI. TC2 

1000  CONTINUE 

I F(M0.EQ.2 IWRITE(6.2)WSU ( I • J ) .WS12 I  I . J  ) .WS13I I tJ) .WS14( I .J) 
lFlNR.EQ.il  GO  TO  1400 
I F( NR.EQ.2 )  GO  TO  1410 
GO  TO  578 
1400  NK=2 


34! 

346 

341 

34£ 

34* 

35C 

331 
352 
352 
35* 

332 
35< 
351 
35£ 
35* 
36C 

361 

362 
362 
36* 
36‘E 
366 
361 
366 
36* 
37£ 
37] 
372 
372 
37* 
37! 
376 
371 
376 


I T1=I1 
I T2*  12 
GO  TO  200 

1410  I F ( VENT )  GO  TO  578 
NR*3 
IT1-I2 
1  T2*M 
GO  TO  200 
578  TQO*  TQ+TCi 

IF  (ST6M.LT .0.)  TQO*-TQO 
THE*0.0 
DO  580  J*1.N 
QQQ( J I *SI N( THE  > 

QQQQ ( J ) *COS ( THE ) 


37? 

38C 

381 

382 
381 
38* 

- 38! 


386 

18! 

38| 


585  TH£«DTH£TA+  tHE  393 

DO  590  t«l»M  39a 

_P  P  <  I  )  *  0  •  0  39? 

PPP(I)*0.0  39£ 

DO  600  397 

DUH-PHHI.J)  39e 

PP(  I ) -PP (  I^+QQQI JJ*OUM  _  399 

t>uo  pppcn-pppi  iT+oubot  J)*dum  <*oc 

PPt  I)-PPm*DTHETA  _  _  _  401 

'  PPP( 1 ) «PPP( I )*DTHETA'  .  . . . '  402 

PPX  ( 1 )  «PP(  I  )  *XX  (  1  )  _  _  403 

590  PPPXtl  )*PPP< I )*XX(I)  404 

FSIN*SUM(PP.M.DELZ)  .  _  405 

FCO£>SUM(PPP.M,DELZ  )  406 

FMSIN*SUM(PPX*M*DELZ)  _  _  J.07 

FMCOS=SUM<PPPXtM*OEL2)  40f 

FMCOS---FMCOS  40? 

FCOS—FCOS  41 C 

WL0AD»FC0S**2+FSIN«*2  _  _  _  _  _  411 

wl6ad*sqrt(wload>  411 

TQO«TQO/WLOAD  _  __  413 

TQ1 =TQO+EP$*FSI N/WLOAD  41< 

WRITE  (6*6)  __  413 

WRITE  (6*7)  REN.EPS.TQI .TQO.FCOS.FSIN.AFLOW.FMCOS.FMSIN  41£ 

6  FORMAT  (112H  REN  NO.  EPS.  TORQUE  _J.__  TORQUE  B.  RADIAL  LOA  413 

ID  TANG.  LOAD  FLOW  COS.  MOMENT  SIN.  MOMENT  )  41t 

7  FORMAT  IF9.2._lX.F6.3,7(lXE13.6_i)  _  _ 

599  IF (MORE  I  GO  TO  10  42C 

STOP  . 

END  422 
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juh  vcur  :ne  c  i  c  1  ( i 

COMMON  MD  I  AO  ,ut  LX  JLti  .DELZS.M.MS  ,MG1  »MG2  .N.NP(I7,15l.Al(2»5.1 

I  <UH2tlVI  tCC  lz!«  19)  lAF  1  I  n.!1}!  .A‘/[  I7.l<}l  .4^1(1  7.1QI  .4FM  17.101.4 
2F  o(  1  7*  1*  )  .AF  n  1  / .  1)  )  .AF  i  I  1  f  »  19  ]  .PHI  (  1  7  •  19) 

DIMENSION  A(l7.17),b(17.17).C<17.17),D(17.17).F(20*17).AK(17.17>.G 
1N  (  1  7  •  1  7  1  .3  I  2  0 . 1  7  1  .  AF  J  I  (1 7  > .H0<  1  7  «  }  7  )  .RF  (  1  7  I  .00' 17.17.'  .DS117)  :(UL7-i 
21 7  > ,GI <  I  7  .  1  7  > .E ( 1 7 , 1 7) 

AN=1.0/DELX  . . 

NC  =  0 

M=MS+MG1  +  MG2  +  1  .  .  . . 

MN 1 =MG1 + 1 

MN2*MNI+MS  . - . 

DX1=AN*0«  3 

DX2=  AN**2  -  -  •  •  - 

203  DO  203  I=1»M 

00  204  II  =  1  iM  ■  - — . 

E  <  I  «  1 1 >=0. 

204  D<  I . m  =0.  - 

f  1 1 1  n  «o. 

203  D  (  I  *  I  )  =  1  •  - - 

DO  300  J=1.N 

WRITE! 7 ) 1 ( t(  1 1 1  I) »0( I  *  I  I )  »l*l.M> • 1 1«1 »M)  - 

IF  (MD.NE.2)  GO  TO  240 

WRITE  <  6.103  )  AF1 ( 1 ,J) »AF2( 1 »J ) *AF3( 1 »J) .AF4( 1 .J  )  .AF5< 1 *J> .AF61 1 ♦ J- 
1)  »AF7 ( 1  .J1 

WRITE  (6.103 )  - . — - 

1  A1(1.1.J).A1(1.2.0).A1(1.3.J).A1(1.4.J) *A1 ( 1 .5.J) .61  (MN1. 

10)  .CC ( MN 1 » J I  - 

WRITE  (6.103) 

1  A1(2,1,J).A1<2.2.4).A1<2.3.J).A1(2.4.J>.A1<2.5.0>»  - 

IB  1 ( MN2  *  J ) »CC ( MN2 » J ) 

240  DZ1*0.3/(DEL2I  --  . 

DZ2=1 ./  (l)ELZ  )**2 

241  D0242  I =  1.M 
DO  242  1  1  =  1  .M 

B(  I  ,  1  I ) =0.  -  - 

A(  I  .  I  1  )  =0. 

242  C  < I .1  I ) =0. 

DO  230  1=1. M 

2U6  I F ( I —MG  1 - 1 )  212.210.212  . . 

209  AF J I ( I ) =-AF7 ( I . J I 
A  (  I  .  I  1*1. 

GO  TO  230 

210  IF (M5)  211  .212.211  - 

211  DZl*u.5/DELZS 

DZ2= 1 » / ( DEIZS ! **2  - - - 

IF (NP< I .J)  1209.232 .209 

232  lF(MGl)  233.212.233 

233  IFINC.EQ.2)  GO  TO  600 

A( 1 . 1-2 )=A1( 1  .1  .J)  . . 

A ( I.I-l ) =A 1 ( 1.2. J) 

A( 1. 1 )=A1( 1.3.J I  -  - 

A (  I  .  1  +  1  )«A1 ( 1 .4, J) 

A<  I  .1+2 )=A1( 1.5.J)  -  -  -  -  - 

6(1.1 )>B1( 1. J) 

CII.I)-CC(l.J)  - 

600  AF J I ( 1 )=-AFTl I.J) 

GO  TO  250  -  -  ■  -■  - . . 

212  IF(NP( ' »J> >209.234.209 

234  I  F ( I -MG  1-MS- 1 )  215.213.215  -  - 

213  IF(MS)  214.215.214 


♦ 


2  14 

0 Z2«  (  1  .  /DEI  2  I  •*? 


63 
..  ** 


IF ( NC • EU. 2  I  GO  TO  601 
A(  ltl*2)>Al( 2  *  1  *  J  ) 

A (  I . 1-1  ) »A1  4  ?  »?,.n 
A ( l *  I '  -A1 ( 2  *  3  * J ) 

A  < 1 . l  +  l >=A1 { 2.4.J) 

A ( 1 . I +2 )«A1 ( 2 ,6, J ) 

B  (  1  •  I  >«ih(2.J) 

C  (I  .1 >  *CC ( 2 ♦ J l 

601  AFJI ( I )=-AF7d  .J) 

GO  TO  260 

215  1FINC.EC1.2I  GO  TO  602 

B  (  1  .  1  -  1  >=AT  3(1. -J  1*0X1*  02  1 
B (  I • I  I = AK  6 (  1  . J ) *OX.  -AF4(  I  . J ) *0X1 

B( I .1  +  1  I *-AF  3 ( 1 ,J 1*0X1*021 

A ( 1 1 1-1  I  *AF1(  I  .  J  1  *02  2-AF  2 ( 1 . J1 *01 1 

A(l.  I)*AF6<  I  .  J  1  -2  «0*  <  AF 1  (  I  .  J  1  *022  +  AF5d  .J  1  *0X2  1 

A  (  I  . t+1 )=AF1 < 1 . J 1 *D22  +  AF 2 ( I »J1*DZ1 

C  (  1  •  I- l  I *-AF 3 ( 1 »J 1*0X1 *021 

C  d  .  J  )«AF4< I *J)*DX1+AF5( I «J 1*0X2 

C (1  •  1  +  1 1  * AF  3 ( 1 .J 1 *0X1*021 

602  AFJI (1 J--AF7I1 .J) 

250  CONTINUE 

I F ( NC.EQ.2 1  GO  TO  603 
00  260  I *1 *M 
DO  260  II-l.M 
AXII.IIdAd.il) 

DO  2 60  1 1  1*1  .M 

260  AX(  I  .1  I  )*AX(  1  .  I  I  1+Bd  •  II  I  )*Sd  II  .  II  1 
IF (M0IAG.LT. 21  GO  TO  262 

WRI TE ( 6 .101  I J 
DO  261  1-l.M 

261  WR I TE (6  •  100 1  ( AX (1 . 1 1  1  *  I  I « 1 .Ml 

262  CALL  MATINVI AX.M.0UM.0.0UM1 1 

603  CONTINUE 

DO  404  l-l.M 
BFd  1  -AFJI  (  1  1 
F  < J+l.I )«0. 

00  404  II-l.M 
E(1  .1 1 1-0. 

BDd.IU-O., 

DO  403  I  I  1-1  .M 

Ed.dl»Ed*4I)-AXd.UI  l*CdIl»i  11 

403  BD(  I  .1 1  1-BO  I  I.  1 1  1+Bd  .11 1  1*0(1  II  .  II  ) 

404  BF(I)-BFlll-B4l*lI.)*ftJ.in 
DO  406  I-l.M 

DO  406  II-l.M  . . . 

DII.1D-0. 

DO  405  111*1 *M 

405  D(  I  .1  1 1-0(1  .  ID-AKt  1.1  II  )*B0(  1 1 1  .  II  1 

AO*  FLJ*l^U*fl4*4.*U^A)U-l-dl-l*B£UIl _ 

300  CONTINUE 

- 00  5051 ■ l.M - 

DO  504  II-l.M 

SOA  00(1. Ill «-JH,UUl _ 

905  DD(l»I)-l. +00(1.1) 

- LFiMOlAfi-LT *2 1  GO— LQ  2*4 _ 

WRITE  (6.101) 

r.  -4KT241  1*1.M - - - - . 

mJWte  i#**oo>  iooii*ii i.u-i.Mi 


65 

66 

6  7 
60 

69 

70 

71 

72 

7  3 

74 

75 
lb 
7  7 
la 
79 
BO 
6 1 
62 

63 

64 
ob 
66 
67 
60 
09 

90 

91 

92 

93 

94 


95 

96 

97 


9B 

99 

100 

101 

102 

103 

104 

105 

106 
107 
100 
109 

no: 

ill' 

112! 
1 1 3  j 

1141 

115 

-  •  11* 
117 
-  110. 
119 

121 

_  ...  122 

123  j 


110- 


I 

I 


i 


call  ha r  > v :  o'; . a , ^ . o « ./um » 

DO  bo  7  I  *  1  * M 

s »  n  .  1 1 * o . 

on  Mt  ?  ]  i  »  i  - m 
ON  Civil  1=0. 

DO  506  III*] *M 

GN ( I .  1  J  ) =  GN  (  I  t  1  1  ) +  DD  ( I  •  1 1  I  1 *E( I  I  1  .  I  I  I 
b  ^0  O  (  1  *  I  1  J  = ON  4  X  *  I  1  I 

bo  7  0 (N, 1  ) =  5<  N« I  )*l>D(  J . 1 1 ) *f I N« 1 • 1  1  t 
J  -  N+  X 

DO  bl2  K*2»N 

WRI TE (8 1  (  Ui  <  i » I  I  i  . I  =  1 .Ml  . I  I  -  1 .MJ 
IJACKGPAO.  7 

KLADt  7  )  (  (L  (  I  ,  I  I  )  ,01  1  .  1  I)  ,  J  -1  ..<i)  .  II  -1  .Ml 

BACKSPACE  7 

I  I-  ( Mo  lAo.LG.2)  wRl  It.  (6  •  100  1  (  (  L  (  J  .  I  I  )  .u(  I  ti  I  1  •  1-i  *!'■)  *  I  1  -l  .Ml 
J  =  J-1 

DO  b09  1=1  •  M 
SlJ-ltJ  )*f  (J.l) 

00  b J V  I  I  «■  i  ,M 
G I C I* 1  I ) =0. 

DO  boa  I  1  I = 1  . M 

Ol  (  I  »  1  1  I  =o  J  l  I  .  1  I  I  -MX  1  , 1  J  I)  *oNt  J  1  1  .  1  1  I  +L  <  I  .  1  1  1  1  *u(  1  I  I  .  1  I  I 
If  (  MIX  1  AO.Nt  ♦  2  )  00  TO  bU8 
bob  CUNTINUL 

bU9  S(  J-l  •  I  IHI  I  .  1 l >*$(  J.  J  1 )  ♦D(  i  .  1  i  I *6  IN. I  1) 

00  bl2  I  -1  *M 
DO  b 1 2  I  1  =  1  «  M 
bl?  G< I  .!  1  l#Ol ( I  till 
<760  DO  bl  1  1  =  1  .M 

DO  b 1 0  I  I  *  1  .  V 
b 1 0  UU(  1  .  1  I  1 -  —  o 1  (  1  .  I  I  I 

b  1 1  uoil.llH.  ♦no  (  1  ,  l  i 

I I  (  Mi)  j  Au.  I.  1  •  2  I  «. j i>  TO  2oo 
WKI  Ti  (6.101) 

no  M>  j  i  - 1 . 

t  6  *'  a'K’  1  T  t  ( (i .  1  .1.  *  (  Oi'  II. 11). 11  =  1.^) 

<.’06  CALI.  MA  1  UnV  (  I '0.  M.OOM.  0,  Oil'll  I 
I  f  (  Ml  >  1  A  u  «  L  I  ./  )  oO  TO  270 
1 

•vKlTL  (O.lUl)i 
DO  27 i  I  =  1  * M 

2  7  J  WRITE  (O.lOo)  (DD(  I .1  1  I ♦ J l *1 .Ml 
2  7u  DO  bib  1  =  1  .M 
PHI  ( I  ,  1 ) =0 
DO  bib  I  1  =  1  .M 

bib  Pul  (  1  .  11=00  (1. 111*5(1.  Hl.PMlll.il 
DO  bl6  J  =  2  »N 
BACK  SPACE  6 

kLAUttt )  < (G( I  .11) . 1 =1 .Ml .1 1  =  1 .Ml 

BACKSPACE  a 
DO  516  1=1. M 
PHI(I.J)=S( J.ll 
DO  516  11=1. M 

516  PHI ( l . J I *PH I  I  I . J)+G( 1 .1  I  l*PHI ( I  1  .11 
REWIND  7 
REWIND  S'** 

IE  (MD.tU.2l  GO  TO  110 
IFIMDI AG.LT.2 )  GO  TO  268 
110  WRITE  <6.1021 
DO  267  i*l. M 


124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 
13b 
1  36 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 
15b 

156 

157 
150 

159 

160 
161 
162 

163 

164 

165 

166 
167 
160 

169 

170 

171 

172 
1  !  3 

174 

175 

176 

177 

178 

179 

180 
184 
182 

183 

184 


1 
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267  WHITE  lfc,109)  < PH  I {  I . J > * J« 1 • N I 
a6*-HS4*X»* 

100  FORMAT ( 1X*1P10E11«4) 

101  -FORMAT <  6MOC1CU 16)  - 

10?  FORMAT  t  1  HO ,10X1 OHF I NAl  PHI  / 1  HO ) 

!0»  FORMAT  UX,1P7E11.4*  - 

END 


165 

166 
167 
166 
169 

190 

191 


112 
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SUBROUTINE  MATINVI A.N.b.M.DETEK ) 

C  MATRIX  INVERSION  WITH  ACCOMPANYING  SOLUTION  OF  LINEAR  EQUATIONS 

DIMENSION  IrIVOii7i*MtxTiiri#t3{x^#i  (.INDEX (17. 2  )«PIV0T(17) 
EQUIVALENCE  ( 1  ROW  »  JROW ) «  (ICOLU  *JCOLU  )*  IAMAX.  T.  SWAP) 
INITIALIZATION  _ 

10  DETER  *1.0 
15  DO  lu  J«1*N 
20  IPIVO  (J)-0 
30  DO  550  I  *  1  * N 

SEARCH  FOR  PIVOT  ELEMENT 

A0  AMAX-O.C 
A 5  DO  105  J  =  1 » N 

50  IF  (IPIVO  (J)-l)  60.  105.  60 
60  DO  100  K-l.N 

70  IF  (IPIVO  (K)-l)  80.  100.  7A0 
80  IF  (ADS  (AMAX)-ABS  (AIJ.K)))  85.  100.  100 
85  I ROWs J 
90  ICOLU  -K 
95  AMAX-AU.KI 
100  CONTINUE 
105  CONTINUE 

110  IPIVO  (ICOLU  )=IPIVC  (ICOLU  )+l 

INTERCHANGE  R0W5  TO  PUT  PIVOT  ELEMENT  ON  DIAGONAL 

130  IF  (IROW-ICOLU  )  1 A 0 »  260.  140 
1 AO  DETER  --DETER 
150  DO  200  L*1.N 
160  SWAP-A (IROW.L  > 

170  A ( I  ROW .L )®A( ICOLU  .L) 

200  AdCOLU  . L )  -SWAP 
2o5  IF(M)  260.  260.  210 
210  DO  250  L  =  1 .  M 
220  SWAP=B< IROW.L) 

230  BHR0W.L)*B(  ICOLU  .L) 

250  B (ICOLU  .L) -SWAP 
260  INDEX!  I . 1 ) *  I  ROW 
270  INDEX! J.21-ICOLU 
310  PIVOT! I )*A( ICOLU  .ICOLU  I 
320  DETER  -DETER  *PIV0Td) 

DIVIDE  PIVOT  ROW  BY  PIVOT  ELEMENT 


330  AdCOLU  .ICOLU  >*1.0 
340  DO  350  L-l.N 

350  AdCOLU  »L  )  -A  ( ICOLU  .Ll/PlVOTil) 
355  IFIM)  380*  380.  360 
360  DO  370  L-l.M 

370  B  ( ICOLU  *L  )  -B  (  ICOLU  .Ll/PIVOTU) 

REDUCE  NON-PIVOT  ROWS 

380“D0  550  Llr-l  »N 

390  IE(L1-IC0LU  I  400.  550.  400 

400  T -A  (LI  .ICOLU  >  . . " 

420  A(L1» ICOLU  1-0.0 
430' DO  450'  L-l.N 

450  A  <  tJjJLL*  A  ( L 1 .  JL> -A  (ICOLU  .L  j_*T 


J 

— -i 
i 

t 

M 

6 
' 

l 
1C 

"n 
1 2 
i: 

u 
1; 
_i_< 

li 

JE 

is 

2C 
21 
22 
~21 
2A 
2£ 
2£ 
~21 
2t 

"2l 

_ M 

31 

_ 32 

3 1 

_ _ 34 

35 

_ 36 

37 

_ 3€ 

3S 

_ 4C 

41 
_ 42 

42 
_ 44 

45 

_ 46 

47 

_ 4f 

4S 

_ 5C 

51 

_ 52 

55 

_ 5£ 

55 

_ 56 

57 

5E 

- 5s 

6C 

_ 
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*35  1 F  ( M t  s«>n.  s’-.n.  iaa 
*60  DO  500  L«L»H 

500  BIU1»L)-B(L1»LI-B( ICOLU  »L)*T 
550  CONfTNyE  "~ 

C  _  _ _ _ _ 

c . Interchange  columns 

c _  _  _ 

600  DO  716  r-l.N 

6_10  L*N+lti  v  _  _ _ 

620'"Tf"  (INDEfXIL.l  )-IND£X(Lt2> )  630.  fToT'630 
630  JROW* I NDEX ( L . 1 ) 

6*0  JCOLU  *  INDEX! L*2 ) 

650  DO  705  K«1.N 
660  5 WAP* A ( K . JROW ) 

670  A(K.JROW)^A(K.JCOLUJ  _ _ 

700  A  (it.  JCOLU  y-SWAP 
705  CONTINUE 
710  CONTINUE 

_ 7*0  RETURN 

END  . 


63 

bH 

63 

St 

6  7 
6E 
6S 
"7Q 
71 
7 ; 

7 

In 


V: 

It 

77 


7  fc 
75 


Jd 
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■  viv  i  ivn  ivv  me  i 

IFire. gt. 100.01  go  to  10 

TC08./RE 
GO  70  100 

10  IFire. GT. 400.0)  GO  TO  2 0 
TC04.175/(R£)#*.86 
GO  TO  100 

20  IFire. GT. 1000.01  GO  TO  30 
TCC=»547/(SE) **• 521 
GO  TO  100 

30  IFIRE. GT. 4000.0)  GO  TO  40 
TCC=.342/(RE)**.4b3 
GO  TO  100 

40  TCC=. 064/IRE)**. 25 
100  RETURN 
END 


4 


e 

~1 

E 

5 

1C 

Tl 


14 

is 

it 

"n 
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function. sum tp.w. ox i 
DIMENSION  Pll'M 
_  K-2  _ 

KX-M-1 

. . .«*■* 

SUM-0.0 

10  DO  2J>  Jj>K*KK»KKK 

20  SUM-SUM+P ( it 

60  _T0  130.40 _j50  )  .(C _ 

30  SUM ■ SUM* DX / 3 • C 
RETURN 
40  K»3 

45  SUM-SUM*2.0 
GO  TO  10 
50  K  ■  1 
kk-m" 

KKK-M-1 
60  TO  45 
END 


FUNCTION  GZCF(REYN) 

RE=REYN 

IFIRE. LE. 70. 0)GQ  TO  20 

I K ( (RE. GT. 70.0). AND. (RE. LE. 4000.0)100  TO  30 
10  IF(  (RE. GT. 4000.0 ) .AND. I  RE . L E . 7 . 0E  +  03 )  I  00  TO  AO 
I F ( (RE. OT. 7000.0). AND. (RE. LE. 2*0 E+OA) )  GO  TO  50 
GZCF=25.6/(RE)**.756 
RETURN 

20  GZCF=1. 0/12.0 

RETURN  . . . . . 

30  GZCF=1.85B  E -09* ( RE ) **2- 1 . 8 76 E-05*RE+ . 06A6 
RETURN 

AO  GZCF-9.62/ (RE )**.652 
RETURN 

50  GZCF  =  1 1 .3/ 1  RE  I »*.67A 
RETURN 
END 


FUNCTION  GTCF(REYN)  2 

RE-REYN  j 

IF (Kt.Lt.70.0100  TO  20  l 

IF  I IRE. GT. 70. 01. AND. (RE. LE. 2000.0) )  GO  TO  30  £ 

10  IFIfRE.GT. 2000.0)  .AND.  (RE. LE»5.5£+03))GO  TO  40  f 

IF( (RE.GT.5500.0 I .AND. (RE.LE.2.0E+04) ■  GO  TO  50  "  *i 

GJCF-20.5/IRE)  **0.784  t 

RETURN  S 

20  GTCF-1. 0/12.0  1C 

"  RETURN  .  Ti 

30  GTCF-.619E-08*(RE)**2-3.465E-05*RE+. 08569  12 

RETURN  O 

40  GTCF-4.90/(REl**.626  U 

RETURN  1* 

50  GTCF-10.35/ ( RE)**.716  It 

RETURN  IT 

END  IE 


...  _ ...  . . -f- _ 

1  0.1400000IE  09  9.31QOOOOOT  i)2  0.310900001  02 

_ Hi! _ LU _ _ _ 

1  0.  0. 24000000* 01  0. 

4L0M4  »l>  _  _  _ 

i  ■  "o'.  '"'(i.'soooooofff  "n  o.ssiSoffcTrsnfs 

MOUT _ t _ _ _ _ 

YON*  0. 


fjHll  fPg»fUJO|  0  [ST4|8UTiON. 


.*42770}  0.71002*0  0 , 7400000  I.DHIII  0,70240917  0,9023(74  (.77*0714  0.7400201  Q,  >110111 

0. *037000  -0. *702090 


1.0121377  1.0047*4*  1,10*30*]  1.1013404  1.1*1*010  1.10*209*  1.1724*00  1.13*70*0  1.0*2*132 


0.3410142  0.341*344  0.3**2073  0.3*1907*  0.391712*  0.1*9*441  0.3*493*9  0.30149*9  1.1(9101*  (.17*09*4 


10.00  0.090  0.19771*0  03  0.l799l3£  03  0.l79 


09*11  01  0.101379E  00  0.9701(((-(2 


COMPUTER  PROGRAM  PN  406 


STATIC  PERFORMANCE  OF  A  SPIRAL -GROOVED.  FLOATING -RING  JOURNAL 
BEARING  OPERATED  IN  THE  TURBULENT  REGIME 

Input 

All  Input  data  appear  in  the  form  of  name  list  except  when  the  characteristics 
of  inner  and  outer  film  are  known  and  need  not  be  computed  within  the  program. 
In  that  case,  two  sets  of  data,  each  containing  9  cards  in  a  specified  form 
as  explained  in  detail  below,  are  required  to  provide  the  information  on  the 
film  characteristics. 

For  the  readers  who  would  like  to  be  familiar  with  the  format  of  namelist,  it 
is  recommended  that  he  read  pages  14  and  19  of  Ref.  (10)  . 

The  choice  of  whether  to  provide  or  to  compute  the  film  characteristics  is 
Indicated  by  the  first  word  of  the  namelist  "NGPUT".  The  preparation  of  input 
for  each  case  is  shown  below: 

Case  I:  The  film  data  generated  within  the  program  - 

the  namelist  contains  two  listings;  these  are  "NGPUT"  and  "INPUT". 

A.  "NGPUT"  includes  the  following  input: 


1. 

INPRD, 

INPRD  ^  1:  The  film  characteristics  will  be  generated 

2. 

NEPS'= 

within  the  program  and  the  information  in  the  namelist 
"INPUT"  must  be  provided. 

Number  of  (inner  film)  eccentricities  ratios  to  be 

i 

3. 

NCASE, 

examined  (maximum  10) . 

NCASE  -  0:  Last  set  of  input. 

2 

4. 

NN  - 

NCASE  £  0:  More  input  follows,  starting  from  the  namelisting 

"NGPUT". 

Number  of  iterations  to  be  allowed  to  achieve  an 

r 

« 

k 

% 

I 

i 

5. 

R2R1  - 

equilibrium  condition  in  the  floating-ring  system 

(recommend  10) . 

Radius  ratio  of  ring  and  journal. 

r 

i 

“ 

R^/R^  Rj  *  radius  of  journal,  in.;  R^  ”  radius  of  ring,  in. 

# 

6. 

ANTT  =» 

Initial  guess,  on  the  speed  ratio; 

u>2  “  speed  of  ring,  rad/sec.;  «  speed  of  journal,  rad/sec. 
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7.  DN  -  Incremental  value  of  the  ring  speed  ratio  during  the 

iteration  (recommend.  .15) . 

8.  KEY  -  Overall  Reynolds  number  under  which  the  bearing  la  to  be  operated. 

m  — — ’  C.  “  inner  film  thickness,  in., 

*  2 

v  “  kinematic  viscosity  in  /sec. 

9. EPIS  •  Eccentricity  ratios  of  inner  film  to  be  examined.  There  are 

NEPS  number  of  eccentricity  ratios  to  be  provided  (maximum  10) . 

10. COR1  -  Inner  film  clearance  ratio 

■  CA 

11. C0R2  ■  Outer  film  clearance  ratio 


-  V«2 

12. FEES  -  Overall  dimensionless  gauge  pressure  measures  at  the  end  and 
the  middle  of  the  bearing.  . 

-  2«F  (pressure,  psig)  / 


where  4  is  the  dynamic  viscosity. 


The  pressure  at  the  end  appears  first. 


B.  "INPUT"  includes  the  following  input: 

1 . . 'RIM ST  ■  The  starting  value  of  ring-speed  ratio  under  which  the  film 

data  will  be  generated.  It  is  noted  that  0  <  RINSP  <1.  Since 
the  film  data  of  each  film  covers  three  different  speeds,  a 
small  number  for  RINSP,  say  .25,  is  recommended. 

2.  DELS?  -  The  incremental  ring-speed  ratio,  DELSP>  should  be 

sufficiently  small  such  that  RINSP  +  2  x  DELSP  <1.  The 
recommended  value  for  CELS?  is  0.15. 

3.  BLOVD  -  Length  diameter  ratio,  L/2Rlt  where  L  is  the  length  of  journal. 

4.  II  -  Number  of  axial  grid  points  in  the  first  region  of  the  bearing 

counted  from  the  initial  end  of  journal  (see  Pig.  2)  .  An  odd 
number  is  required  and  also  it  must  be  at  least  2  less  than  12. 
The  minimum  permissible  value  for  II  is  3. 
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®-  5.  12  ■  Number  of  axial  grid  points  counted  from  the  initial  end  of 

the  journal  to  the  middle  of  the  journa  1  •  lilt:  uiuuiuuM 
permissible  value  for  12  is  5  and  the  maximum  is  17.  Again, 

12  must  be  an  odd  number.  The  number  of  grid  points  in  the 
second  region  is  I2-I1+1  which  gives  (12-11)  intervals. 

6.  N9  ■  Number  of  grid  points  in  the  circumferential  direction, 

maximum  permissible  value  for  N9  is  19. 

7.  EPS1  ■  An  array  of  three  values  of  for  which  inner  bearing  film 

data  are  to  be  generated.  The  range  of  EFS1  should  be  wide 
enough  to  cover  anticipated  operating  eccentricities. 

8.  BETA  *  Groove  angle  in  degrees.  BETA  must  be  specified  for  the  first 

region  of  the  bearing  (first  value)  and  the  second  region  of 
the  bearing  (second  value) .  For  a  pump-in  design,  with  grooving 
in  the  first  region,  and  a  smooth  seal  in  the  second  region, 

BETA  should  be  read  in  as  an  obtuse  angle,  the  same  value  fcrr 
BETA  being  read  in  for  the  second  region  ae  for  the  first 
region.  For  a  pump  out  design,  with  grooving  in  the  second 
region  of  the  bearing  and  a  seal  in  the  first  region,  BETA 
should  be  read  in  as  acute  angle  with  the  aame  value  of  BETA 
being  read  in  for  the  first  region  as  for  the  second  region. 
Never  set  BETA  ■  0. 

9.  DEP  -  Groove  recess  ratios  in  the  two  regions  of  Journal  with  the 

first  value  referring  to  the  first  region,  etc. 

»  ft/c ,  where  5  ■  groove  recess,  in. 
and  c  ■  nominal,  radial  clearance,  in. 

To  impose  the  condition  of  a  smooth  beering  (no  grooving)  in 
either  region,  set  DEP  -  0  for  that  region. 

10.  ALPHA-Fractlonai  groove  width  in  the  two  regions  of  a  journal  with 

the  first  value  referring  to  the  first  region  ■  a  /  (a  +  a  ) , 

g  g  r 

where  a  and  a  are  the  widths  in  Inches  in  the  groove  and 
ridge  portion  respectively.  To  Impose  the  condition  of  a 
smooth  bearing  in  either  region,  set  ALPHA  *  0  for  that  region. 

11.  EPS2*  An  array  of  three  values  of  ^  for  which  outer  bearing  film 

data  are  to  be  generated.  The  range  of  EPS2  should  be  wide 
enough  to  cover  anticipated  operating  eccentricity. 
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Case  II:  In  the  case  of  known  film  characteristic*,  only  the  namelistlng 

"NGPUT"  la  required.  In  addition  to  that,  IB  caras  must  be  Tullowsd 
which  contain  the  film  data. 

The  content  of  namelisting  "NGPUT"  la  essentially  the  tame  as  those 
of  Case  I.  provided  INPRD  must  be  equal  to  1. 

The  inner  film  data  appear  first  corresponding  to  three  eccentricity 
ratios  al  three  different  Inner  film  Reynolds  numbers. 

The  data  Include  (these  symbols  are  defined  below) 

«.  TJt  Tb,  8,  * 

which  are  punched  on  one  card  with  the  format  of  (IX,  F8.1,  FS.2, 
4E13.6).The  first  three  cards  are  for  the  same  Reynolds  number,  the 
smallest  of  the  three  valuea.  Each  card  refers  to  a  different 
eccentricity  ratio.  Again,  the  order  of  the  eccentricity  ratio  is 
ascending.  The  next  three  cards  correspond  to  a  higher  Reynolds 
number  but  with  the  same  set  of  eccentricity  ratios.  In  total, 
there  are  9  cards  for  the  inner  film. 

Following  the  inner  film  data,  there  are  9  cards  for  outer  film 
data  arranged  in  an  order  similar  to  that  for  the  inner  film. 

The  six  quantities  used  as  input  for  the  inner  and  outer  film 
are  defined  below. 

1.  Re  *  Reynolds  number 

-  ■  -  fr  -—  (for  inner  film) 

“  (for  outer  film) 

v 

2.  €  •  Eccentricity  (see  Fig.  3) 

■  ej/c^  (for  Inner  film) 

■  e^/c^  ( f°r  outer  film) 

3.  T,  ■  Dimensionless  Torque  of  Journal 
(for  inner  fllm)}  where  V  -  load,  pounds 

(for  outer  film) 
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Case  II: 


4.  T 


Dimensionless  Torque  of  Bearing 


(for  Inner  film) 


(for  outer  film) 


5.  S  “  Sommerfeld  Number 

Rj  2 


(for  inner  film) 


6.  t 


Rj  2 

“2“»2  <cj> 


(for  outer  film) 


where  N^  end  N^  ere  speed  of  journal  and  ring  in  rev. /sec., 

.  ■  ,  ,  '  ..  respectively. 

■  Attitude  angle,  degree  (see  Pig.  3) 

■  ^  (for  inner  film) 

■  (for  outer  film) 


Output 

The  output  appears  under  the  title  of  "Floating-Ring  with  herringbone  journal". 

1.  Main  program  input:  it  prints  out  the  title  of  "Main  program  input". 
Immediately,  there  follows  the  title  of  "Namelist  MGfUT"  and  the  entire 
contents  in  that  namelist.  At  the  end,  it  prints  out  "end  namelist 
NGPUT" . 

2.  Input  for  subroutine  '"HERNB":  it  prints  out  the  title  of  "Herringbone 
Bearing  Input"  and  a  title  of  "Namelist  INPUT".  Then,  the  entire  contents 
of  that  namelist  are  printed  out.  At  the  end,  it  prints  out  "end  namelist 
INPUT". 

3.  Single  film  data 

A.  For  Inner  Film: 

There  are  three  eccentricity  ratios  at  three  different  Reynolds 
numbers  corresponding  to  three  different  speed  ratios.  The  output 
starts  with  s  title  of  "inner  film  data"  and  then  the  headings 
REYNOLDS  NO.,  ECCENTRICITY,  INNER  TORQUE,  OUTER  TORQUE,  SOMMERFELD 


NO.,  ATT.  ANGLE,  FLOW,  on  one  line.  Immediately,  there  follows 
9  lines  of  data.  Each  line  contains  seven  quantities  under  the 
appropriate  heading.  The  first  six  of  these  seven  quantities  are 
defined  above  in  the  input  list  for  NGPUT,  Case  II.  The  flow  is 
defined  as 


B.  For  Outer  Film: 

After  the  output  of  the  inner  film,  there  are  a  set  of  output  referring 
to  the  outer  film  Just  like  those  for  the  inner  film.  The  output  com* 
prises  the  title  of  "Outer  Film  Data",  the  heading  and  9  lines  of  data. 


Final  performance  characteristics  of  the  floating  ring  bearing  at  the  steady 
atate  equilibrium  condition.  Under  the  title  of  output,  there  are  four  values 
in  a  line.  These  are: 

ECCENIK/Cj'  ■  e/C^  -  the  overall  eccentricity  of  the  journal  at 

equilibrium  position  divided  by  the  nominal 
clearance  of  the  inner  film  (see  Fig.  4). 

»  C2/Cl  and  R^/R^  aa  previously  explained. 

Next,  there  is  a  table  of  output  referring  to  the  inner  film,  the  outer  film, 
and  the  overall  bearing.  There  are  seven  values  in  a  line  and  a  total  of  3 
lines.  Each  line  contains 

(0)  -®  )R  0 

REYNOLDS  NO.  -  -* — k — k-L-  (for  inner  film) 


13b 


(for  outer  film) 


»,K,c, 


"1V1 


(for  overall) 


ECCENTRICITY  -  e^^ 


(for  inner  film)  -  e, 


(see  figs.  3  U.iU  J 


62  ’VS 


(for  outer  film) 


€  -  e/C^ 


(for  overall) 


TORQUE  ■  Dimenaionleaa  form  defined  aa 


^  |  "  2  ^l^l^tl  (fo^  inner  film) 


[ 


WcT  +2S2e2?t2  (for  outer  film) 


Ll 

W\ 


(for  overall) 


where  Ffcl  and  Ft2  are  defined  on  page  139. 


SUPPLY  IRES  “  Dimenaionleaa  aupply  preaaure 

P 


h 

(nx+n2) 


(for  inner  film) 
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uN 


2 


(for  outer  film) 


(for  overall) 


(for  outer  film) 


(for  overall) 


ATI.  ANGLE  ■  Attitude  angle  in  deg.  (see  Fig.  3) 

■  0^  (for  inner  film) 

»  (f°r  outer  film) 

0  (for  overall) 
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TANG.  FORCE  “  Dimensionless  form  of  tangential  force 

*J| 

rr 

i— ■ 

Ftl 

(for  inner  film) 

WS1 

F 

t2 

Ft2 

(for  outer  film) 

ws2 

where 

■  Tangential 

force  of  the  inner 

film,  pound 

F 

t2  "  Tangential 

force  of  the  outer 

film,  pound 

A  Fortran  Hating  of  program  FN  406  la  provided  in  the  next  few  pages, 
listings  of  input  and  output  are  also  given. 


Typical 
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0370*7*9 


-  V.  v  i  «  i  ui 


■  »ruru9 


—  te*  in 


FLOAT INS-P !NG  IN  CONJUNCT f ON  WITH  HERR INGBONF  JCLRNAl  IN  TMtBULENT 


DIMENSION  P(6, 3, 51,0(6, 3, 3)  , A ( 4  «  3 . 3 ) «  B ( 4  »  3, 3  f  *A  Pi  4 V3 , 3)  *  B  FT*  *3,3)  33 

_ OJMeNSJ.CN  RiV(l),C2Cm),EPISU0),BA<6),BC(6),AA<4,3>»BBM,3)  _  54 

DIMENSION  APII4,3),HPU4, 31, Al(4, 3), 81(4,3). SAM), SBI4)  -  ---  . 55 

_ DIMENSION  EPS1 ( 3 ) , BETA ( 3 ) , OEP ( 3 ) , ALPHA ( 3 ) , PRE S ( 3 ) ,5 PS2 ( 3 )  56 

900  FORMAT ( IX , F8. 1 ,  F5. 2»4E 12.6  )  .  .  ~~  57 

525  FORMAT(6tlP612.5) ) _ 58 _ 

C  INPRO-1 »C AT  A  FOR  HERR INOBONE  BEAR  I N  G  READ  IN  62 

_ NAH6HST/NGPUT/INPRC.NEPS.NCASE.NN.R2P1.  ANTT,DN,RSY,  63 

1EPIS,CUR1,C0R2,PRES/ INPLT/P INSP ,0ELSP,BL0VD,I1, Ii,N9. 

1EPS1, BET  A, OEP, ALPHA i EPS 2  65 

52  NRE-1  ’  '  66 

GAM-0.0.  67 

BOOT -0.0  . .  "  ‘  '  68 

MUST-0 

PHOOT-O.O  69 

_ KOIAG-O _  TO 

NSES1-0  71 

_ NSES2-0  72 

R¥aO 1 5 i N  GPUT )  80  6 

IF (MUST ,NE«1)  GO  TO  53  73 

KOI  AG- 1  74 

__  _  NSESV-1 _ ; _  75 

NSES2-1  76 

53  NCC-i  77 

MR  ITS (6,765 )  *  78  14 

^  505  JF0RMAT(4lH|  FLCATINO*R  ING  MI7H  HERRINGBONE  JOURNAL  ) 

C2C1 ( 1) -0082/ CCR 1*R2R 1  81 

_  MR ITE (6, 210)  82  15 

PCS- MGS  (21  .  83 

210  FORMAT ( 20H  RAIN  PROGRAM  INPUT  )  84 

"  MRITLI6*,N6PUT)  85  16 

IF(lNPRC.EQ.l)  GO  TO  202  86 

”REA0(5, INPUT)  87  20 

_MRITE(6_,212)  88  21 

’  212  FORMAT ( 2Th  HERRINGBONE  BEARING  INPUT  )  ’  89 

HR  lTLtb, INPUT )  90  22 

RPAS-l  ~  91 

_REO-REY(J)  92 

Tall*  me  r  n  b  (  blc  v  cytoRl  ,  cgr  2  ,’Pft  ESV  1 1,  i2,n9',eoot,fhcoT,epsi  , gam,  bet  a,  93 

10EP, ALPHA, RINSP, REC, R2R1.DELSP, MPAS, KCTAG)  94  25 

RlWtND  9  . . .  “  .  . . 95  "26 

00  1213  K-1,3  96 

DO  1215  J-a,3  '97  ‘ 

_ READ (9_>  <P(1,J,K),P(2,J,K),P(3, J,K),P(4,J,K),P(«,J,K),P(6,J,K) )  98  31 


ERROR  MESSAGE  NUMBER  1 _  _ 

1213  CONTINUE  *'  ~  ~  99 

_ REMIND  9  41 

GO  TO  304  106 

-  20J.JLEAO(5J5CO)  ( (P  >1,  J,K )  ,P(  2,  J.K )  ,P(3, J.K)  ,P(4.  J,  R),P(5,J,K»  ,P(6, J,K  101 

1 1  ♦  J*I,3 1  ,K-1 ,3)  102  43 

304  (F (MUST «EQ«1)VRITE16,3GG)((P(1,J,K),P(2,J,K),P(3,J.K) ,P(4, J,K> ,P(5  103 

l,J,KI,P(6,J,K),J-l,3),K-i,3’i  104  59 
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MA I h _ -  6FN  SOURCE  STATEMENT  -  J^FK f S )  - 


_104  _I  K*1  _ _ _ 

00  l  1-1*6 

_ 00  i  J-1,3. _  _ 

DO  1  K*i *  3 

JL_C  <  I  *  J*KI»PII,  J,K» 

2  DO  5  J- 1 > 3 
DO  5  I-3j6 

-  *(W,Jfl)«T.C 
AII-2, J. 31-0. C 
BCI“2.J*n-0'.'C 
B ( 2-2*  J  *3 1-0. C 
AP(I-2* J.l 1-0.0 
APU-^,  J,3»»O.C 
BPti-2%  J,i)«0.C 
8PU-2*  J.3J-O.C 
Y2*Q<I,J,iT 
Y3-QII*  J*2) 

X2»g  ii*  J.l>  "  ' 

X3«U(1, Jj2) 

DO  5  K-2,2 

Yl-Yi 

Y2-Y3 

X1-X2 

X2-X3 

Y3-Q ( I • J i K  +  l  I 

X3*0(if J.K+ll 

0T1-X2-X1 

CT2-X3-X2 

C1-Y2-Y1 

C2-Y3-Y2 

C3-X3-X1 

C4-DT1/CT 2 . 

C24-C2-C4 

A<  Hi,J,KI-tC2A-Cl)/(CTl*C3) 
BI  I*2.J*M-(C24+C1/C4)/C3 
GO  TC  (3. 5). IK 

3  AP(I-2_,.I*M-A<  I-2,J,K> 

B  P  <  I  -  2 ,  J ,  K  » •  B I  I  -  2 .  J  *  K  > 

3  CONTINUE 

GO  TC  17*91* IK 


115 

116 


117 

118 
119 

-  120 
121 

'  122  "" 

_ 123 

124" 

125 

- m  -- 

127 

128  - 
129 

13U 

131 

132 

133 

134 

135 

136 

137 

138 
13«* 

140 

141 

142 

143 

144 


T  IFUNPRC.cQ.I)  GO  TO  208 
MPAS-2 

CALL  HERKB(BLCVD.CCP1,C0R2*PR6S, II* I2,N9>E00T. PHCOT ,EPS2 ,GAM,8ETA, 
1CEP* ALPHA  *RINSP*REU*R2R1* OE ISP*  MPAS  *KC  TAG  V 
REWIND  10 

DO  1220  K-l.l 


145 

146 

147 

148 


149 


UO 


III 

132 


DO  1228  J-l*3  151 

PE  AO  ( 13  J  I  fl  (1 ,  J  ,k » ,  QfJ't  J.K),C(  3*J*K  )  ,0'(  4*J*KJ*  Q  t « *  J  *K  1*0 1 6  ,  J  ,IO  >  . 152  . 137 


i 

V 

% 

V 


ERRCR  MESSAGE  NUMBER  2 
1229  CONTINUE 

REWIND  10 
GO  TO  309 

209  REA0t5|500>HC(l»JiK)tQ(2»U*K)*a(3iJ*K)fQ(  4*  j*'Kl#Q(  5iJ*K)*C(6*J*K) 

.  lfJ»l»3J/Rn*3> 

309  IFIMUST •E0»1IWRITEI6* JOCI  .  . . 


•up- 


HAIR 


EFK  SOURCE 


1  IIU1.J.K).UIZ.J.R1.0(3. 

l.J-1.3) ,K«1,3) 

209  I  M2 

. GC  TC  2 

_  9  00  231  I«-l.hR£ 

RE-REYC It) 

_  00  231  ie»l,NCC 

CC-C2C1UB) 
oo  201_JC»1.,N1PS_ 

s P*’i PIS  (  tC)  . . 

_  ANRT-ANTT 

KD-J 

NA-J 

SO  REi-RE-U.-ANRT) 

RE2-R2Rl«CC*R:l/(l.a/ANPT-l.C J 
C1«P( 1,1.3) 

IF(RlI-CI)  11.10.1j 

10  KC-J 

GO  TC  19 

11  00  13  K-2.3 
Cl-P(i.l.K) 

IFCRE1-C1)  12.13.13 

12  RC-R 

GO  TC  19 

13  CONTINUE 

19  00  18  1-1.9 

co  ia  j-i,3 

AR-APCI.J.KC) 

ei»-opit.j.KC) 

CR-P ( 1+2. J.RC I 
XR-PI1. J.KC) 

RL-RC-i 

AL-APtl.J.KL) 

BL-3PII.J.KL) 

CL-PM+2.J.RL) 

XL-Pfl.  J.RD 
Cl-Rcl-XP 

C1-CR+C1«(BR+C1*AR) 

C2-RE1-XL 

C2«CL*C2*(BL+C2*ALI 

IFIRC-2H5.15.16 

15  C2-C1 

GO  TC  ia 

16  IF  13— RC )  17.17.18 

17  Cl-C* 

18  AA( I » J)  «tCl*C2l/2« 

IFINSES2.NE.il  GO  TO  88 

8T  WRITE(6*325HIAA(I.JItJ-l»3).I-1.4) 

88  00  19  1-1.4  _ 

API! 1.1 )*0.0 

_  AP1II, 31-0.0 _ 

‘BPICI.  31-0.0 
_  BPII I .1 )«0.0_ 

X2-PI2. i.ii 

_ X3_-P  ( 2  »  2.1) . . 

fP (1-4)101. 102. 101 


HA  IK 


SFN  SOURCE  STATEPENT 
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IFMS1  - 


101  Y2*A  A ( I , 1 ) *X2 
Y3« A A ( 1 , 2 ) *X  3 
CO  TC  1J3 

102  V2*AA(1,1) 

Y3-AA(t,3> 

103  OO  i9  J-2 ,i 
Y1*Y  l 
V2«Y3 
Xi*Xi 
X2*Xj 

X3*PI2. J4l.ll 
IF  ( I--4I  104,105,104 
1C4  Y3«AA(I,J41)*X3 
GO  TC  ICC 

105  V3*AAiI,j4i) 

106  CT;«X2'X1 
CT2-X3-X2 
C1-Y2-Y1 
C2»YJ-Yi 
C3-XJ-X1 
C4-0T1/CT2 
C24-C2*C4 

API  C  I  ,  J  1  *  I  C 24“  C I  )/  (  DT  I  *C3  ) 
BP1I  I,JI»IC244C*/C-il/C3 

19  CONTINUE 
C1«P(2,3,1) 

I F ( 2 P-C 1 1  21,20,20 

20  KA-3 

GO  TC  24 

21  OC  23  J-2. 3 
C1-PI2, J, 1 ) 

IF(6P“C1>  •2.23,23 

22  KA»J 

GO  TC  24 


23  CONTINUE 

24  CO  1C9  1*1,4 
AP-APKI  vKAI 
BP-BPK  I.KAJ 
XB»P (2 « KA  » 1 1 
CP«AAii,rA) 

IFU-41  131,132,132 

ill  CR-CMXR 
132  KB-KA-1 

li.-APni.K8T  " 
BL-0P1J  IjKB) 
XL*P(2,Ke,ll 


CL»AA(1,KB> 

IFi 1-411X3. 134.434 


_ JAL£JL”tt*Ui _ 


134  Cl-EP-XR 

Ci-CR*C1  « ( 6R  4C 1 « AR 1 
CX»E  P-XT 

C2-CL4C2  « <  BUC  2  *AU 
ItTkS-2 1  25  tit  ,26 


mfcsi 


GO  TO  2S 


116 

iif 

218 

219 

220 
221 
222 

223 

224 

225 

226 
22T 
228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 
269 


T«X. 


MA  U 


03/04/61 


ECU 


I  r-M  >1 


26  I F  <  3-  KA  I  27,27,28 

27  C1-C2 

...2.8  IFII-4)  J07,1C8,107 

107  SAIt I«(C!aC2)/(2.*cP> 
go  rc  ici 

108  SAIt I • ( Cl*CZ I /2 * 

101  CONTINUE 

RF*SA(l) 

_ $1*SAJJ ) 

RCFNa  (S  A  ll  H$/{  2  )  )/2, 

ATU*SA(4) 

FT* ( SA(  1 l»S  A  ( 2 ) ) / ( 2P*S 1 ) 

Cl*R2Rl*«3 

Ci*CC*CC 

Si»S:*Cl/(C2* (!,♦!./ AN RT  )) 
C1«Q(1,1,3) 

IF ( As2'*Cl I  3i, JO, 33 
33  JC-J 

GO  TC  34 

31  00  33  K-2,3 
Cl-OU.l.KI 

IF  |Rt2“Cl I  32,33,33 

32  JC-K 

GO  TC  34 

33  CONTINUE 

34  CO  36  1*1,4 
00  Jb  3*1,3 
AR*A 1 1 , J, JC ) 

es-aii,  j,  jo 

CR-Q(I«2,J,JC) 

XR-g(i,j,jct 

JL-JC-* 

4L*4(I, J.JLJ 
8L*8 ( I , J , JL I 
XC*0(i , J, JLI 
CL*Q(I*2,J,JL) 

Cl*Rt2'.*R 

C1*Ch,Ca*<8R«C**AR) 

C2*RC2*’Xl 

C2*CLaC2* I BL*C2*AC I 
IFIJC-2)  33,32,36 
33  C2*C1 
GO  TC  39 

36  IF  I 3“JC I  37,37,38 

37  C1*C2 

38  BB(I,JI*(Ca-»C<)/2, 
lFIN8eS2.NE.ll  GO  TO  90 

89  MRirfcI6,323)((BBM,JI,J*l,3),I*l,4) 

90  00  59  J*l,3 
NB*4“ J 
C3»0i2, J,l) 

Ci*C3*C3 
BCINei-.Cl/Cl 
BA(NBI*BS (3, J I 

J?9  CONTINUE 
*1 (3,1)»C. J 


270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 
241 

292 

293 

294 

295 

296 

297 

298 

299 
3  30 

301 

302 
333 
3  04 
305 
396 
307 
■  08 
309 
31) 

311 

312 
213 

314 

315 

316 

317 

318 
219 
i  20 

321 

322 

323 

324 

325 


358 


HA  I K _ -  fPM 

_ All3,3l-C^.y  _  „ 

Bits, 11-0.0 
61(3,31-0.0 
X2-0AU) 

X3-BAI2 ) 

Y2-BC(i> 

_ V3-BC(2) 

00  39  J-2,3 

_u«y* 

Y2-V3 

X1-X2 

X2-X3 

_  X3-BAI  j«:i ) 

Y3-3C  I  JO  1  ) 

DT1-X2-X1 

0T2-X3-X2 

C1-Y2-Y1 

C2-Y3-Y2 

C3-X3-XI 

C4-0T1/CT2 

C24-C2*C4 

*1 (3  » J>« (C24-C1 1/(071*031 
81<3,J)-(C244C1/C4I/C3 
39  CONTINUE 

_ C1-8A(3I 

lF(S2-Ci» '4iYO,4J 
AO  NC-3 

GO  TC  44~~ 

A1  00  A3  J-2,3 
C1-8A( J I 

IF(S2-CH  A2,A3,A3 
*A2  VC-J  * 

GO  TC  AA 
A3  CONtlNUg 
44  AR-Ai(3,AC» 

BR-Bl ( 3  *80 
_ CR-BC ( NO  I 

xf»«e#<Ncy . 

KL-NC-1 
At-A.  (3,80 
0L-Bi(3.8L) 

*  CL-BCI NL  I"" 

_XL«BA(NU _ 

Cl-Si-XR 

C1-C6+.C1*  (  BR*C i  AAR ) 
C2-S2-XL 

C2-CI+C2* ( BL*C2*AL ) 
IF(NC~2V  45,45. 4b 
A3  C2-CI 
GO  TC  48 

46  IF ( 3*-NC  J  47,47, A8 

47  C1-C2 

48  C3«(Cl*C2l/2. 

C2-1./C3 
EP2-.14SCRT(C2J 
DO  49  1-1,4 


cniiare 
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—  i  m  ,  .  i 


_ »i. 

327 

_ 328 

329" 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 
359 

359 

360 

361 

362 

363 

364 

365 
>66 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 
300 
381 


413 


-  SFN  SOURCE  STATEHENT  -  IFMS1  - 
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C1-EP2-*R _ 

Ci«Cfl+Cl*I8R«Cl*AR  » 

C2^=P2-XL _ _ 

C2-C  l*C?T(  8L+C2**L  I 

iMl.C-2)  59.5i.56  _ _ _ _ 

':2»Ci 

<50  TO  58 _  _ _ _ _ 

. i-  !  <*‘LC  |  57,57,58' 

ca»c«  _  _  _ _ 

rFit‘i'1  ii  8 .11 9,1X8 

SBm»jfCl*C2l/(2.*fP2) 

CO  f C  12C  '  -  . - . . 

SBM  )«ICl*.2>/2. _ _  __ _ _ 

CONTINUE 

_ATT2«S8U)_ 

PJ'3. 1415 926536 
BAD-PI/100. 

Ci«EP2*fC  . 

_CZ*160.-mi*ATT2 

IFCC2-90. 1  150,150,151  .  . . 

C2«  10O.-C2 

C2-C2BRAC . 

r_CA— C0SIC2I 

GO  TC  132  . . . 

_C2-C2*R<C  _ 

CA-CCS(C2»  . .  . * 

CSA-SIN<C21 

e pp-s q*  r i cp*e  p* ct*ct-2 .Pci*ep«tii 

CSB-C1*CSA/EPP  . 

ccB-$o»T(i.-csi*cia> . . 

ere-csa/cce  _  _ 

TB-AfANrffai  “  *' 

TP-ATT14RAD-TB 

'YOPtB/P'PC'  *  ~~  . 

EPS-&PP/ ( i.+CC  I 
RCF2-IS6 ( 1 1+SEI2) 1/2, 

FT2-  (SBIU“SJB_<2J_1/(S2*6P2I _ 

"Cl  ■  ( 1 »  -  A  a  R  T~ (Til,  ♦ANRI) . 

IP  UNPRO.K£.l)  Cl-1. 

C2-RCPH*C1 

Cl-Sl*FT*EP/2. 

RCFI-C2-C1 

Ci-sP2*S2*FT2/2. 

Rcro-acpzvci . .  . 

Rin-PCFI/RCFC  _ 

ER-R 10-CC . 

IFINSESl.NE.il  GO  TO  06 

WR  172(6, 523 lREl»Rci,  ANRtiCC  ,‘EP  . 

WRITE (6,523)RCFI,RCF0,ER,S1,S2 

NR  I  Ti.  ( b « 5 251EF 2 , FT » Ft  2 , RCFH ,RtF2' . 

WRITE ( 6 , 325 1  SBtl ) ,SB t 2 1 , SBI 3 > . SB( 4 1, ATT! , ATT2 

IFIN6-1I  70,71.72 

01-ER 

Cl-AKRT* 

ANRT-OUCN 

NA-1 


_y* _ 

439 

460 

441 

442 
‘443 

444 


44  a 

449 

450 

451 

452 

451 

454 

.  435  ' 

45* 

4jjy 

438 

459  " 

_ 

461 

491__ 

462 

. 4*3 . 

494  ' 

464 

49* 

"  443“ 

497 

466 

■'467 . 

49| 

46B 

- - 

”4  *9 

470 

. 471 

472 

-ros 

_4«6 

304 

‘  307 

4B0 

SOI 

490 

492 

493 

Ikt 


-\  ecn  source  statement  -  misi  - 
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"n  If <01*eB»  73,14,74  - 

73  NA«2 

M«AMf  . .  .  •—■ 

03 -EB 

*NRT«(Ql«:J3l/2. 

CO  TC  30 

74  C1*AN«T 

- OL«i<* _ 

4NRT-014CN  - - - — - - — 

. _.hp«NC*i 

IF  <K0~6M  SC, 700, 700 
700.  WRITE  (6,702 )  NN 

T0Z  ClVERciO  AFTER  ,13.6H  TIMES) 

_  GO  TO  201 
72  IF(NA-S)  76,7*, 79  v 
..7*  IFI01*ERI  78,77,77 
77  Cl-ANRT 
01«ER 
NA«3 

„  AW»T«(Ql*Q3)/2. 

GO  TO  SO 
...71  Q3«A6RT 

03-ER . . 

. .  01*031 /2. 

NA-3  ' 

_ ...GO  TO  30 

75  Q2«A6(ir "" 

C2-1R 

RA»4  . 

OTA-fll-Cl 
C1»0TA*C7  A 

C2-(CI*C1-2,*C2I*2./C1 
C3-<C3~C1)/CTA 
IF(C2I  00,81,80 
•1  C4—02/CJ 
GO  TO  82 
10  0TA«.5*C3/C2 

C4»$CRT (CTAAOT  A-02/Ci ) 

1FI0TA)  83,04,84 

83  C4--C4 

84  C4-C4-0TA 
02  ANRT-Q24C4 

GO  TC  SO 

79  S3«S1/I1*4ANRT) 

!F(  PDS)  91,92,91“ 

92  FDi-O.v) 

F02«u,0 
CO  TC  93 

91  POl  ■  POS  /  ( 1 , 4  AE  AT) 

WRITE(6,9O0) 

8RlTE<6,600)EPP,ANftT,CC»R2ftl 
..  «ITE(6,S02) 

UPITt(6,522) PEI ,EF«RCPM,PC1 ,S1, ATT1,FT 


-  EFN  SOURCE  STATEMENT  •»  1FN(S)  - 
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SDBSCUTlNE 'HEANBLBLOVC,  C0Ri,CdR2,PRE5.  11 » I 2*N4 , SCOT , PHOOT ♦  E FSl, G*M 
l.BETA.OEf  .ALPFA.RIN!lP,RFn,  R2#l, nFl_$P- MPAS sROIftC  > 

j^IHU'jNGECNE  JOURNAL  V CORRECT  FOR  TORQUE  "Ik  CAVITATED  REGION, 

1-16-68  _  .  _. _ ... _ 

kTTHLARGc  ECCENTRICITY  ANO  MISALIGNMENT 
IN  TURBULENCE  REGIME 

COHMGN~PCIAGi £ £  LX  »  CEL2 ,  CELZ  S«  M, MS, HG1 , MG2  ,N *NP  1 17 , 1 9  ) *  51 (2,5,1 
19) ,B1( 2,19) *CC (2,19) ,AF1< 17*19) , AF2(17*19 I «AFj I IT  ,19 ) ,  AF4  (17,19)  »  A_ 
*  "'256(17 ,191  ;aF1( J, 7Vi« I.AFJ ( 17, 19), PH  1  fl7 ,19  1 

OIMENSICN  BET A(3),CEF(3l,ALPFA(3)«PFtX(3)«  00Q( 19 » ,M9 (17„ 

T,  1 9  i  ,oqc  c  (19  r;  pp i  a  up i»p  c  it  r,  ppx  ( 17  >  •  ppp  x  1 17 1 ,  a  xm  ? ,  1 9 ) ,  a  x  2  ( i?  ,1$ 

2), 6X3(17,19 1,8X6(17,19) , AX7( 17, 19) , AX 6(17, 19),AX9 (17,19), MS11(17,1 
3  9) ,WS 12Tl 7, 1 9 ) ,NJl3( 1 7 , 19), WS14 (17 • 19  >  t  X  X 1 1 7 ) , A  XL  1 2 , 1 9 1 ,C2(2,19) 
OIMENSICN  £PS1 ( 3 ) , PRES ( 3 ) 

- fling  NSTCk  TtOTYTiT,f9  j  . . 

LOGICAL  FPOUT, VENT, TORO 

“I'raiMAtmR'Tr™-1  . 

4  FORMAT)/ ItMjX.F  11.7))  ..  . 

)  FORMAT  (8 9H  REYNOLDS  NC.  ECCENTRICITY  INNER  7CRQUE  OUTER  TORQUE 
1  SOPMERFELC  NO  ATT.  ANGLE  FLOW  1 

r’PWfrtAT ( 3iH”  fllT  ANG'COiT  FORCE  COMPONENTS  •  S14.7,2H,  R14.7) 

T  FORMAT (J1H  SIN  AND  COS  FOMENTS  ABOUT  2*0  S14.7.2F,  E14.7) 
■nrTflRMATriWoFiVAL  PRESSURE  DISTRIBUTION.  //) 

IF  ( KOI AG.EQ.l 1  WRITE  (6,21  BLDVD, CORl ,C0R2,PRSS(1) , PRESI2) , PRES (9 
rrr  EDBTTp HW, SPflT  1 ) ,  EP SI ( 2» , EPS1( 3  ) , EPS1  (4 1 ,  EPS!  ( 8 1 , GPS1  < 6) 

l,CAM,BETAIl),eETA(2),EETA(S),0EPIll,DEP(2),DEP(3),ALPMAm,ALPHA(2 
~ TTTftFHim  ,ETKSP,Afc"0,R2Rl.  OEL  SP 

_ IND3«K01AG _ _ 

MOi AG»0 


MPASS-FFAS 

_ IF  (MPAS.EQ.il  WRITE16.J) 

“TToAmat  (ImD'  . . . 

_  IF  (MPAS.efl.il  WRITE  (6,1) 

IF  (MPAS. cQ. 2 )  WRITE  (6,1GI 

_ 1  FORMAT  (19X171-  INNER  FILM  DATA  ) 

Id  F0RMAT~riTxi7F“QUTsR'FlLM  DATA  I 

_ WRITE  (6t3i_ 

~GaMDQT  *0.0 

_ MOO  _ 

FGl-U-1 

AMQaPGl 

- Fi?i2-n — . — . 

VENT-.TRUE. 

. IF ( VENT F  CO  tC'i02 

M-I2+I1-1 

BMMl-M-1 

M$2-MG1 

. GO  T C  101 '  * 

102  HmI2  _  _ 

•  ■  bmmi>  m-1 
BNM1“BMF1*2. 

MG2-0 

101  IFIMPAS.EQ.2)  GO  TO  35 

00  2U  1*1, M 


150' 


i 

} 

t 
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_ 00^20  J«1.N 

N>(I,4)«0 
20  CONTINUE 


00  30  J*1,N 
__  NP(1  » Jl"l 
NP(M,J(»i 

_30  _C0NT  INU| _ 

35  KK-1 
OC  AN-N 


_ DTHETA»2.*PJ/AR  _ 

0T HE2*6«S?DTN£T  a 

_  PFIX(l)-PReSll) _  _ 

SRRITO-RtNSP 

PFIX(3)»PRESm 

fr  TPPASS.Tc.if . do  TO'idJ 

_ .80  TQ._J.08 

IQS  SPR IT0«2«*0BLSP+RIN$P 

. . . . 

00  TO  107 
)08  DSLS I"  CELSP 

107  DO  '999  L 6*1  *  . . . 

„  IPIMPASS.EQ.2)  60,  TO  120 
'  *  ASNi-SPPifO-1.0  . . . 

_  anp;«sprito41.o  _ 

.  COR-COR1 .  . 


__  SI8N-ANP1/ANP1  . 
PPIX(2)-PRBS(2»*iWi 


_  AL0V0-6LCV0 

R6N-PeO*ANHl«I-i.I  . 

. _JSLTC..ilO  ... 

120  SIGN--;. 

COR-COR2 

PPIX<2)-PRES(2)*(C0R2/CGA1  »•* 
ALOVDaBLC VQ/ftiRi 

“  REN-«Ea*<R2Rll*+2*COR2/CORl*SPRITO 
130  ELVD-ALCVD 

VAT L  0n2 «  0  * A  L  OVC  . . 

DO  992  L7-1.3 
EPS-EPSHL7) 

_ TPS«2.*PI*SIGR 

ConI-6.*Tp5 .  "  . . 

R API  AN* >01749329291994 3 
SSi«  gSmpr'AO  I AR 

_ E0T*0.0 _ a _ 

CONdaCOKl/REN 

_  TORQ-. FALSE. _ 

’Tp5i*tp571 .0 

_ lRflslPJ.A»Mfl _ ; _ 

DELX-DTHETA 

_ 0LL2-RATLD/BNP1 _ 

t  psi -V2V*  p  n*«2«27*ELvc 
AMS-MS  _ 


I**2/SPRIT0 


97 

98 

99 
100 
101 
102 

"IDS' 

-m- 

106 

■"i'flf 

108 


fftnarr  ft  a  tcucut 


INQ51 _ 

6P02»ls 

LI  «I_" _ 

~  L2-U 

39  00  40  J-l.N  _  _  . . . 

CC"<LW  J)»6.0 

_AXL(Ll,J)-0.  _  _  _ 

cz(Li7Jr-oV 

Al<U<.2,JJ-O.C  _ 

U(U,i7f)-O.C 

40  AF7<L2.J»-0.0  ,  _  _ 

IFUKO.EC.2J3C  TO  48  ' 

JF ( VENT I  GO  TC  48 _ 

Ll-2  " 

L2-12 _ __ _ 

I  NO*  2 

GO  T0_38  _  _ _ 

48  NR-1 

fn-mg; _  _ _ _ 

‘MM- IX 

iss-i  _  . 

“  OEL-CeL’Z'” . 

„oz-c.s/c6iz  _ 

:s-i . 

...  im.o.o  _ _ _ .... 

00  BET>8£T4  4NR>*1T4i:!AN 

06 PH- PEP (NR )  '  _ 

alph-alpVJTnFT . 

ALM1-1.-ILPH 
ALTAI  i-"APH*ifLMl 

SINB-SIMBETJ  _ _  _  _ 

COSB-CaS<3tT) 

SI NB2^$ 1KB*S IK  B  _  _  _ 

COSB2-CCSB*CCSe  ‘ 

COT2-COSB2/SIRB2 

CLZ-CEL*CTHETA 

IF(TCRQ)  GO  TC  100i_  _ _ _ 

00  2J30  !-‘IS. PM 

xxm-z  _  _ _  _ 

ZCOR-’Z/CtR 

EPZGH-6PS*GAl*ZCOR  _  __  _ 

1 0  ZR- EOC  T ♦ G  A  P  C  C  T*Z  COR . 

EPZG-EPZGH*PHIMC_  _ _ _ _  _ 

"OO  2001’ J-X  *K 

Sl-SINI AAG)  _  _  _  _ 

CO-CCSI ARG) 

H«i.4C0*EPZGM 

H9U.JJ-P 

H3-H*H*h  . 

HG-H+OEPh 

HGMR-HG/P 

AX9(l.J)-(ECZR*CO4cPZG*Sn*SINB*Pi<2A,0 
HG3-HG**3  _  .  .  .. 

HGHR3-HG3/H3 

S1-1./SIRB2  _ _ _ 

$2—  COS8/SIN62 


03/04/69 


03/04/69 


-  EFN  SOURCE  ST6TSHENT  -  IFMS)  - 


v _ CSnruTc  Tnc  uirffcKeMCiS 

3100  f ilK 

00  4010  J-l.N 

- im»7T!-c.3 - 

_ AF2(  IjlJI’0.0 _ 

lF3(i»J)*0.0 
_ 45*41 I»Ji*0.0 

»5ir,Jl-0.0  " 

AF6( I 

fft'l  ItJ) *-PFI XI NRI 
4010  CONTINUE 
20fff  lST"lSU 

_ FHWH-i 

00'  4000  l-fSTilNl . 

IO*I-i 

- ypjj-j . . - 

00  4000  J*I,N 

“ - 5fH-0fHll  . . 


IN  XS  AND  CQfF.FtC.IENT 


IFU.EQ.l.OR.J.CQ.M  00  TQ.4CQ4 


OO  tC  4050"' 

4004  IFU.EQ.ll  00  TO  4004 

"jifrN-i . 

- Oti'k  4tf«8 . 

_ 

Jl-1 

4000  AFlt 2.J t*SIN646X?l l.J) 

“  ~~W2(I.ji*tAXati,Jll-6X2(l,jO»4CO$a4tAXTU.Jl)-AX7(I,JOIll*OTH 

_ I  ♦  5IM4(6X7(  17.J>-6X1(  IQ.  J)  >402 

-WJTTTjT-Wt  fi  J 1 ♦COSB*6X7(  1 . J  »  ♦SIN0*AX4(  I ,  JJ 
6F4( I.JI*(AX1( It  41 1— 6X1 1 1. JOI<COSB*( AX6( ! .J1 1-6X6 ( I .JO) I  MOTH 
1  VSlNf  *(  6X6  t (7»  J  t-AXt (  I0.JIX02 

_ AF5( ltJ>*AXl(I.JI.COSe*6X4( I.JI 

AF6(  l.jW.O  ’ 

AFT(  I  »4  l»(AX3(  1. 4H-6X3U. 40  t*C0SB4(AXa<  1.411-6X1(1  .J0t»»*0TH 
I  IT *4  I— AX 0 if 1 0  <  J 1 1 A’dl* 6X9 (  !•  J  I 

IF  (NO, tq. 2  .AND.  J.  EC.  It 

X  .  WRITE  I  A. 2)  AFlt I.J  t . AF2( I.JI.AF3(I.J) »AF4( I. JI.AF5 i 

1I.JI.AFTU.J) 

Awjrwsffwuf" . . 

IFINF.EC.lt  6C  TO  4020 

- WNtfTfC.ir  CC  TC  404C - -  — . 

II  “2 

TE-ia  •  .  . 

_ GO  TO  4090  _  _ 

4020  I  f-1 

_ Ifc-S.F.  . .  . . .  . . 

GO  TO  4990 
4Q40  II-l 

IE-IS  "“  ■  " 

4040  00  4100  4*1  »N 
Aim. i.4i*o.a 
61 ( 1 1 »S • J)*O.C 
'  TM-AXTIIE.JI 


5F/0*  ’  9 


-  EF6  SOURCE  STATEMENT  -  I FN I  Si  - 


_ TSSZ«TSS*0Z*2.U  _ _ _  _ _ 

Al(ii,4,j)— TRSZ 

_ Aim, 2,  J)-AUJI,2,J>-TR$Z*ER02 

D0T«AX6IIE,JI*CTHE2 
_ ccm.ji«czm,j)-ccT 

CZ(  II, JI-CCT 
BUlUJIt-CCIll.JI 
Al(ii»3,J)«Ai(II,3,J)*TASZ 

AF7<  IE. JI-AXL (1 1  »J 1-3X81 te,J) _ 

AXL ( 1 1  ,  J  l-AXE ( IE,  J  ) 

~  AF2  (  IE,  J~)-d,0 

AF3( 1E.JI-0.0 _ 

AF4(  16, jl-J.O 
__  AF5(_IE_.  Ji-Oj.0 
AF6(1E, J)-0.0 

_ IF1MC.E0.2.AND.J.E6.A I 

*  1  WR  IT  €<6~ *51  AKII,?,  j|,Ailil 

_  11  .,CC  ( II ,  J I ,  AF7  (  lEtJt 

4100  CONTINUE 

IF(NR.EC.l)  CC  TO  4102  _  _ 

GO  TC  41 SO 

4102  IS-I1 _  __  _ 

NR«2 

_ K£i.I2 _ .  .  - _ ■  _  _ 

ER02-0. 

OEL-CSLZS 
'C2-J.5/CELZS 
GO  TC__20Q 

4 1 6 0  ~ I  F  I  VE N T  »  GO  T C  4 1 9 0 

_ IF(KK.EC.l)  GC  TO  4lJC_ 

IS- 1 2  . . 

MM-M  _ 

NR«J  ' 

ER02-0. 

‘DEL-CELZ 

CZ-3,S/C6LZ  _ _  _ _ 

GO  TC  200 

4180  IE-12  _  _ 

ER02-1. 

KK-2 
II  *2 

GO  TC_  4090 

Also  IF(NP,£C.2 )  GC  T'0  4l'«C . 

4190  MTSR-MO 
M0  -FTSR 
CALL  CIC1  (MCI 
MC-0 

_ _  PPOUT-.  FALSE. _  _  _  _ 

TF  IK0IAG*NE,1 I  GO  TO  571 
IFUN03«EQ.1.CR.INC3.£0.2)  PPOUT-.JRUE. 
IF(PPOUT)  WRITE (6,9) 

571  00  5  75  1-7  ?  M__ 

|F( PPOUT I  WR I f  E ( 6,  4 I  ( P  H  I  ( 1 » J  I , J -1 , N I 

_ &P-.S.7.S-JJ.1.1.N _ _ _ 

KUPT ( I , J 1-0 


280 

281 

282 

~2B3 

284 

'285 

286 

287  - 

288 


,ji,Armv3Vj),Biui,j 


248 

280 

249 

300 

301 

302 

303 

304 

303 

306 

307  ' 

308 

309 

310 

311 

312 

213 

314 

315 

316 

317 

318 

314 

320 

321 

322 

323 

324 

325 

326 

319 

327 

328 

324 

330 

331 

327 

332 

333 

332 

334 

335 

153 
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IFCPMI  I.JI.Ge.O.O)  GC  TO  37!  ..  . . 

PHI  (  1 ,  J  1*0.0 
KUPT 1 1 , J 1*1 
575  CONTINUE 

C  DIMENSICKLiiSS  FLOW 
IP- 1 1-1 

4204  AFLOh-J.C 

00  4200  J-l  ,N 

IFU.EQ.i.CR.J.tQ.M  GO  TO  4220 
VFl-IPHl  (IP,  JM  l-PHI  1P.J-1 J  )*CTHE2 

4210  AFL!JN-( VF 1*4X6 1 1P,JI|-»AX8(IP.J l+AXTI IP. J I* (PHI (  IP*1,JI“PHI 1 1  P-1, J) 
1 1  *DZ-*  AFLCW 

IF(  MCI  AG. EC. 2  I  WMTEIb.Z  lYFi.AFLOW 
'  GO  TO  4230 

4220  IFIJ.EQ.M  GC  TO  423C 

VFl-IPHl<IP,2l-PHIt IP.Nl I40TF62 
GO  TC  4210 

4230  VF1- (PHI ( I P» 1 l-PHl I  IP* N- 1 1 1 4CTHE2 
GO  TC  4210 
4200  CONTINUE 

AFLQk- AFLOW*DIF IT  A/12 • 

C  TORQUE  CIVICEC  BY  MUXNXRXRXRXR/C 

MR-1 

T0RQ-.TR1.2. 

TC-3.0 
IT1*1 
TT2-U 
TCi-C.O 
GO  TC  23C 

1001  CO  1000  I-IT1.IT2 
00  iuOJ  J-1,N 

IFIKUPTU,  Jl.tC.lI  GC  TC  li-CC 
IFII.eQ.ITl.CR.I.ca.IT2|  GO  TO  1S0C 
0P02-IPHK  Ml.JI-PFK  I'l.JIl/CZ 
AFTR-1.0 

1200  IFIJ.EQ.1.0R.J.2Q.NI  C-G  TO  ItOO 

CPOT-IPHII  1,J«1  l-PHM  l.J-ll  M0THE2 
1240  H«M9<|,J) 

BQ-H*ALT  ILi*C£PH 

AC«(-BCPC0SB*hS13( I ,J I-FPALPI-WS12 I I,J I-AIPH9MS 11( I , J l*IM+CEPHII • 
10P0T 

A0*AC*BC»(-CCN2*MS14( I ,J I-NS13 (1, J I*DPD1* SINBI 

TQ-TC*AQ*OL2*C.5*AFTR 

GO  TC  10U 

1600  IFIJ.EQ.M  GC  TO  1610 

0P0T-(PHII»2I*'PHI(  I,NM40TFE2 
GO  TC  1240 

1613  OPOT«(PHIl,l)-PHI(  l,N-l)»*CIHt2 
GO  TC  1240 
1803  AFTR-,5 

IFII.cQ.IT2)  GC  TO  16 10 

0P04-CPF II IT1+1 , JI-PHI 1T1 ,  J )  I  /CELT 

GO  TC  1200 

1810  OPOZ-IPhtl IT2,J)-PHI( IT2-1, J) l/OELZ 
GO  TC  1200 
1010  RE-R£N*H 


336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

330 

331 

332 

333 

334 
333 
336 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 
38B 

389 

390 

391 


364 
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392 

432 

A£-R£N*|H-OEPM 

34i 

TC1*TC1+; dC2*T CC TRE  )*AL PH  )*TRQ*DLZ*AFTR 

394 

433 

IP  M-F Eil!  EG  TO  10Q£ 

3w 

IF(MC.EC.2.ANC.I.EU.IT1> 

394 

1  WRITEI6*2 )AO*TO»RE,TCl»TC2 

397 

43r 

1COO  CONTINUE 

398 

1F( MC • CC,2) WR I TE ( 6,2 )WS 11 ( 1 » JltWSl 2( I »  J  ),WS13< I,JI ,MS14tt,JI 

"349 

442” 

IF(NR.EC.l)  GC  TO  1400 

400 

IF ( NR,£C .2 )  GC  TO  1410 

Tfll 

GO  TC  578 

402 

1400  KR*2 

403 

iri-ii 

404 

IT2* 12 

403 

GO  TC  203 

406 

1410  I F ( V  E  NT  1  GO  TC  578 

40T 

MR*3 

408 

ITl* 12 

409 

IT2*F 

410 

GO  TC  20G 

~  411 

578  TOO*  (TO'tSTCl  J 

412 

TH?*U»0 

413 

00  380  J*1»N 

414 

000( JI-SINITHE I 

41S 

468 

COQOI J)-COS(TFE) 

416 

470 

580  THE-CTHcTA+THE 

417 

00  597  1*1, F 

418 

FPd  l-'J.C 

419 

pppu  )«!.:» 

420 

00  6U0  J *  1 » N 

421 

CUM*  PH I ( I , J | 

422 

PP<1  )«PPd  >*CCC(J>*CUF 

423 

600  PPPtd-FFPdMCCQQUMCCM 

424 

PP( I )*FP( 1 1 *CTF wT A 

425 

PPP(  IJ-FFP1  IMCTHETA 

426 

ppx(  n*pp<i»4xxui 

427 

590  pppxii »*FPpm*xxi  n 

428 

FS1N»SUP(PP,M,CELZ ) 

429 

501 

FCaS*SUF(PPP,F,DELZ) 

430 

502 

fhsik-Sup(ppx,p,oelz) 

w  - 

303 

FMCOS-SUPI PPPX,H,OELZ 1 

432 

504 

FCOS—FCCS 

FNCOS— FPCOS 

1 T  (  Xl5 lT57Nf .  rr  GO  TC  "i94 

43$ 

MRIT£(6,6)FSIF, FCOS 

434 

507 

WRITE (6, 7 IFMS IN, FNCOS 

433 

308 

594  M.0AC-FCCS**2*FSIN**2 

436 

MLOAC-AaSlNLCATl 

437 

MLOAO*SORT(MLCAO) 

438 

510 

SOHER*2.0O*ELVC/WLCAC 

PHtE*ATAN2 ( FS  IN, FCOS ) 

440 

311 

phIe-phHTracian 

441 

TOO* T 00/ h LOAD 

445 

TQI*TQ04EPS*FS IN/HLOAC 

446 

WRITE  (6,2)  REN,EPS«TGI , TOO, S ONER, PHEE,AFLOW 

447 

312 

lFIMPAS.EO.il  GO  TO  593 

448 

137 
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SOURCE  STATEMENT  -  1FMS)  - 


CICI 


EENSOliHCE  STATEMENT 


IFfcISl 


SUB A CUT IAE  Cl  Cl (MO) 

CPMMCN  FClAC.CEi.X.CSLZ,egi.ZS.K.WS,HGl,HG2,N,NFIlT,19) >11(2.5.1 
14),Bl(i»l,*r.Ci(2.19).AFl<17.1'9V,iF2(lTf19»,AFi(lT,l9r,A^(lT,l9r,T' 


2F4(17,19) , AF7 (17,19), AF5I  1T>  19  >  ,PH1(  17,19)  _ 

01MENS1CA  A(17, 17) ,B( 17,17) ,C ( 17. 17 ) , C ( 1 7 , 1 7 1 , F(20 .17 ) , A(f(  1 7,17 ) ,C 
1NI17,17),S(2Q » 1 71.  AFJIJ 17), BC( 17.17), BF(17),0D( 11,17) ,PS( 17)  & C 17, 
217), GK17, 17), £(17,17) 

AW41.Q/CEIX _ 

NC-0 

M*MS+M61*MC24l _ 

MN1-MG141 

MN2-MM+HS _ _ _ _ 

0X1*AN*0,5 
0X2*  AN««2 


209  00'  205  1*1  »M 


2 

3 

- 5 - 

_ 5 _ 

6 

_T _ 

3 

To 

11 

12 

_13__ 

14 

15 

"  t* - 


_ CO  2u4_U*ijM _  _  __  _  _  _  17 

ETi,ii)*o.'  .  .  "ft . . 

204  01 1,  ll)«C. _ _ ; _ _  _ _ _ 19 

FlTTn^C.  ~  *  """  ~  20" 

205  0(1,11*1.  __  _  _  __  21 

GO  300  j*i,lT  *  . .  .  .  21 

WHITE (7 ) ((E(I,11),0( 1. II), 1*1, H), 11*1, M)  23  23 

IFTSoiTSm  GO  TO"  240  2* 

_ WH 1TE  (4,1031  AF  IJ  i  ,  J ) ,  A  F2J  1 ,  J ) ,  AF  2 J1 ,  J  ),AF4ll,J  I.AF5  (1,J  )*AFA  (_1 ,  J  _ 25  _ 

1 1,  AF7(  1,  J  I  "  *' . .  24  34 

WHITE  (4,103)  27 

1  TiTITiiJ),Ai»l,4VJ),Al(l,3,J),4i(i*,4.j),Al(i,5,J),Bl(MNl,  2B 

1J)  ,CC  (MM,  J )  29  42 

WHITE  (6,103)  * .  .  "  30  . 

1  A1(2,1,J),A1(  2,2,4  ),A1(2|3,<J)  ,  A1  (2  ,4,  J I  «A1  (2  ,S,  J) , _  31 

mTMN2“i  jTTCCT ,  j  f  """  "*■  '  ~  32 . 50 

240  0Z1*J.5/(GU.Z)  33 

'022«1.7TC¥lT)«*Z  34 

241  00242  I-ljM  35 

00  242  lI*i,F'  34 

B(1.I1)«S. _ _  _  37 

Alf,  1 1)*^  ~  '  . .  38 

242  CU.IJLL'C*  _  _  _  _  39 

DO  250  r*Tv>  "  40 

206  IF ( I*MG1“ 1 1  212,210,212  41 

20  9  AF  J I  (Tl  AF71  f  ,’jT  "  42 

A(I,I)-1.  43 

Gfi  tt  HiC  '  . .  . .  "44 

210  IF(MS)  211,212,211 _ _  _ _  45 

Til  ozi*o,s7celzE  '  .  .  44 

0Z2*1./(CELZSI*«2  47 

IT(  N  P"(T,TTf21i9  iljl  ,TB9  .  "  48 

232  IF(MCl)  233,212,233 _  „ . 49 _ 

233  IF(NC,EC«2)  GC  TO  600  50 

All,  I“2 1  •  A1  ( ljl_,  J )  _  _  _ _  51 

A(I,r-il*'AiTl,2,j1  '  '  . 52 

A(I,I)*A1U,3,J)  _  _  _ _ _  _ 53 

A ( I ,  !♦!  )*  Al  1 1 ,4  ,.J  I  -  —  -  .  -  TJ4" . . 

_ AJ  I  (  I'A’jl  > *  *4J  lj 5 » J )  _  _ _ _  55 _ 

6(1,11*81(1. J)  54 
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CICl.  -  f.FX  SOURCE  STATEMENT  -  IFK<  5 1  - 

C41* i*-CC«i* J*  57 

600  AFJUIJ— AF7(  l.J)  58 

GO  TC  25C  59 

212  IF(NPI1»JII2U9,234*2C9  60 

234  If (l'MG;-F$-i>  215,213, II*  61 

213  IF ( MS )  214,215,214  62 

214  CZ.-'J.S/CGLZ  63 

0Z2- ( 2 « /C  £(.Z )**2  64 

IF  (NC.tl2.at  GC  TO  jJi  65 

4(1, 1-21*11(2, i,Ji  66 

4(  I, I— 11*42(2, 2, J)  67 

~  4(1,0-41(2, 3,  J>  68 

411, 1411*44(2, 4, J»  69 

A(1,142)*A1(2,5,J)  70 

B<l.n-81U,J)  71 

C(  1, 1  l-CC (2,  J1  72 

6C1  AFJtdl  —  4F7(I,J»  73 

GO  TC  25C  74 

215  IF(NC,tC«2)  Gt  TO  «k»2  75 

B(I, 1— il*4F3(l,J )* CXI* Oil  76 

B(t,ll«AF5(I,4)*CXi»Af»(I,J)*0Xl  77 

8(I,-|4l|-“AF3(I,JI*DXl*CZl  78 

4(I,J-H*4F1(  l,Jt*C22-4F2C  I«J)«C21  79 

4(l,n*4F&(I,J)-  2.  j*(IFm,J)«C22  *AF£( 1 ,  J  1-0X2  I  80 

411,  14U-4FK  I,J|*C22«4F2(  I,J)*DZ1  81 

C(I,i-n  —  AF3Cl,J)*DXi*CZl  82 

CU,I*-AF4(l,J»*DX.4AF5(ItJI*0X2  83 

C(I,l4ll»4F2M,Jl*CX**CZi  84 

602  4FJKII  —  4F7-II.JI  85 

25C  CONTINUE  86 

IFINC.CC. 2)  CC  TO  i j 3  97 

CO  itO  1-1, K  88 

00  it.<  1 1- i  ,M  89 

4K(I,U)*A(  1,111  90 

00  260  111-1,8  91 

260  AKt!,in*AK(i,in«im,iH)*6<iii,in  92 

|F(K(:I«3.1T.2>  GO  TO  it!  93 

MP1T£(5,:jUJ  94  138 

OO  ifc 1  1-2,8  95 

261  MF  IT  &  (6 « (AMI, 111*11-1, 81  96  191 

262  CALL  H4T1NV(4K,M,CUM,C,CUM11  97  198 

613  CONTINUt  98 

CU  424  1-1,8  99 

t)F(l  f-AFJ!  1 1)  100 

F(J4i,l)*0.  101 

00  404  11-1,8  102 

r(l,llt-c.  103 

80(1, 111-0.  104 

CU  4o3  111-1,8  105 

€U,II»-f  ( 1.111*  AM  l,lin*C(  Ill,  III  196 

433  aDii,iii«RC(i,ii)4e(i,iin*c(iu,in  iot 

404  BF(n-8Fm-e<i, m*f( j,iii  10s 

00  4 06  1-1,8  109 

CO  4C6  1 1-2,8  110 

ou.m-3.  m 

CO  4J5  111-1,8  112 


SOURCE  STATEMENT 


IFMS) 


4fl»  od»in*cd.in-A}<(  i,nii»ecmi,in  _ 

60s  FIJU.d-d.i-ii.Ii+ANti,  ii  i*SH  1 1) 

300  CONTINUE _ _ _ 

BO  503  1*1 »* 

__ _ CO  SO»  1 1-1.  P _ _ _ 

HoOpii  ,  1 1  >»— di  I ,  i  n 

503  CPU, 11*1. *0011,11 _ _ 

IRIN01AG.lt. 21  GO  TO  264 

_ _WR  ITE _ ( 6. IPX  I _ _ _ 

DO'  263  l*I,M  - 

263  ijlTi  J 6 . 100 1  (  C0(  I >11  >  d  1*1 , NJ 
"264  'CALL  MT INV  fCCYP,  OUM.O  ,  CUMl > 

DO  507  I*1,H _ _ _ 

$<N,I)*Q. 

00  307  II-i,M 

GNI  Mlf-O.  '  "  . 

00  306  Ilj*l,F 

GN 1 1 , 1 1 I*  GN  ft ,  II  Rfffi (1 VI 1 1  t*E(  111,  HI 
306  G(l,II)»GNU,m 
5CT  S(N,T 1*$  IN*  I  MCD(  lV|  I  J*FiN*i»"Il  I 
J»N*1 

DO  512  K*2,N 

HR  1  TelB  I  KGd,lI),I*l,PI»II*l«MI 

. BACKSPACE  7 . 

R£A0|T1_  1 1  Ed,  II  I, Oil,  II  1,1*1, Ml,  1 1*1, Ml 
"  BACKSPACE  7 

IFIMCIAG.EQ.2)  HR1TEI6.100I  (IE  1 1, 1 1 1.Ct I , I II , 1*1, N| ,11-1, P| 
J*J-i 

DC  309  I«l,K 
SlJ-i,U*F(J,I» 

00  509  11*1, M 

Gin.TilJo.  ~  "  . . 

CO  508  111*1,1 

Gi(i7in-irnT,nv6ciT;rrii*GKiiii,iii*5rriniT*Bi!n,n) 

IF  I  MCI  AG_*NE*  2 J _ GO  TO  SOE 

308  CONTINUE 

308  SI J-l , I)*S(J-1, I )♦ c I  I , I II  *  S  U ..111*0(1, IJ I^S  I Mil  n _ 

00'5T2  1*1, M 

CO  5*2  II-l.P 

Ell'  g<  i,  it  i*g  tiivii  r . 

280  DC  511  1*1, P 
CO  310  II-UM 

310  OOII.HI— GM.YIII  _ _ 

line'll,  tT*lY»C£l  iTTT  ' 

_IF|M01AG.lT„2I  GO  TO  266  _ 

HRiTE'urrpiv .  ■* 

00  263  1*1, K 

imnrm-'nf.iMT  rcai  r,Tir,ii-i,M i -  -  - - 

166  CALL  MAT1NVI CC,H,QUN,C,CUN1 1  _ 

iPIMClAG.LT.2l  GO  TO  27C 


1*0 _ _  _ 

WRITE  (6,10111 
00  273  1*1, M 

12TTflnfT'T671MTTcBn ,  f  nTTTi'lTri 


I  4111,11*0 


Mar  «Tf  v 


-  EFN  SOURCE  STATEMENT  -  lFN(S)  - 


_ MLIttJIXsi*" _ _ _  . 

sis  PHm,i)«Dou»ii  i*su»m+pn<  i,i> 

_ 00  Si*  Ja2,  N 

~  SACK  SPACE  8 

_  READIB1  ((Gtt,lI>,Ial,M>,tI«l,M J 

Sack space  a 

DO  Si*  Ial,M 

DO  SX6  llal,M _ _ 

si*  phi(i,j)*phi(  i,ji-»G<i,in»PFi(  iuYT 
REWINO  7 
REWIND  8 

_ GC  TC  no  _ 

1FIMCIAG.LT.2)  GO  TO  266 
110_WRJTE  (6,1021 

00  2*f~Tal,M 

2A7  WRITE  (6 tiUO)  (PjHK  I,J  l.J-l.N  ) 

268  RETURN  " 

100  FORMAT (1* « 1 PIC £11. 4  ) 

101  FORMAT ( 5H~ 1C I C 1* I  S  I 

102  FORMAT! i>0( 30K10HF INAL  PHI  MHO) 
ici  format  c fx\  iPTtii.M i 


9 

0 

1 

2 

441 

3 

442 

4 

450 

njnn  non  non  oon 


MATK. 


6FN  SOURCE  STATEMENT  -  I FK C S> 


03/04/6* 


10 

15 


20 

30 


SUBRCUT IKE  MATINVI A»N,B.M,0ET5R I 

-ATP  IX  JKV5RS ICN  WIT}-  ACCOMPANYING  SOLUTION  UK  LINEAR  EQUATIONS 
0IMEMSI CN  IP  IVC( 171, A 1 17,171.8(17, l  >,!NCEX(17,2  I.PlVOTIlfl 
EOUl VALEKCE  I  IRON, JR0H I •  (ICCLU  »JCOLU  I,  UMAX,  T,  SNAJM _ 

initial izaticn 

ODTER  -1.0 

OP  2m  J-l.N _ _ 

IPIVC  (jl-0 
00  550  1-1  ,N 

SEARCH  FCR  PIVCT  ELEMENT  . . 


AO  AMAX-O.C 
AS  00  lu5  J-l.N 

SO  IF  IJPIVC  Ul-ll  60»  ION  60 
60  CO  100  Ml  (ft 

70  IF  (IPIVC  ( K  I“1  I  60*  100,  7 AO 

80  IF  (ABS  (AMAXI-ABS  (A(J.KII)  65*  ICO,  100 

85  I  ROW"  J 

90  ICOLU  * K 

95  AHAX-A(JiK) 

100  CONTINUE 
105  CGNTINUc 

110  IPIVC  I ICOLU  I-IPIVC  I ICOLU  )+l 


130 

1A0 

ISO 

160 

170 

209 

205 


INTERCHANGE  ROMS  TO  PUT  PIVOT  ELEMENT  ON  DIAGONAL 
IF  (IROW-ICClU  I  lAJ*  260.  140 

OETcR  —CETER  _ 

DO  2u0  L*1,K 
SWAP-A(IWOW.L) 

A(  IRCW.LI-AI ICCLU 
A I  ICCLU  .LI-SMP 
IF (HI  260.  260.  213 
DO  25U  L»i.  M 


.LI 


4 

. 5 

'  6 

7 

8 
9 

nre~ 

11 

12 

13 

14  ~ 

15 

16 

17 

18 
1* 
20 
21 
22 
23 
2A“ 

25 

26 
27 
25 
29 

50 

51 

32  " 
S3 
34" 
S5 

36 

37 


220 

230 

250 

263 

SWAP-BI I *OW»L 1 

B( tRCU.L  )*B( ICCLU  .LI 

B(  ICCLU  .LI-SNAP 

INDEX ( I . 1 1* IRCW 

-  35 

39 

- - '  . . 40 

41 

270 

310 

INDEX ( I . 2 1* ICCLU 

PIVOT! I)*A( ICCLU  .ICOLU  ) 

. .  42 

43 

320 

OETcR  -CETER  * PlVOT III 

44 

45 

DIVIDE  PIVOf^RCH  5v~PiVot  fL'IWfliT 

" . -  -  4* 

47 

330 

340 

A( ICCLU  .ICOLU  1-1.0 

DO  350  L-l.N 

‘  48 

49 

350 

355 

A(  ICCLU  .U-AIICOLU  .Ll/PIVOTII) 

IF (Ml  380.  360.  360  . 

30 

81 

360 

370 

00  370  L-l.M 

61  ICCLU  f  L I  —  6  ( ICOLU  .ll/PIVOTUI 

82 

S3 

_ RgOUCg  NCN-PIVCT  ROWS 


35 


56 


MATH. _ "  EFN .  SOURCE  STATEMENT  -  IfkISI  . "T  "  * 


310  DO  590  If  1»K . 5? 

J90  I MLl-ICCLU  I  400,  950.  400  58 

400  TjUl  llj  ICOLU  )  59 

“iJV'AIUtlCClU  1-0.0 .  ~  60 

430  00  490  L«1jN  61 

"W  *  (U  .1 » -At LI ,U-A»  ICCLO  »  L  )*T  62 

499  1F(H>  55C »  590,  46b  63 

440  00  500  L»1 > A  64 

_  900  8(tl,l)*e(Ufll«BUCGLU  ,U-T  65 

590  COKtfNUE .  66 

C _  _  .  _  67 

C  INTERCHANGE  "CClUKr.3  60 

£ _ _ _ : _ ; _ , . .  69 

400  00  710  I«l,N  70 

_  610  L*N*1-I  71 

'62'0 "Tf  ( INffEM  L»1  t-~ INDEX (L. 2 ) )  620,  71C,  630  .  72 

_  430  JNON-lNCEXtUll  73 

440  jftOLU  ■  INOiXU.il  ■  74. 

630  DO  703  _M1|N  -  79 

460  SWAI»«A<K,JftGW)  7 1 

670  A(*,JBON>-A(K,JCOLU  I  ,7 

100  MK,  JCQLU  (■SNAP  78 

709  CONTINUE  '  79 

“flO  CONTINUE  80 

740  AEfURN  91 

T  "END  '  .  .  82 
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gtcf.  -  pffc  so<j°es  ST*T£“Er;?  -  :rnjs;  - 


FUNCTION  GTCF(REYN»  2 

RE-REYN  _  _ 3 _ 

IF (R£*LE •70*0 ) GG  TO  2C  4 

IFMRE.GT.70.C).ANC.(RE.LE.2COO.O)J  GO  TO  30  _  S 

10  IF<(RE.GT.200C.OI.AND.<RE.LE.S.SE*03>lGd  TO  40  * 

IF((«E.GT.550C.O).flNC.(RE.LE.2.0e-l04M  GO  TO  50  _  T 

GTCF-2J •  5/ ( RE  1**0 *784  8  “  14 

RETURN  _  9 

20  GTCF-1. 0/12.0  10  “ 

RETURN  11 

30  GTCF«.619E-08*(RE)**2-3.4«5E-0S*RE*.0856S  12 

RETURN  _  13 

40  GTCF-4.9C/(RE)**.628  14 . 18 

RETURN  15 

50  GTCF»10.35/<RE }**.716  16  20 

RETURN  17 

END  18 
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SZCF.  -  Egg!  SOURCE  SUTEWCNT  -  IgjMSj 
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funStTon  czcfTHynT 

"C-M™ 


!F<RS.tE.7i). 0)80  TO  20 


If ((RE.CT. TO. C)»4NC<(RE«LE. 4000.0) )C0  TO  30 
10  lFI(RE.6T.4000.0n*NC.IRE.LE.7.0E403>>  CO  TO  4C 
lFMBE.0T.7OOC.O).4NC.lRE.Lt.7.0E-»Q4))  CO  TO  80 


MC#»E5 

RETURN 

iO  GZCF-1.1 

RETURN 
30  GZtF-l.i 


RE>»*.tS6 


£-C9*UJEI«*2-i.a7«E-d5*RE^.a346" 


T 

l 


T 


fffll 


6FN  SOURCE  STATEMENT 


SUM. 


FUNCTION  SUM ( F i H,DX 1 

_ DIMENSION  PU7) _ 

K-2 

KK-H-1 _ _ 

MKK-2 

SUM-0.0 _ 

10  00  20  I-K.KK.KKK 
20  SUM- SUM4F  < II 

60  ft  (3?(>O.90).K 
30  SUM- SUM4CX/ 3. C 
RETURN 

40  K«3 _ _ 

45  5UM-SUM42.0 

go  rc  io_ 

30  K«t  . 

KR-M 
KKK-M-l 
CO  TC  45 
ENO 


ICC.  -  EFN  SOURCE  STATEMENT  -  IfMS) 


fjMCTiUK  rCCtKtt 

_ IF(RE.GT.lOO.C»  GO  TO  10 

- IfCC-t./M 

GO  TC  100 

~WTM*fe.GT.400.C>  GO  TO  20 

_ TCC«4#  175/ 1  (*£>**•  86 

GO  TC  IOC 

_20 _ IFIAc.GT.1000.0 )  GO  TO  30 

tC0.547/Ue  l**.521 
GO  TC  IOC 

30  IFIRE. GT.4J00.G)  GO  TC  40 

_ TCQ* • 342 / ( R£J  4**453 

GO  TC  IOC 

40  TCC-.064/IREI**.2S 
100  tf TURN 


MCJTljtg-MNO  KIT*-  t-EMINCB0U6  jCkllNM. 
HU  M0G*»*  f*?uT' 


t'wi’V  1 

CNMC 

• 

0.  "  *  '  ■ 

NEM 

■ 

"  » 

NC»»E 

* 

1, 

NN 

• 

10. 

*2*1 

■1 

o.iioooocoi  oi» 

ANTI 

* 

c.2»dooocot-oo. 

ON 

C*69949999f-qo» 

ACT 

■ 

C. 099949991  02* 

EPIS 

• 

c . iooodooot-co .  o. ioc  ooooos- co. 

-C. OOOSOOOOE-19 ,  -0. OOCOOOOOE-19, 

0.500000006  00, 

- O.OOOOOOOOE-19, 

COM 

m 

8.009990991-01. 

_ £P8I _ 

u 

C *999994991-03 • 

_ ?SfL._* _ Ei _ _  .  _  .  Oi«tmiME  ci.  -a. 


HilRl^lONl  BiARiMi  INPUT 

IIMUT 

- 

-  . 

AINSP  • 

c.nojuooot-oi. 

ofiir  • 

C. 099999991-00 • 

•tCw  • 

C. 0999999 91  01. 

ci 

3* 

12  « 

7. 

N»  •* 

"  "  . .  12." 

mr  • 

(.200000007-00. 

ll!i  « 

0.109000006  03. 

p"  ■ 

0*  • 

lt»H*  • 

-c.  » 

Mil  r 

0.300000001  do. 

TBc 


o.iocoooooe-co. 

“ff.lOTWii  ci, 

o.limodiE  il," 

o.liidooooe  to. ~ 
t,']ntnlK  "to," 


INNER  P ILM  CATA 

REYNCLOS  NO.  ECCENTRICITY  INNER  TOCOUE  OUTER  TOROU5  SOHNERFELD  NO  ATT.  ANGLE  FLOW 

3. 550C00u=  01  C.2OUC00DE-O0-9. 1273569=  0 2-0. 1 2 8 5 3 75E  C2  0.2129888E  01  0.3618759E  02-0 .72 1  8  548  E  01_ 
0. 5500000=  01  C • 300CC  OOE-OC-O .8305493E  0 1-0. 848 3987E  Cl  7.13699355  01  0.36511BIE  02-0.756541CE  01 

3.5500003=  01  0. 5O0CC00E  00- 0 . 46 7 8 2 5 8E  C  1-0. *98 5645E  Cl  3.7334086E  CO  0.3793541E  02-0. 8 67991 3 E  01 

0. 6  5QC0  00E  01  C.200CC005-OC-0.1531C37E  C  2-0.  1  54  8261 E  02  0.1993225E  01  0.3424303E  02-0. 6979 1 8eE  01 
0.65CCCO0E  01  C.jOOCOOOE-OO— C.1C01766E  0 2-0 . 10 1 875 1 E  C2  3.12793455  01  0.3448471E  02-0 . 7299 742E  01 

3. 65CC003E  01  C.5jOCCjOE  30-0 . 56 354 32E  0 1-0. 5920835E  Cl  0.6808451=  CO  0.3564804E  02-0. 83 30470 E  01 

0.  7500000E  01  C.20CCCU0E  •00-0.1795C31E  0  2-0 . 1 80  56  62E  C2  3.184742SE  01  0 . 321 11  24E_02^0 . 6772657_E_01 
0.750C003E  01  0  ■  300  C  v.  OOc  -00—0 .116E741E  02-O.U64764E  C2  3.118325CE  Cl  0.3228307E  02-0. 7066473  E  01 

0. 750C003E  01  C.530C000E  00- G .656 5S 3 5E  0 1-0 . 6 8 39612E  C!  0.6259595E  00  0.3322680E  02-0 . 80 1 2 1 16 E  01 

CUTER  FILN  CATA 

REYNOLDS  NO.  ECCENTRICITY  INNER  TORQUE  OUTER  TORQUE  SONNERFELD  NC  ATT.  ANGLE  FLOW 

0. 36CCOOO;  01  C.5QOCCOOE  OC-C .60 3 45 2 o£  0 1-0. 611 6465E  Cl  0 . 34 803 82 E-00  0.93856C4E  01-0 . 2 802687 i  02 

0.  3600703=  01  C.bJOOCuOE  j > 3 . 5C2 44 80=  C 1-0 . 51 2 4692E  Cl  7. 2 722008 E- 00  0.135B461E  02-0 .3048 566 E  02 

0. 360COOJ=  01  O.SOOCGOOc  00-0.4138659=  0  1-0.  4  34  9845E  C!  3 . 1  781  508  E-00  0 . 1  53  3646E'~02-6. 3678925T_02_“ 
0.5040000=  01  0 . 50CC=G0t  Ji- 7 . 71 3 5 * 47E  0  1-0.  72 5 7334E  Cl  0. 4060655E-C0  0.1280149E  02-0 . 20 57270 E  02 

0. 5C4C0035  01  C  .600:C*70E  00- 0.5860299=  0  1-0. 6  02  5  070E  Cl  0 . 3147292  E-03  0.1396234E  02-0.22319846  02 

0. 5C4C003E  01  C.dUJCCj3£  30-u .470 37 3 1 E  0 1-0. 4960753E  Cl  0 . 2008 538 E-00  0.1874024E  02-0.2681073E  02 

0. 64  8C  700=  01  0.500CC30E  OC- 3 . 79 1 1 4 38C  0 1-0. 804 2 8 33E  C!  0. 446 54 74 E- 00  0.1523563E  02-0 . 164 3149E  02 

0. 6480303=  01  C.oOOCCOOE  30- J.6466  842E  C  1-0 . 6  63  7  5  70E  Cl  3.343929CE-00  0.  1633265E  02-0. 17J8327E  02_ 

0.6480300=  01  C  •  6  jOC  COOc  .70-  ■) .  5074  26  3E  C  1-0. 536  17C2E  Cl  3 . 21  572  08  E-00  0.2104944E  02-0.21  267 11 E  02 

OUTPUT 


ECCEKTR/C1 

N2/M 

C2/C1  P2/RI 

0.60C1E  OC 

0. 39416- 

00 

0.12C0E  01  0.12C0E 

01 

RE  YNCLCS 

NO 

tCCENTRIC  TORQUE 

SUPPLY  PRE 

s 

SONNFD  NO 

ATT.  ANGLE 

tang,  force 

INNER  RING 

0*6)59g 

01 

J.20CUE-CC  -0.1427E 

02 

C.3C126 

01 

0.2C55E  01 

0.3512E 

02 

0.2800E-00  ^ 

OUTER  PING 
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